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ABSTRACT 


Theories  of  mechanical  behaviour  of  saturated  soils  are 
reviewed.  The  discussion  spans  the  classical  approaches  to  quasi¬ 
static  and  dynamic  response  of  deformable  porous  media  as  well  as 
theories  of  mixtures  and  their  application  in  geomechanics.  Different 
approaches  to  description  of  motion  of  superposed  continua,  derivation 
of  balance  laws  and  development  of  constitutive  relationships  are 
compared.  The  notion  of  motion  of  the  mixture  as  a  single  body  is 
questioned.  Alternative  concepts  of  behaviour  of  liquid-solid 
mixtures  are  considered.  The  role  of  thermodynamic  principle  in 
development  of  constitutive  relationships  is  examined. 
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NOTATION  AND  LIST  OF  SYMBOLS 


Both  the  standard  indicial  notation  and  matrix  notation  have  been 
used  in  this  report.  Indices  in  Roman  characters  take  on  the  range  of 
values  1,  2,  and  3.  Repeated  indices  imply  summation  over  the  range 
of  index  values.  Pairs  of  indices  in  parentheses  represent  the 
symmetric  part  of  the  tensor,  while  square  brackets  denote  the  anti¬ 
symmetric  part.  A  superposed  dot  is  used  to  indicate  differentiation 
with  respect  to  the  time  variable.  A  subscripted  comma  indicates 

differentiation  with  respect  to  the  coordinates  defined  by  the 

subscripts  following  the  comma.  6-.  is  Kronecker's  delta  and  e^i, 

denotes  the  permutation  symbol.  Indices  used  elsewhere  as 

mathematical  symbols  and  summation  parameters  have  significance 

evident  from  the  context. 

For  matrix  notation,  {•}  denotes  a  column  vector  and  [•]  a 
rectangular  matrix.  Diagonal  matrices  are  indicated  by 
Superscript  T  is  used  to  denote  'transpose'. 

Generally  speaking,  symbols  have  been  defined  as  they  appear  in 
the  text.  The  list  given  below  is  not  exhaustive, 
a  Constant 

a"  Constant 

a^  Constants,  Initial  compressibility 

av  Coefficient  of  compressibility 

A  Surface,  Helmholtz  free  energy 

Aj ^  Components  of  the  matrix  of  phenomenological  coefficients 

A1jkl  Soil -water  interaction  tensor 

b  Constant,  total  body  force 

bj  Constant 

b-  Components  of  the  body  force  vector 

B  Body 

c  Constant,  Mass  production  factor 

Cj  to  Cg  Constants 

cv  Coefficient  of  consolidation 

[C]  Viscous  damping  matrix 

Cj 1  Components  of  the  flow  resistivity  tensor 


Components  of  the  compliance  tensor 


Cijkl 

C-ijkl^*  Components  of  the  intrinsic  compliance  tensor  of  solid 

dW  Reversible  work  of  fluid  phase. 

d^j  Components  of  rate  of  deformation  tensor 

D  Phenomenological  coefficient 

D.j  Components  of  the  diffusive  force  vector 

D;  .•  Components  of  the  spatial  rate  of  deformation  tensor 

^  ic 

Components  of  the  intrinsic  rate  of  deformation  tensor 
e^j  Infinitesimal  strain  tensor 

ekk  Volumetric  strain 

E  Young's  modulus 

Ep  One-dimensional  rebound  modulus 

Ekl^j  Components  of  the  elasticity  tensor  for  dry  solid 

f.  Components  of  the  acceleration  vector 

F  Deformation  gradient 

Fj  Prigogine  forces 

F^j  Components  of  the  material  deformation  tensor 

F^n  Total  force  function 

h.j  Components  of  the  heat  flux  vector. 

H  Constant 

Constant 

J  Jacobian,  Bulk  volume  measure 

K  Equilibrated  inertia.  Bulk  compressibility,  Bulk  modulus, 

Consolidation  factor.  Constant 
Kj  Intrinsic  bulk  modulus  for  solid 

K2  Intrinsic  bulk  modulus  for  fluid 

K$  Intrinsic  compressibility  of  solid 

f  1 1* 

Kp  '  Bulk  modulus  of  rock 

1 ^  Compressibility 

1  External  equilibrated  body  force 

Components  of  the  spatial  velocity  gradient 
L-jj^  Scalar  coefficients 

m1  Momentum  supply  density 

M  Material  parameter,  Constant 

-vi i i - 


(k) 


{P(t)l 


(k) 


(k) 
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ij 
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at 
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lk 
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Components  of  the  material  parameter  tensor,  Components  of 

the  skew-symmetric  matrix 

Diagonal  mass  matrix 

Porosity,  Volume  fraction 

Volume  fraction  of  constituent  s^k^ 

Initial  porosity 
Constant 
Fluid  pressure 
Excitation  vector 
Axial  compressive  load 

Components  of  the  interaction  force  or  diffusive  resistance 

Surface  flux  of  energy 

heat  flux  vector  of  constituent  s^k^ 

Constant,  Unitary  operator 

kth  constituent  of  mixture 

Entropy  per  unit  mass  of  mixture 

Components  of  the  equilibrated  body  force  vector 

Components  of  the  surface  force  tensor 

Total  traction 

Total  stress  tensor 

Effective  stress  tensor 

Partial  stress  tensor  for  constituent  s^^ 

Intrinsic  stress  tensor  for  s^Hh  constituent 

Total  kinetic  energy,  Temperature  of  mixture 

Displacement  vector 

Velocity  vector 

Acceleration  vector 

Diffusive  velocity 

Molecular  energy.  Specific  total  internal  energy  of  mixture 

Weighted  average  of  energies  of  constituent 

Relative  fluid  velocity 

Components  of  the  velocity  vector 

Components  of  barycentric  velocity 

Partial  volume  of  constituent  s^) 

Components  of  the  nominal  relative  velocity 
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Surface  interaction  force  of  constituent 
Constant 

Components  of  internal  force 
Constant 

Components  of  the  vorticity  tensor 
Differential  change  in  volume 
Incremental  pore  pressure 
Relative  strain 
Change  of  water  content 
Interaction  force  of  constituent  s(k) 

Constant,  Lame's  constant 
Constants 

Components  of  entropy  flux 

Coupling  tensor  due  to  interaction 

Difference  in  velocity  tensor  for  a  binary  mixture 

Bulk  shear  modulus.  Chemical  potential.  Constant,  Lame's 

constant 

Constants 

Poisson's  Ratio 
Bulk  fluid  pressure 
Intrinsic  fluid  pressure 

Mass  density  of  the  saturated  soil.  Density  of  the  mixture 
Intrinsic  fluid  mass  density 
Initial  mass  density 

Mass  coupling  between  solid  &  fluid 
Mass  density  of  constituent  s>*' 

Intrinsic  fluid  density  of  constituent  s^k^ 

Potential 

Volume  interaction  term  for  constituent  s^ 
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SECTION  I 
INTRODUCTION 


Study  of  response  of  fluid-saturated  porous  solids  to  applied 
loading/unloading  or  fluid  injection/withdrawal  is  of  considerable  and 
growing  interest  not  only  in  geotechnical  engineering  but  also  in 
underground  disposal  of  liquid  wastes,  geothermal  energy  extraction, 
petroleum  reservoir  engineering,  and  to  evaluate  the  hazard  associated 
with  a  seismic  or  other  significant  dynamic  event  (e.g.,  a  blast  or 
explosion).  It  is  necessary  to  investigate  the  motion  of  the  soil 
mass  in  the  vicinity  of  the  source  of  the  disturbance  and  near  any 
significant  structures  whether  above  or  under-ground.  This 
investigation  is  made  difficult  by  the  presence  of  fluids  in  the  pores 
of  the  soil  mass.  In  the  case  of  water-bearing  deposits,  loading 
results  in  a  build-up  of  pore  fluid  pressures.  Subsequent  dissipation 
associated  with  gradual  fluid  flow  leads  to  time-dependent 

deformation.  Incidence  of  earthquake  motion  or  other  dynamic  loading 
on  fluid-saturated  geological  deposits  results  in  changes  in  the  fluid 
pressure  as  well  as  the  stress  field  in  the  solid  matrix.  This 
phenomenon,  in  some  cases,  can  lead  to  instability  of  the  soil  matrix 
resulting  in  soil  liquefaction  and  disastrous  consequences  to  the 

structures  that  may  depend  upon  that  soil  for  support.  Also,  the 
energy  dissipation  associated  with  the  relative  oscillatory  motion 

between  the  fluid  and  the  solid  may  introduce  significant  attenuation 
of  the  propagating  wave.  This  effect  would  be  in  addition  to  the 
energy  dissipation  in  any  inelastic  deformation  of  the  soil  matrix. 

Inasmuch  as  the  transmissibil ity  of  motion  through  soil  layers  depends 
upon  the  soil  characteristics ,  the  soil  deposit  may  act  as  a  selective 
fi lter/ampl i fer.  The  importance  of  transmission  characteristics  and 
stability  of  the  soil  matrix  in  the  study  of  earthquake  hazard  and 
risk  of  structural  damage  has  long  been  recognized.  For  seismograph 
installations  located  on  soil  deposits,  the  records  include  local  site 
effects  which  have  to  be  allowed  for.  For  reliable  prediction  of 
stability  and  time-dependent  deformation  of  the  system  and  for 
evaluation  of  various  remedial  measures  to  control  or  alleviate  any 
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adverse  conditions  that  might  arise  due  to  natural  causes  or  certain 
engineering  or  other  human  activity,  it  is  necessary  that  adequate  yet 
simple  mathematical  simulation  of  events  be  available.  Various 

theories  and  methods  of  analysis,  including  numerical  procedures,  have 
been  developed  to  calculate  the  surface  motion  response  and  histories 
of  stress  in  fluid-saturated  geologic  deposits. 

Terzaghi  [187]  was  the  first  to  study  the  problem  of  development 
and  dissipation  of  pore-water  pressures  in  embankments  and  in 
foundation  soils  in  the  vicinity  of  structural  foundations.  His 
theory  of  one-dimensional  consolidation  was  extended  by  Gibson  [75]; 
Schiffman  [162];  Barden  [10,11];  and  Fukuo  [68]  to  viscoelastic 
soils.  Gibson  [77]  developed  a  one  dimensional  consolidation  theory 
admitting  finite  strains.  To  use  Terzaghi 's  one-dimensional  theory 
for  studying  three-dimensional  problems,  semi -empi ri cal  relationships 
were  proposed  by  Skempton  [179,180].  This  approach  was  inapplicable 
to  cases  of  repeated  loading  and  needed  modification  [113].  Moore 
[129]  developed  a  theory  applicable  to  one-dimensional  consolidation 
of  soils  consisting  of  spherical  or  cylindrical  particles  which  are 
themselves  fluid-saturated  and  consolidate  radially.  Treatment  of 
consolidation  on  the  basis  of  rate-process  theory  was  proposed  by  Wu 
[199],  De  Jong  [50]  proposed  a  statistical  theory  of  consolidation. 
Various  extensions  of  these  basic  approaches  to  allow  for  nonlinear 
material  properties  etc.  were  proposed  by  several  investigators  [e.g., 
48,121].  In  all  this  development,  primary  interest  was  in  surface 
settlements . 

Biot  [15]  was  the  first  to  propose  a  general  theory  of  three- 
dimensional  deformation  of  porous  fluid-saturated  solids.  First 
proposed  for  linear  elastic  isotropic  soils,  this  theory  was  extended 
in  a  sequence  of  papers  [18-20,24,25]  to  anisotropic  elastic  and 
viscoelastic  soils  and  to  systems  under  initial  stress  and  finite 
strain.  This  theory  has  been  extensively  used  [e.g., 

16,17,76,119,163]  to  obtain  solutions  to  specific  problems.  A 
simplified  theory,  known  as  the  Terzaghi -Rendul ic  theory,  based  on  a 
direct  generalization  of  Terzaghi 's  theory  to  three  dimensions,  has 
had  some  advocates  because  of  its  simplicity.  These  theories  have 


been  documented  in  geotechnical  engineering  literature.  As  an  example 
one  may  cite  Zaretsky  [201]'s  monograph,  and  Schiffman's  review  paper 
[164].  Biot  [21-23]  extended  his  theory  to  wave  propagation  in 
saturated  geological  materials.  His  work  has  been  the  basis  for 
several  exact  [e.g.,  21,23,38,52,71]  as  well  as  numerical  [e.g.,  72- 
74,148]  solution  procedures. 

Balance  laws  for  mass,  momentum  and  energy  and  constitutive  laws 
for  the  constituents  had  been  traditionally  proposed  on  essentially 
intuitive,  ad  hoc  basis  to  explain  phenomenological  observations. 
Typical  of  this  approach  of  writing  the  governing  equations  'by 
inspection*  are  the  several  theories  which  are  used  widely  in  soil 
mechanics.  Theories  of  mechanical  behavior  of  mixtures  based  on  basic 
principles  of  mechanics  were  developed  by  Truesdell  and  Toupin  [189], 
Truesdell  [190]  and  Green  [81-88]. 

Truesdell  [189],  Adkins  [1],  and  Green  and  Adkins  [81]  proposed 
the  equations  of  linear  momentum  and  continuity  as  postulates.  Green 
and  Naghdi  [83],  following  Green  and  Rivlin's  [82]  work,  derived  the 
equations  of  balance  of  momentum,  mass  continuity  and  heat  balance 
from  the  frame  invariance  of  a  rate  of  energy  equality.  The  theory 
was  restated  [87]  and  compared  with  the  classical  viewpoint.  Crochet 
and  Naghdi  [45]  specialized  the  theory  to  the  case  of  fluid  flow 
through  an  elastic  solid.  Green  and  Naghdi  [83],  postulated  an  energy 
equality  for  each  constituent  admitting  at  first  all  dependent 
variables.  Local  equations  of  mass  and  momentum  balance  were  derived 
by  repeated  use  of  Noll's  [137]  principle  of  frame  indifference. 
Beatty  [12]  showed  that  Green  and  Rivlin's  method  is  essentially  based 
upon  Noll's  [137]  principle  of  frame  indifference  of  mechanical  power 
and  does  not  require  the  energy  equation  at  all.  Essentially,  both 
the  concepts  are  correct  and  closely  related.  Tabaddor  [184,185] 
followed  the  line  of  thought  presented  by  Green  and  Naghdi  [83]  and 
Crochet  and  Naghdi  [45].  Shi,  Rajagopal  and  Wineman  [177]  applied 
Green  and  Naghdi 's  theory  to  flow  of  a  fluid  through  a  nonlinear 
elastic  solid.  Schiffman  [165]  used  Green  and  Steel's  approach  [84] 
to  develop  equations  of  soil  consolidation  and  to  define  the 
components  of  stress  of  a  porous  solid.  In  application  to  soil 
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mechanics,  Sandhu  [156,158]  sought  to  introduce  Terzaghi's  effective 
stress  concept  and  to  regard  the  soil  as  buoyant  in  the  fluid.  It  was 
noticed  by  Green  and  Naghdi  [83]  that  even  though  the  total  stress  had 
to  be  symmetric  in  the  absence  of  body  couples  applied  to  the  mixture, 
the  partial  stress  for  individual  constituents  need  not  by 
symmetric.  Other  theories,  based  essentially  on  the  pioneering  work 
of  Truesdell  and  Green,  have  been  proposed  in  recent  years  [e.g., 
Schiffman  [165],  Morland  [130-132],  Aifantis  [4],  Prevost  [147-148], 
Vardoulakis  [194],  Katsube  [102],  among  others].  An  excellent  review 
of  various  theories  was  recently  presented  by  Bedford  and  Orumheller 
[14]. 

The  existing  theories  may  be  divided  into  two  groups.  Within 
each  group  there  is  considerable  divergence  of  opinion.  One  group 
uses  equations  of  balance  and  constitutive  relationships  based  on 
purely  phenomenological  observations.  The  other  class  of  theories 
would  use  basic  principles  of  mechanics  and  thermodynamics  to  develop 
the  governing  equations.  Biot's  theory  is  the  one  most  widely  used  in 
soil  mechanics.  However,  it  is  based  on  several  ad  hoc  assumptions. 
Theories  of  mixtures,  proposed  by  various  investigators,  are  based  on 
different  assumptions. 

Reliable  analysis  of  the  motion  and  stability  of  fluid-saturated 
solid  deposits  subjected  to  dynamic  loading  involves  three  steps: 

i.  Formulation  of  the  equations  of  dynamic  equlibrium. 

ii.  Correct  representation  of  material  behavior. 

iii.  Exact  or  approximate  solution  of  the  problem. 

In  view  of  the  apparent  diversity  of  opinions,  postulates,  and 

assumptions  made  in  setting  up  various  approaches  to  the  problem,  it 

seems  necessary  to  examine  the  available  theories  to  determine  their 

theoretical  foundation  and  the  assumptions  made  in  their 
development.  All  the  three  components  listed  above  need  to  be 
addressed.  The  present  report  contains  a  review  of  the  theoretical 

formulation  of  the  problem  of  quasi-static  as  well  as  dynamic  behavior 
of  fluid-saturated  soils.  In  Section  II  we  review  earlier  theories  of 


flow  through  porous  deformable  media  including  dynamic  effects. 
Sections  III  and  IV  cover  theories  of  mechanical  behavior  of 
mixtures.  The  basic  assumptions  of  these  theories  are  reviewed  in  the 
context  of  their  application  to  the  problem  of  fluid-saturated 
solids.  Section  III  contains  the  various  approaches  used  to  describe 
the  motion  of  superposed  continua  and  to  set  up  the  equations  of 
balance  of  mass,  linear  momentum,  angular  momentum,  and  energy. 
Section  IV  is  devoted  to  review  of  the  procedures  for  settng  up 
constitutive  relationships  for  fluid-saturated  porous  media  including 
phenomenological  approaches  as  well  as  those  based  on  theories  of 


mixtures. 


SECTION  II 
EARLIER  THEORIES 


S 


ft 

I 

I 


ft 


i 

3 


I 

Si 

* 

3 


“V 


* 


// 


2.1.  SLOW  FLOW  OF  A  FLUID  THROUGH  A  SOLID. 

2.1.1.  d'Arcy's  Law  and  its  Generalizations. 

Linear  dependence  between  nominally  one-dimensional  water  flow 
and  pressure  gradient  was  observed  by  d'Arcy  £47].  Classical  approach 
to  analysis  of  fluid  flow  through  porous  media  uses,  as  the  basic 
equation  of  flow,  a  generalization  of  this  observation  to  two  and 
three  dimensions  as  well  to  coupled  flow  and  diffusion  phenomena. 
Symbolically,  in  one-dimensional  flow,  d'Arcy  equation  is: 


v 

x 


(II-l) 


where  v  denotes  the  macroscopic  or  nominal  fluid  velocity  and  <j>  is 
* 

the  gradient,  along  the  direction  of  the  mean  flow,  of  the  potential 
causing  the  flow.  The  constant  of  proportionality  K  has  been  found  to 
vary  with  the  size,  shape,  distribution,  and  orientation  of  pores  and 
with  temperature,  viscosity  and  velocity  of  the  fluid.  Also, 

nonlinearities  in  the  relationship  have  been  noticed  and  empirical 
attempts  made  to  devise  formulas  fitting  the  observed  data. 
Comprehensive  reviews  are  available  in  the  literature  [e.g.,112]. 


In  extention  of  d'Arcy's  observation  to  three-dimensional  flow, 
the  constant  K  was  enlarged  to  have  the  nature  of  a  second  rank  tensor 
transforming  the  potential  gradient  to  the  flux  vector.  Thus,  using  a 
cartesian  reference  frame. 


( 1 1-2) 


Here,  and  in  the  sequel,  the  indicial  notation  is  used.  Summation  on 
repeated  indices  is  implied  except  where  indicated  otherwise  by  an 
explicit  statement  to  the  effect  or  by  enclosing  each  of  the  indices 
in  parentheses.  An  index  following  a  subscripted  comma  indicates 
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differentiation  with  respect  to  the  spatial  coordinate  represented  by 
the  index.  The  range  of  indices  is  one  to  three  representing  spatial 
coordinates  except  when  indicated  otherwise  by  the  context. 


The  permeability  tensor  satisfies  the  symmetry  requirements 


*u  ■  kji 


(II -3) 


Nonhomogeneity  of  the  solid  and  spatial  variations  in  fluid 
properties  have  been  allowed  for  by  assuming  the  components  to  be 
spatially  varying.  Their  tensorial  character  also  permits 
consideration  of  hydraulic  anisotropy.  In  case  where  the  solid  as 
well  as  the  fluid  are  in  motion,  d'Arcy's  law  is  taken  to  apply  to  the 
relative  velocity  of  the  fluid  with  respect  to  the  solid  matrix  [e.g., 
18,77]. 


In  the  so-called  thermoelastic  theory  of  consolidation,  Schiffman 
[166]  wrote  the  following  equations  for  coupled  mass  and  heat  flow 
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Mokadam  [125-127]  studied  the  thermodynamics  of  d'Arcy’s  law 
under  multicomponent  flow.  Setting  up  an  experession  for  rate  of 
increase  of  entropy,  Mokadam  wrote,  for  the  flow  of  n  fluids  through  a 
rigid  solid: 
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where  ,mi ,v.j ,  Oi  denote,  respectively,  components  of  the  heat  flux, 
the  diffusive  flux,  the  mass  velocity,  and  the  diffusive  force  vector. 
iim  is  the  chemical  potential  of  the  mth  constituent.  For  isothermal 
flow  in  the  absence  of  chemical  reactions,  this  equation  merely 
indicates  the  temperature  dependence  of  the  permeability  tensor. 

The  d'Arcy  fluid  flow  equation  has  been  generalized  still  further 
[79,80],  and  indeed  forms  part  of  the  general  phenomenological 
equations  of  force-flux  given  by  Onsager  [143].  These  express  the 
effect  of  simultaneous  presence  of  fields  of  mechanical  pressure, 
electric  potential,  temperature  and  chemical  concentration.  The 
general  linear  relationship  is  expressed  as 


where  Ji  are  the  fluxes  viz.,  mass  flow,  heat  flow,  electric  current, 
chemical  diffusion,  and  Fj  are  the  Prigogine  forces  [79,143].  are 

components  of  a  positive  definite  symmetric  tensor.  Onsager  [143] 

expected  the  relationship  to  be  symmetric  i.e.,  =  Lj^.  Evidently, 

the  Onsager  equation,  proposed  originally  for  small  perturbations  on 
an  equilibrium  state,  is  a  restricted  type  of  relationship,  assuming  a 
quadratic  form  for  the  entropy  rate  function.  In  general,  Ji  can  be 
treated  as  functionals,  and  may  depend  linearly  or  nonli nearly  on  the 
spatial  gradients,  of  all  orders,  of  the  potential  fields  and  their 
history.  This  has  been  discussed  by  Coleman  and  Mizel  [44]  in  the 

case  of  heat  conduction.  Same  type  of  reasoning  would  apply  to  other 
flux  phenomena,  and  the  only  thermodynamic  restriction  is  that  the 
scalar  product  of  flux  and  force  be  non-negative  [190].  Truesdell 

[191]  discussed  the  limitations  of  the  Onsager  relationships. 

To  allow  for  compressibility  of  the  fluid,  d'Arcy's  rule  was 

further  generalized,  to  the  form 

n  p  v.  =  -K<j>  .  (1 1-7) 

■  f  ' 

♦ 

where  p  is  the  intrinsic  fluid  mass  density  and  n  is  the  porosity 
equal  to  volume  fraction  of  the  fluid. 


2.1.2.  Mass  Continuity 


For  conservation  of  mass  in  arbitrary  fixed  volumes,  in  the 
steady  state  of  saturated  flow  in  rigid  porous  media, 

(p%i)ji  +  n  p*  =  0  ( 1 1-8) 

where  a  superposed  dot  signifies  partial  differentiation  with  respect 
to  the  time  variable.  Substituting  for  velocity  from  the  generalized 
d'Arcy  equation  ( 1 1 -2) , 


For  isotoropic  flow  of  an  incompressible  fluid  through  a  homogeneous 
solid,  the  density  is  invariant  over  time  and  space  and  also 


In  this  case  (II -9)  reduces  to 


( 1 1 - 10 ) 


♦  f1i  =  0  (11-11) 

This  is  the  Laplace  equation  governing  d'Arcy  isotropic  flow  of  an 
incompressible  fluid  through  a  homogeneous  rigid  porous  solid. 

2.1.3.  Analysis  of  Deformation. 

2. 1.3.1.  Scheidegger's  correction. 

To  extend  (II -8 )  to  the  case  of  compressible  media,  Scheidegger 
[160]  introduced  a  correction  treating  porosity  as  a  function  of 
time.  The  corrected  equation  was  referred  to  the  deformed  solid. 
Thus 

(p*  vi},i  =  "  ft  [p*  n(1+k)]  (II_12) 

where  (1+k)  is  the  'consol idation  factor'  for  the  porous  medium.  If 
the  solid  phase  is  relatively  incompressible  and  the  change  in  volume 


is  entirely  due  to  change  in  pore  volume,  in  the  above  equation  we  can 
write 

n( 1+k)  =  nQ  +  k  (11-13) 

where  ng  is  the  porosity  at  time  zero  and  k  is  the  volume  strain  of 
the  solid.  Scheidegger  was  interested  primarily  in  the  flow  of  the 
fluid,  and  made  the  simplifying  assumption  that  solid  deformation 
could  be  defined  by  a  single  variable  k. 

2. 1.3. 2.  Terzaghi's  theory  of  one-dimensional  consolidation. 

Drainage  and  consolidation  of  a  saturated  column  of  soil  under 
constant  load  was  investigated  by  Terzaghi  [187],  Assuming  d'Arcy 
flow  and  linear  elastic  saturated  solid,  he  arrived  at  the  equation 

c  it  ■  n  (no  sum)  (11-14) 
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governing  dissipation  of  fluid  pressure  u.  Here  cy,  the  coefficient 
of  consolidation,  is  defined  as 

cy  =  K/(ay(l-n))  (11-15) 

where  av  is  the  coefficient  of  compressibi  1  ity  and  K  the  permeability 
constant. 

The  only  deformation  of  the  porous  medium  that  interested 
Terzaghi  was  the  displacement  of  the  top  surface  of  the  vertical 
column.  This  was  obtained  by  defining  the  degree  of  consolidation  as 
a  linear  function  of  the  excess  water  pressure  over  the  equilibrium 
(steady  state)  value,  and  integrating  it  over  the  height  of  the 
column,  after  obtaining  the  solution  to  the  diffusion  equation  by 
Fourier  type  analysis.  Terzaghi's  theory  was  extended  to  the  case  of 
varying  permeability  and  compressibi 1 ity  by  several  investigators 
[e.g.,  9,48,77,100,121].  The  governing  equation  was  written  in  terms 
of  uniaxial  strain  ez  as 


( c „  e  )  =  e  (no  sum) 


(11-16) 


In  this  manner,  a  non-constant  cy  could  be  considered.  Actually 
variations  in  permeability  and  compressibility  become  significant  only 
if  the  void  ratio  changes  and  strains  are  appreciable.  However,  all 
these  theories  assumed  small  strains.  Gibson  [77]  extended  the  theory 
to  finite  strain.  Davis  [48]  and  Simons  [178]  extended  the  theory  to 
nonlinear  compressibility. 

A  noteworthy  aspect  of  Terzaghi's  work  was  his  definition  of 
fluid  pressure  and  effective  stress.  For  the  static  case,  the  fluid 
pressure  in  a  saturated  porous  medium  is  measured  by  means  of 
piezometers  and  the  measured  pressure  is  assumed  to  act  over  100%  of 
the  area  of  internal  surfaces.  Water  pressure  measurements 
[5,33,92,103,111,188,200]  in  soils,  rocks,  and  concrete  support  this 
assumption.  Indeed,  stability  calculations  allowing  for  uplift  in 
hydraulic  structures  are  based  on  this  assumption.  It  essentially 
amounts  to  regarding  the  solid  particles  as  having  only  point  contacts 
with  one  another  and  treating  the  solid  as  buoyant  in  the  fully 
connected  fluid.  This  assumption  is  not  valid  for  very  dense  rocks 
where  the  fluid  may  not  completely  enclose  each  particle  of  the  solid. 

2. 1.3. 3.  Biot's  theory  of  three-dimensional  consolidation. 

Assuming  water  to  be  incompressible  and  d'Arcy  flow  through  an 

isotropic  porous  soil,  Biot  [15]  wrote  the  equation  for  mass 

★ 

continuity  of  the  fluid  in  terms  of  the  intrinsic  fluid  pressure  n  as 

K  **u  =  9  (II-17) 

where  9  is  the  change  of  water  content  in  a  fixed  volume  of  the  soil 
expressed  as  the  difference  from  the  initial  state,  and  equals  the 

inflow  rate  given  by  the  left  hand  side  of  the  equation.  The 

superposed  *  indicates  an  'intrinsic'  rather  than  a  'bulk'  quantity. 

Biot  [15]  regarded  the  total  stress  t^  and  the  fluid  pressure  as  the 
mechanical  variables.  Assuming  that  the  hydrostatic  fluid  pressure 
causes  volumetric  strains,  in  the  soil  mass,  which  can  be  superposed 
on  the  elastic  equilibrium  strains  for  no  'excess'  pressure;  and  that 
the  water  content  per  unit  volume  of  the  soil  would  depend  upon  the 
hydrostatic  components  of  stresses  in  the  solid  and  the  fluid,  Biot 


[15]  proposed  the  following  relationships  for  homogeneous,  li,,aar 
elastic  isotropic  soil  undergoing  small  strains: 


„  (  1  +  V  )  *  V*  c  ,  1  _*  c 
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(11-18) 


(H-19) 


Here  e^j  are  components  of  the  strain  tensor  for  the  soil,  H,  ,  and 
R  are  material  constants,  E,  v,  are,  respectively,  Young's  modulus  and 
Poisson's  ratio.  Effectively,  this  amounted  to  introducing  e  as  an 
additional  kinematic  quantity.  Further  assuming  the  existence  of  an 
energy  function  for  the  saturated  soil  in  terms  of  e^j  and  9  ,  Biot 
[15]  set  H^=H.  The  inverse  form  of  (11-19)  was,  for  that  case. 
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(H-21) 


where  \,  p  are  Lame's  constants  and  a  is  related  to  0»  R,  and  H  by  the 
following  equations 


J  "  3(l-2v)H 

and 

1  1  a 

IT  =  ~K  H~ 

Combining  (11-17)  and  (11-21),  Biot  [15]  obtained 


Kn  ..  =  9  =  a  e .  .  +  4- 
,n  n  0 


(11-20),  (11-21)  can  be  rewritten  as 


t..  =  2  ue .  .  +  (\  +  a  Q )  e ,  ,  6-.  -  a  Q66  •  . 


(11-22) 
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Physically,  the  constant  1/R  represents  change  in  fluid  content  per 
unit  'excess'  fluid  pressure  for  no  total  stress  on  the  system.  To 
define  a  Biot  [15]  considered  the  deformation  of  a  uniaxial ly  loaded 
system  under  lateral  restraint  i.e.,  e^  =  e32  =0*  initial 

condition,  e  =  0  gave 


u  =  -  aQe. 


(11-27) 


and  the  total  stress 


-Pq  =  2^e33+(\+a  0)233 


=  (\  +  2^)e 33  +  a  Qe33 


(11-28) 


Here  PQ  is  the  axial  compressive  load  on  the  unit  column.  The  initial 
compressibility  a^  of  the  soil  column  is  given  by 


1  -  0  -  *  .  o  .  2  n 

—  =  — -  =  X  +  2\x  +  a  Q 
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(11-29) 


At  the  equilibrium  condition  when  complete  consolidation  has  occurred, 

★ 

the  excess  pore  water  pressure  vanishes,  i.e.,  %  =0.  For  this 

case,  9  =  ae^  and 


-PQ  =  (\  +  2n)e33 


1 1-30) 


Thus,  the  final  compressibility  a^  of  the  soil  column  is  given  by 


1  -  0  -  ,  +  o 
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(II-31) 


Subtracting  (11-31)  from  (11-29)  to  eliminate  Lame's  constants,  and 
rearrangi ng 


1  +  a  afQ 


(11-32) 


Biot  [15]  combined  the  equations  of  equilibrium,  mass  continuity  and 
the  stress-strain  relations  to  set  up  Navier  type  equations  governing 


consolidation  of  linear  elastic  soils  undergoing  small  deformations: 
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Biot  used  this  analysis  to  obtain  solutions  for  the  standard 
consolidation  test-constant  axial  load  on  a  saturated  laterally 
restrained  column  of  soil -and  specialized  it  for  the  case  of  saturated 
clay  taking  a^  very  small  in  comparison  with  a^,  the  final 
compressibility  of  soil.  This  implied  Q — >  ®,  a—-*  1  and  H  =  R  = 
v  The  analysis  was  later  used  [16,17]  to  find  solutions  for 


2v)  The  ana^sis  was  1ater  used  [16,17]  to  f 
two-dimensional  cases  of  strip  loading  of  half-space. 


While  extending  the  analysis  to  compressible  fluids  and 
anisotropic  elastic  [18]  or  viscoelastic  [19]  solids,  Biot  introduced 
the  volumetric  strain  of  the  fluid  as  the  additional  strain  parameter 
instead  of  the  change  in  water  content  used  in  the  earlier  theory. 
Biot  [24]  later  explained  that  the  two  variables  were  essentially  the 
same.  Also,  the  reference  to  the  stationary  state  used  in  the  earlier 
theory  was  omitted  and  the  force  balance  equation  written,  ignoring 
inertia  effects,  as 


'jt.j'1’  *\l  *  pt>i  ’  0 


(11-34) 


Here  bi  are  components  of  the  body  force  .ector  acting  on  unit  mass  of 
the  saturated  soil.  Superscript  1  refers  to  the  soil.  The  fluid 
pressure  here  is  the  'bulk1  value  equal  to  n^  n  .  Superscript  2 
refers  to  the  fluid.  The  flow  equations  were  based  on  d'Arcy  rule 
applied  to  relative  velocity  of  the  fluid  with  respect  to  the  solid. 
Symmetry  of  the  permeability  tensor  was  based  upon  the  assumed 
existence  of  a  dissipation  function.  The  resulting  equations  were,  in 
the  isotropic  case,  the  same  as  those  obtained  earlier.  The  theory 
was  further  extended  [19]  to  the  case  of  viscoelastic  solids  and  to 
allow  for  initial  stress  [24].  For  linear  elastic  solids  [18]  the 
constitutive  relationship  was  written  as 
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where  ekkv  '  is  the  volumetric  strain  in  the  fluid  and  j ,  C^j,  M 
are  material  parameters.  In  [24]  the  quantity  e^  '  was  again 
replaced  by  a  fluid  content  variable.  The  same  concept  was  extended 
to  the  case  of  finite  elastic  deformations  [24],  The  inflow  into  an 
elementary  volume  was  regarded  as  a  vector  and  it  gradient  regarded  as 
the  quantity  conjugate  to  the  fluid  pressure  in  making  a  contribution 
to  the  energy  of  the  system.  This  would  correspond  to  the  velocity 
gradient-fluid  stress  pair  used  by  Adkins  [2]  and  Green  [81]. 
Westmann  [196]  discussed  this  noticing  that  while  writing  mass 
continuity  relations  care  must  be  exercised  because  the  volume  of  each 
continuum  phase  is  not  the  same  as  the  true  volume  of  each  material. 
For  the  case  of  initial  stresses,  equilibrium  had  to  be  satisfied  in 
the  initial  as  well  as  the  deformed  configuration  realized  after 
incrementation  of  stress.  This  necessitated  introduction  [24]  of 
incremental  stress  and  body  force  intensities  in  local  coordinates 
after  deformation.  Constitutive  laws  for  incremental  stresses  and 
strains  were  expected  to  have  the  same  form  as  (11-35),  (11-36). 


2. 1.3. 4.  Gibson's  Theory. 

Gibson  [77]  removed  the  limitation  of  small  strains.  The 
relative  velocity  was  expected  to  follow  d'Arcy  rule.  The  theory 
allowed  for  nonhomogeneity,  rate  dependence  of  materials,  and 
compressibility  of  the  pore  fluid  and  the  solid.  Consolidation 
associated  with  the  self-weight  of  the  soil  was  allowed  for.  In  the 
one-dimensional  theory  developed,  the  soil  had  a  reference 
configuration  and  motion  of  the  fluid  was  considered  relative  to  the 
current  configuration  of  the  solid.  The  equation  of  equilibrium  of 
vertical  forces  associated  with  the  reference  volume  of  the  solid  in 
the  reference  configuration  is 


(11-37) 

where  X,  x  are,  respectively,  the  reference  and  the  place  coordinates, 
s  is  the  vertical  total  stress  and  p  is  the  density  of  the  saturated 
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soil  in  the  current  configuration.  This  density  is  related  to 
densities  of  the  constituents  by  the  relation 


nU)  0(D*  +n(2)  0<2>*  (1)  +0(2) 


(11-38) 


Here  n^^  is  the  volume  fraction  of  the  constituent  indicated  by  the 
superscript.  As  the  description  of  the  solid  phase  is  unaffected  by 
deformation,  the  equation  of  mass  continuity  for  the  solid  is  simply 


(D*  n  (1)  _  (1  )*n(  1 )  3x 

i  nn  -  p  n  tv 


n  (D  -  Jl>  ax 

p0  "  p  aX 


(11-39) 


where  the  subscript  0  refers  to  the  initial  value,  in  the  reference 
configuration,  of  the  quantity  subscripted.  The  equation  of  mass 
continuity  of  the  fluid,  equating  the  inflow  into  an  elementary  volume 
of  the  solid  with  the  increase  in  fluid  content,  is 


n(2)*vi  +5-  (J2)  J2)*  d*\  -  n 

dl  (n  p  v)  +  at  (n  p  ax}  ~  0 

a  /  (2)  ,  a  ,  (2)  ax.  _  n 

ax  (p  v)  +  at  {p  ax’ *  -  0 


(11-40) 


(11-41) 


Here  v  is  the  velocity  of  the  fluid  relative  to  the  soil.  Allowing 
for  compressibility  of  the  pore  fluid,  Gibson  [77]  used  Scheidegger's 
formulation  of  d'Arcy  law  ( 1 1 - 7 ) .  In  this  equation,  the 
differentiation  is  with  respect  to  the  'place*  or  the  current 
configuration.  Keeping  these  considerations  in  view,  the  flow 
equation  has  the  form 


(2)  (2)*  ax  _  r  an  (2)*  ax  , 

n  p  vax'_KCax  +  p  ax] 


(11-42) 


The  permeability  coefficient  K  may  be  a  function  of  porosity  and 
location  of  the  point  under  consideration. 

Gibson  expected  the  effective  stress  to  depend  upon  the 
deformation  or  the  rate  and  history  of  deformation  of  the  solid 
skeleton.  The  fluid  pressure  was  expected  to  depend  upon  the  fluid 


density  in  an  isothermal  system. 


2.2.  WAVE  PROPAGATION  IN  SATURATED  SOILS. 

2.2.1 .  Prel iminaries. 

Presence  of  water  in  soil  strongly  influences  its  static  and 
dynamic  response  as  well  as  its  strength.  In  addition  to  the  effect 
on  dynamic  amplification  factors,  continued  oscillatory  motion  can 
result  in  undesirable  build-up  of  fluid  pressure  in  the  pores  leading, 
in  some  cases,  to  "liquefaction"  and  "soil  mobility".  This  has  led  to 
interest  in  development  of  methods  to  predict  liquefaction  and  soil 
mobility  during  seismic  and  other  dynamic  events,  and  to  define  the 
liquefaction  potential  of  soils.  Several  different  formulations  of 
the  problem  and  solution  procedures  have  been  used.  Herein  we  review 
some  of  the  approaches. 

2.2.2.  The  Engineering  Approach. 

Seed  and  Idriss  [96-99,170]  used  methods  of  structural  dynamics 
to  solve  the  problem  of  shear  wave  propagation  in  soils.  They  used 
the  finite  element  method,  originally  developed  by  Clough  and  Chopra 
[42]  and  extended  to  oblique  incidence  of  seismic  waves  by  Dibaj  and 
Penzien  [58],  to  examine  the  effect  of  the  nature  of  the  soil  deposits 
upon  surface  motion.  These  methods  had  the  capability  to  account  for 
irregular  surface  geometry  as  well  as  nonlinear  material  behavior. 

To  allow  for  pore-water  pressure  build-up  during  dynamic  loading, 
Seed  and  his  co-workers  developed  a  methodology  and  examined  several 
case  histories  [107,110,153,167-176].  Essentially,  this  approach 
consisted  of  a  finite  element  analysis  of  the  system,  under  an 
acceleration  history  applied  to  a  number  of  points  on  the  system,  to 
evaluate  the  stress  history.  This  was  followed  by  a  laboratory  study 
of  the  material  behavior  under  cyclic  stress  conditions  equivalent  to 
those  determined  from  the  finite  element  analysis.  This  approach  was 
used,  with  apparent  success,  by  Seed  and  his  co-workers  to  analyze  a 
number  of  case  histories.  However,  the  method  was  deficient  in 


several  respects.  It  could  predict  the  initiation  of  liquefaction  at 
different  points  in  the  system  but  could  not  predict  the  distribution 
of  pore  water  pressures  after  the  initiation  of  liquefaction;  nor  was 
it  able  to  estimate  the  extent  of  the  resulting  displacements  in  the 
system.  Furthermore,  in  laboratory  testing  under  cyclic  loads,  pore 
pressures  were  observed  to  accumulate.  However,  in  the  initial 
analysis  used  to  determine  the  stress  history  under  input  excitation, 
variations  in  pore  water  pressure  were  ignored.  Dissipation  effects 
could  only  be  included  through  appropriate  choice  of  damping  ratios. 
The  soil  properties  were  assigned  at  the  beginning  of  the  analysis  and 
assumed  to  represent  soil  behavior  throughout  a  wide  range  of  pore 
water  pressures.  Setting  up  an  "equivalent"  number  of  cycles  to 
reproduce  the  effect  of  a  non-cyclic  stress  history  was  based  on 
"judgment"  and  would  be  suspect.  Subsequent  work  has  been  concerned 
with  improvements  in  the  method  leading  up  to  the  formulation 
described  by  Finn  [66].  This  is  essentially  a  periodic  updating  of 
material  behavior  to  allow  for  the  strain  history  as  well  as  pore- 
water  pressure  build-up  and  dissipation.  In  this  approach,  the 
equations  of  motion,  after  finite  element  discretization  and  lumped 
mass  idealization,  are 

M{u)  +  [C]{  u)  +  [K]{u>  =  {p(t  )>  (11-43) 

where  [M]  is  the  diagonal  mass  matrix,  [C]  the  viscous  damping  matrix, 
[K]  the  nonlinear  stiffness  matrix,  and  p(t)  the 
excitation.  { u > ,  {u},  {u}  are,  respectively,  the  displacements,  the 
velocities  and  the  accelerations  of  the  masses  CMJ.  The  stiffness  was 
assumed  to  be  a  function  of  the  volumetric  strain  and  the  effective 
stress  in  the  soil  [66,117].  Generation  of  pore  water  pressure  was 
assumed  to  be  related  to  volume  changes  in  water  and  soil.  Assuming 
water  to  be  incompressible,  an  incremental  relationship  was  written  as 

An  =  Er  Aekk  (11-44) 

where  An  is  the  incremental  pore-water  pressure  and  Ep  is  the  one¬ 
dimensional  rebound  modulus  of  the  material  at  the  effective  stress 


applicable  to  the  increment.  For  one-dimensional  idealization  of  a 
horizontally  layered  system  subjected  to  shear  wave  at  the  base, 
dissipation  of  the  pore  water  pressure  was  expected  to  be  governed  by 
the  equation 


*  *  Er  Hr\z>, z  +  Er 


(11-45) 


where  w  is  the  unit  weight  of  water  and  K  the  permeabi  1  ity .  Finn 
called  this  the  coupled  theory  of  liquefaction.  The  sequence  of 
occurrences  is  assumed  to  be  as  follows:  shearing  stresses  cause 
volume  changes,  volume  changes  result  in  pore-water  pressure  changes, 
pore-pressure  dissipation  follows  (11-45),  pore-pressures  determine 
effective  stresses  and  effective  stresses  along  with  the  cumulative 
shearing  strain  define  the  effective  shear  modulus  to  be  used  for 
determination  of  displacements  and  stresses  for  the  next  time  step. 


2.2.3.  Biot's  Theory. 

Biot  [18-23]  extended  his  theory  of  consolidation  of  soils  to 
wave  propagation  in  saturated  soils.  For  a  statistically  isotropic 
saturated  material,  the  kinetic  energy  was  expected  to  have  the  form 


T  .i(pUv  (*>,  <»  .Zp12,'1*,.'2)  *  „22v.(2V(2) 


1  1 


(11-46) 


where  p  is  the  mass  coupling  between  the  fluid  and  the 

f  k )  I  k  ) 

solid  v..'  denotes  velocity  of  the  constituent  sv  .  For  no 

relative  velocity, v^  v  1- 2 ;  2T  =  (p11  +  2  p12  +  p22)  v^^v^1). 

Hence 


11  .  ,  12  22 
p  +2p  +p  =p 


(11-47) 


Biot  also  showed  that,  assuming  gradient  of  fluid  pressure  is  the  only 
force  causing  acceleration  of  the  fluid  for  no  relative  velocity, 
Langrange  equation  applied  to  T  gives 


(2)  22x  12 

1=0  +0 


(11-48) 
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12  ,  n 

p  <0 


(11-50) 


and  the  effective  mass  of  soil  in  motion  is 


11  -  J1)  n12 

p  -  p  -  P 


(II-51) 


Similarly, 


22  _  (2)  12 
p  -  P  -  P 


(11-52) 


The  above  formulation  was  assumed  to  be  non-dissipati ve.  For  the 
dissipative  case,  a  dissipation  function  quadratic  in  relative 
velocity,  was  introduced.  This  led  to  the  equations  of  equilibrium  in 
the  following  form 

-  (pV11  W21)"  -  «,W>'  (11-53) 


«,(  ■  (P12ut(1>  ♦p».,(*))~-b(ut»)  -u/2’)-  (..-54) 

Here  u^^,  are  the  components  of  the  displacement  vectors 

associated  with  the  soil  and  the  pore  water  respectively,  and  the 
superposed  dot  denotes  material  time  derivative.  The  coefficient  b 
was  related  to  d'Arcy's  coefficient  of  permeabi 1 ity .  For  low 
frequencies,  it  was  [21,69,71] 


.(2)  ,2 


(11-55) 


where  p  is  the  fluid  viscosity  and  K  the  permeability.  Biot's 
formulation  of  the  coupled  problem  was  the  basis  of  theoretical 
solutions  developed,  among  others,  by  Biot  [21-23],  Oeresiewicz  [52], 


Garg  [71],  and  Chakraborty  and  Dey  [38]. 

While  developing  finite  element  solution  procedures  for  the 
problem,  Ghaboussi  [73],  following  Biot  [19,22,23],  introduced 
relative  volumetric  strain  in  the  formulation.  Defining  relative 
strain  as 


(11-56) 


Biot  [19,20]  wrote  the  constitutive  equations  for  the  flow  of  a 
compressible  fluid  through  a  porous  saturated  linearly  elastic 
anisotropic  medium  as 


'ij  '  Wn(1)  *  M1j  5  (”-57) 

**  *  Mijeij(1)+  M  C  (11-58) 


Introducing  d'Arcy's  law  for  relative  velocity  of  the  fluid  with 
respect  to  the  solid  explicitly  in  the  form 


„<2V2>  K,*, 


(H-59) 


and  assuming  that  the  order  of  temporal  and  spatial  differentiation 
can  be  interchanged,  i.e.. 


C  «  [n<2>(«i#i<2>  -  u.  >i(1))]‘  (11-60) 

=  [n(2)(ui(2)  -  u.(1)  )*],.  (11-61) 

Biot  obtained 

i  =  K*  f1i*  (11-62) 


8iot  [20,22,23]  assumming  =  a  M  6^,  wrote  (11-58)  in  an 
alternative  form  as 


=  a 


M  e 


kk 


(1) 


+  M  C 


(11-63) 
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where  a  is  a  measure  of  compressibility  the  solid  particles.  Further 
•  *  > 

assuming  t  •  •  =  t..  +  ait  6..,  where  t.  .  is  the  effective  stress  and 
•J  TJ*  TJ 

substituting  for  n  ,  he  got 

»1J  ■  EkUj  skk(1)  *“  ekk(1)  ♦<=>  ( 11-64) 

This  form  was  used  by  Ghaboussi  and  Wilson  [73,74]  and  others  for 
development  of  finite  element  procedures.  Garg  [71]  wrote  the 

constitutive  relations  in  the  form 

t.  =  a  e..^  6..  +  c  e.,^  6..  +  2  n(e.  -  i  6  - .  e, 

U  kk  ij  kk  ij  ij  3  ij  kk  1 

(11-65) 


*  =  n(2)  =  C  ekk(1)  +  b  ekk' 


(11-66) 


Here  p  is  the  bulk  shear  modulus  of  the  porous  solid  and  a,b,c  are 
functions  of  the  volume  fractions  of  the  constituents,  the  bulk 
modulus  of  the  porous  solid,  and  the  intrinsic  bulk  moduli  of  the 
fluid  and  the  solid.  Specifically  comparing  Garg's  [71]  equation  of 
equilibrium  along  with  the  constitutive  relations  (11-65),  (11-66), 
with  Biot's  equations  (11-53),  (11-54), 


b  -  M(n(2))2 


c  =  M(a  -  n^ )  n^2' 


a  =  \  +  2/ 3m-  +  M(a  -  n(2)  )2 


(11-67) 


(11-68) 


(11-69) 


where  \,  u  are  Lame's  constants  for  the  porous  elastic  solid.  Of 
course,  the  interpretation  of  the  constants  is  different.  In  Garg's 
theory,  the  constants  are  related  to  the  properties  of  the 
constituents  and  the  volume  fractions  as  follows: 


a  =  n(1)  Kj(l  +  bx) 

b,=  a.  n(2)  (n(2)-  (n(1)K„/K, )  a, f1 


(11-70) 


( ii-7i : 
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Here,  K^,  <2  are  the  intrinsic  bulk  moduli  of  the  solid  and  the  fluid; 
and  K,  |a  are  the  bulk  and  the  shear  modulus  of  the  porous  solid.  In 
Biot's  theory,  the  coefficients  are  phenomenological  coefficients 
directly  measured  in  the  laboratory. 

Ghaboussi  and  Wilson  [74]  stated  Biot's  balance  equations  as 


‘ju  *p>i  =p'V1|-p(2)r<V1,-V2,>i-’<1V1)  *’(2V2) 

( 1 1 1-75) 


-  P(2>  .*  (V1’  -u!2>)  (I 


(11-76) 


where  K  is  a  measure  of  permeability  of  the  soil.  Combining  (11-60), 
(11-62)  and  (11-63),  Ghaboussi  wrote 


/  K<"*  ,11  -  *C0+aUj,j(1)  0  (I.-77) 


The  momentum  balance  equations  (11-75)  and  (11-76)  were  written  in 
terms  of  the  six  displacement  components  and  (u.^-u.^)  . 

There  is  no  dissipative  term  in  (11-75).  For  obtaining  numerical 


solutions  to  the  problem,  Ghaboussi  [73]  introduced  an  arbitrary 
viscous  damping  term  proportional  to  u^^.  A  Rayleigh  type  damping 
was  assumed.  Regarding  (11-75)  as  the  equation  of  motion  for  the  soil 
and  (11-76)  as  the  one  for  relative  motion  of  the  fluid,  we  may 
consider  the  damping  in  (11-75)  to  represent  intrinsic  damping  of  the 
soil  and  the  relative  velocity  term  in  (11-76)  to  represent  the 
damping  associated  with  relative  motion. 


'vV'. 

Jl ’J1  . 
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Garg  [71]  used  Biot's  equations  of  motion,  (11-53],  [11-54]  in 
the  following  form: 


'ij.i'1’ =  p<1)  V1’  +b(“j(1>  ■  “j(2,) 

,  .  =  p<2>  u.<2>  -u.(2>) 

*  J  J  J  J 


(11-78) 

(11-79) 


The  second  term  on  the  right  hand  side  of  these  equations  represents 
the  viscous  coupling  between  the  solid  and  the  fluid.  Ghaboussi  and 
Wilson  presented  variational  formulations  for  the  problem  along  with  a 
finite  element  solution  procedure. 

2.2.4.  Fukuo's  Theory. 


Fukuo  [68]  proposed  a  dynamic  theory  of  deformation  of  a  fluid- 
saturated  solid.  He  assumed  the  fluid  to  be  Newtonian,  the  fluid 
pressure  to  be  the  thermodynamic  pressure  determined  by  its  density 
and  temperature,  d'Arcy  law  to  hold  for  relative  motion  of  the  fluid 
with  respect  to  the  solid,  and  the  fluid  to  be  incompressible.  The 
solid  was  assumed  to  be  linear  viscoelastic.  Values  of  displacement, 
stress,  velocity,  were  averaged  over  elementary  volumes.  Strains  in 
the  solid  were  assumed  to  be  infinitesimal.  In  writing  the  equation 
of  motion  for  the  solid,  the  effect  of  pressure  of  the  fluid  as  well 
as  the  'drag'  force  of  flowing  viscous  fluid  upon  the  solid  was  taken 
into  account.  Equations  of  mass  continuity  were  written  considering 
reference  volumes  of  the  solid  in  the  reference  or  original 
configuration.  This  is  parallel  to  the  development,  proposed  by 
Gibson  [77],  for  the  one  dimensional  quasi-static  problem.  The  theory 
is  essentially  similar  to  Biot's. 


2.2.5.  Comments 

We  note  that  the  'Engineering'  approach  is  cumbersome.  The 
physical  properties  can  only  be  determined  as  a  function  of  the 
complete  number  of  cycles  of  stress  at  a  certain  amplitude.  Thus,  in 
an  explicit  type  solution  scheme,  it  is  possible  to  reduce  the  time 


interval  for  updates  to  only  as  low  as  the  time  period  of  vibration. 
Implicit  schemes  can  use  larger  intervals  and  are  expected  to  be  more 
reliable.  However,  it  has  been  difficult  to  show  convergence  of 
results.  Under  harmonic  loading,  the  time  to  liquefaction  and  the 
location  of  the  layer  that  liquefies  have  been  correctly  determined  in 
some  cases.  However,  the  acceleration  and  displacement  response 
calculation  is  not  reliable.  Post-liquefaction  distribution  of  pore 
pressures  cannot  be  correctly  determined  in  this  theory.  The  approach 
has  only  been  used  for  one-dimensional  wave  propagation.  It  cannot  be 
extended  to  two  and  three-dimensional  cases.  Use  of  this  theory 
requires  considerable  experience  and  'judgment',  in  addition  to 
extensive  laboratory  testing  program,  to  get  useful  results. 

Considering  Biot's  theories,  we  note  several  different  forms. 

Adding  (11-53)  and  (11-54),  noting  (11-48),  (11-49),  and  adding  the 

body  force  on  the  mixture,  we  get  (11-75).  However,  (11-76)  is 

different  from  (11-54)  not  only  in  that  (11-75)  includes  body  force 

effects,  but  also  in  that  the  inertia  terms  are  different.  Of  course, 

the  coefficient  b  in  (11-54)  may  be  identified  with  1/K  in  (11-76). 

Again,  (11-78)  and  (11-79)  can  be  identified  with  (11-53)  and  (11-54) 
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only  if  the  quantities  p  ,  p  vanish;  i.e.,  there  is  no  inertial 
coupling  between  the  fluid  and  the  porous  solid. 

In  writing  the  kinetic  energy  of  the  mixture  (11-46)  Biot 
included  a  term  involving  product  of  the  velocities  of  the 
constituents.  This  implies  that  the  kinetic  energies  of  the 

consti tutents  are  not  additive.  This  assumption  will  be  examined 
further  in  Chpater  III  in  the  context  of  Truesdell's  [189-191]  and 
Green's  [83-88]  work. 


SECTION  III 
THEORY  OF  MIXTURES 


3.1.  INTRODUCTION. 

In  this  section  we  outline  various  schemes  used  to  describe  the 
motion  of  a  mul ti -constituent  material  and  to  set  up  equations  of 
balance  of  mass,  linear  momentum,  angular  momentum  and  energy.  The 
relationship  of  quantities  associated  with  the  constituents  and  with 
the  mixture  as  a  whole  is  reviewed.  Truesdell's  [189,191]  and  Green's 
[83-88]  theories  are  compared.  Other  relevant  work  is  mentioned  as 
needed.  Some  of  the  related  work  is  summarized  towards  the  end  of  the 
section. 

Truesdell  and  Toupin  [189]  wrote  the  local  form  of  the  equations 
of  balance  of  mass,  linear  momentum,  angular  momentum,  and  energy  for 
the  constituents  and  for  the  mixture.  They  introduced  the  notion  of  a 
diffusive  force.  They  also  wrote  the  local  form  of  an  entropy  inequa¬ 
lity  for  the  mixture.  Kelly  [104]  presented  integral  forms  of  the 
balance  equations  and  derived  the  local  forms  from  them. 

In  developing  a  rational  theory  of  mixtures,  Truesdell  [191]  laid 
down  the  following  "three  metaphysical  principles". 

1.  All  properties  of  the  mixture  must  be  mathematical  consequen¬ 
ces  of  properties  of  the  constituents. 

2.  So  as  to  describe  the  motion  of  a  constituent,  we  may  in 
imagination  allow  properly  for  the  actions  of  the  other 
constituents  upon  it. 

3.  The  motion  of  the  mixture  is  governed  by  the  same  equations 
as  is  a  single  body. 

To  ensure  satisfaction  of  these  principles,  the  total  stress  had  to  be 
so  defined  that  it  did  not  equal  the  sum  of  partial  stresses  in  the 
constituents.  Green  and  Naghdi  [83,85,87]  proposed  an  alternative 
approach  assuming  that  the  total  stress  must  equal  the  sum  of  partial 


stresses  in  the  constituents.  A  deficiency  in  the  theory,  pointed  out 
by  Muller  [136],  was  corrected  [88]  to  admit  interacting  partial 
stresses  in  the  constituents  which  may  not  contribute  to  the  equations 
of  motion  or  the  total  stress. 

Other  theories  of  mixtures  have  been  proposed.  Our  interest  is 
primarily  in  mixtures  of  immiscible  constituents  so  that  the  results 
might  be  applicable  to  the  problem  of  fluid-saturated  soils.  Of 
immediate  interest  would  be  the  so-called  volume  fraction  theories. 

3.2.  KINEMATICS. 

3.2.1.  Prel imi nari es . 

To  apply  the  principles  of  continuum  mechanics,  it  is  customary 
to  regard  a  fluid-saturated  solid  as  superposed  continua.  We  note 
immediately  that  averaging  of  various  quantities,  kinematic  as  well  as 
mechanical,  associated  with  the  various  constituents,  is  inherent  in 
this  assumption.  The  mixture  is  defined  by  the  current  coincident 
configuration  of  the  constituents.  It  is  assumed  that,  in  the  current 
configuration,  each  point  of  space  is  occupied  by  a  particle  of  each 
of  the  constituents.  This  necessitates  the  introduction  of  'bulk' 
description  of  the  material  instead  of  the  'instrinsic'  description 
which  would  apply  if  the  material  were  the  single  constituent  of  a 
body.  In  any  theory  of  mixtures,  it  would  be  necessary  that  as  the 
volume  fraction  of  one  of  the  constituents  approaches  unity,  and  the 
remaining  constituents  disappear,  the  theory  for  a  single  constituent 
be  realized  as  a  limiting  case. 

3.2.2.  Density  of  Each  Constituent  and  of  the  Mixture . 


If  n 1 


is  the  volume  fraction  of  constituent  s’ 


in  the 


mixture,  its  effective  (we  shall  also  refer  to  it  as  partial  or  bulk) 
density,  when  it  is  regarded  as  one  of  the  continua  occupying  every 
point  in  the  spatial  region  of  interest,  is  related  to  its  intrinsic 

( k  )* 

density  o  as 


treated  as  a  continuum  as  it  does  not  satisfy  the  impenetrability 
postulate  except  in  the  case  of  no  relative  motion.  It  should, 
therefore,  be  kept  in  mind  that  the  total  density  Introduced  above  is 
a  purely  mathematical  entity  having  no  physical  interpretation  except 
in  a  special  case. 

3.2.3.  Motion  of  the  Constituents. 

Several  approaches  have  been  used  to  describe  the  motion  of  the 
constituents  of  a  multicomponent  mixture.  One  is  to  refer  the  motion 
of  each  constituent  to  its  reference  configuration.  Another  is  to 
refer  the  motion  of  all  constituents  to  the  reference  configuration  of 
one  of  them.  Yet  another  is  to  refer  all  motion  to  the  current 
configuration,  which  in  the  same  for  all  constituents.  Superposition 
of  relative  diffusive  motion  of  the  constituents  upon  the  mean  motion 
of  the  mixture  as  a  whole  is  also  used.  It  has  been  pointed  out  in 
the  previous  paragraph  that  the  mixture  consisting  of  coincident 
current  configurations  of  all  constituents  has  only  momentary 
existence  as  a  continuum.  It  does  not  deform  or  move  as  a  body,  in 
general.  For  a  binary  mixture  of  a  solid  and  a  fluid,  some 

investigators  describe  the  motion  of  the  solid  with  respect  to  its 
reference  configuration  but  the  motion  of  the  fluid  is  described  as 
relative  to  the  solid.  Another  procedure  refers  to  a  material  region 
consisting  of  the  same  set  of  particles  of  one  of  the  constituents  so 
that  the  bounding  surface  of  this  constituent  varies  with  time.  If, 
in  this  case,  the  description  of  the  material  particles  is  the  same 
during  motion  as  in  the  initial  configuration,  this  procedure  would  be 
the  convected  coordinate  description  for  the  reference  constituent. 
The  following  quantities  and  notations  are  often  employed.  Using  a 
fixed  rectangular  cartesian  reference  frame,  let  x^k^  (x)  be  the 
place  occupied  at  time  x  by  the  material  particle  described  in 

the  reference  configuration  for  the  constituent  s^.  Then 

xt(k)  (x)  =  x.(k)  (Xi(k),  x)  for  -  .  <  x  <  t  ( 1 1 1-6) 

For  the  position  in  the  current  configuration,  i.e.  for  x  =  t, 


following  Green  [83],  we  use  the  notation 


(k)  (k)  ... 

xi  =  xi'  (t) 


(III-7) 


Then,  the  position  at  any  time  x  can  be  expressed  in  terms  of  the 
current  position  in  the  form 


(x)  =  x.^  (x/k>  ,t,x) 


( 1 1 1-8) 


Material  deformation  F^(k)is  defined  as  the  partial  derivative  of 
the  place  coordinates  x^^'with  respect  to  the  material  coordinates 
XjOO,  i.e., 


F.  . 
1J 


(k)  _  a 


axj 


wr  Aj 


(k)  /y  (k) 


x'*'  (X /*\t) 


( 1 1 1 -9 ) 


For  the  deformation  to  be  possible  in  a  real  material, 


(k) 


J  =  det  [F..v  ']  *  0 

*  J 


(III-10) 


At  a  given  time  t,  consider  the  mixture  consisting  of  particles  of  all 
constituents  having  the  same  configuration,  i.e.. 


(k)  (t)  =  x  ^ 


V  '  (O 


(III-ll) 


(k) 


The  velocity  vectors  at  the  point  x.  ;  are  defined  as  the  material 
rate 


(k)  _  D(k)  x  (k)  _  a  (k),  (k) 
i  ~  Ot  x i  "  aT  xi 


( 1 1 1-12) 


f  k) 

Here  the  differentiation  is  carried  out  keeping  X.v  fixed.  If  the 

«J 


position  x^  is  treated  as  function  of  the  position  in  the  current 
configuration,  the  material  rate  operator  is  written  as 


D(k)  a  .  (k)  a 

"Dt-  -  at  +  vi  ~-“0<T 


( 1 11-13) 


9xi 


The  acceleration  vector  at  time  t  for  constituent  s^k^  is  denoted  by 


V  =-DT-V  =ft  vi  (Xj1  J.t)  (IH-14) 


Components  of  the  rate  of  deformation  tensor,  for  the  material  s(k^ 
at  time  t  are  defined  by,  [45], 


d  (k)  _  D(k)  F  (k)  -  (k) 

1  j  "  ~Bt  hi  j  '  Dt  xi  ,j 


( 1 1 1-15) 


Interchanging  the  order  of  spatial  and  temporal  differentiation. 


(yi(k)) 


=  L. 


im 


(k) 


(F.  (k))_1 
v  jm  ' 


( 1 11-16) 


where  L..^  are  components  of  the  spatial  velocity  gradients.  These 
*  J 

relations  hold  for  each  constituent  but  not  for  the  mixture  as  a 

whole,  because  the  mixture  does  not  have  a  reference  configuration. 

The  symmetric  part  of  L..^  is  the  spatial  rate  of  deformation  and 

I J 

the  antisymmetric  part  is  the  vorticity  tensor  whose  components  are 
given  by 


( 1 1 1-17 ) 


where  the  square  brackets  around  the  pair  of  indices  denote  ‘anti¬ 
symmetric  part'.  The  spatial  rate  of  deformation  is 


( 1 1 1-18) 


Most  investigators,  following  Truesdell  [189]  and  Green  and  Naghdi 

[83]  also  introduce  a  barycentric  velocity  for  the  mixture  and  a  < 

material  rate  for  the  mixture  as  a  whole.  The  velocity  is  defined  as 

< 

Vi  =  l  P(k)Vi(k)/p  ( 1 1 1-19) 

1  k=l  1 


where  p  has  already  been  defined  by  ( 1 1 1 -2 )  as  the  density  of  the 
mixture.  A  material  rate,  executed  on  the  mixture,  is  defined  by  the 
following  argument: 


V  0(M  D(k)  .  y  M  r  a  +  (k)  5 

A  p  T5t - ,A  p  c^t  +  vi 


=  p  at  +  p  vi  ax, 


D 

‘  p  Dt 


( 1 11-20) 


We  note,  however,  that  this  material  rate  has,  in  Atkin  and  Craine's 
[7]  words,  'no  particular  physical  significance'.  It  is  so  because 
the  rate  is  executed  not  on  a  material  particle  but  on  a  center  of 
mass.  Another  quantity,  found  convenient,  is  the  diffusive 

velocity.  This  is  defined  as  the  velocity  of  a  constituent  relative 
to  that  of  the  mixture,  i.e., 


(k)  (k) 

U.  =  V .  -V- 

1  1  1 


( 1 11-21) 


Clearly  the  sum  of  u-j^  over  all  the  constituents  vanishes.  The 
relative  velocity  of  constituent  with  respect  to  constituent  s^k^ 


w  Ik  (1)  (k) 

V.  =  v. v  '  -  vl  ' 


( 1 11-22) 


Biot  [22],  for  the  case  of  binary  mixture  of  a  solid  and  a  fluid, 
introduced  a  nominal  relative  velocity 


.  „C2>  (,.<2>  . 


( II 1-23) 


This  and  its  variants  were  used  by  Ghaboussi  and  Wilson  [73,74], 
Krause  [109]  and  Kenyon  [105,106],  among  others.  It  can  be  easily 
shown  that 


-“'p./*1  -P(I)  «, 


( 1 1 1-24) 


n(2)  p  U.(D  .  .  p(2)  w. 


( 1 11-25) 
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3.2.4.  Some  Measures  of  Deformation  of  Saturated  Soils. 


Noting  that  the  rate  of  deformation  d^j^  is  uniquely  related  to 

the  deformation  gradient  F^1^  and  the  spatial  velocity  gradient 

and  can  in  turn  be  expressed  in  terms  of  and 

*  J  /  k  \  *  J  '  J 

r.  ■  ,  these  last  two  quantities  can  be  used  to  completely  define 

*  \J 

rates  of  deformation.  Superposed  rigid  body  angular  velocity  Q  would 
change  to  +  Q  leaving  the  other  quantities 

unchanged.  Thus,  for  material  objectivity  to  be  satisfied,  the  vorti- 
city  tensors,  in  binary  mixtures,  must  occur  as  the 

( i  \  •  ( 

difference  A-.  =  r. .  -  r..  . 

*  J  *  J  •  J 

The  deformation  gradient  requires  a  reference  configur¬ 

ation.  For  single  constituent  solids,  this  quantity  is  generally  used 
to  introduce  an  objective  measure  of  strain.  For  a  fluid-saturated 
configuration  of  the  solid.  Garg  [69],  Garg  and  Nur  [70]  and  Morland 
[130,131]  among  others,  used  as  measures  of  deformation  of  the  solid, 
the  quantity 


(1) 


( 1 11-26) 


where  u^^=  Thus  e.-.-^  is  the  symmetric  part  of 

( 1 1  1 J 

F(ij)  "  6ij  '  "^is  measure  deformation  characterizes  the  linear 
theory.  In  the  general  case,  other  measures  of  deformation  are  used 
[e.g.  130].  The  deformation  gradient  F^j  satisfies  the  polar 

decomposition  theorem.  Using  direct  notation,  this  means  that  there 
exist  unique  second  rank  tensors  M  and  R  such  that  F  =  R  M  where  M  is 
symmetric  and  R  is  a  unitary  operator  such  and  RTR  equals  identity  and 
the  determinant  of  R  equals  one.  For  an  initially  isotropic  solid, 
the  isotropy  group  would  include  R,  and  it  is  possible  to  define  as  a 
measure  of  deformation  the  quantity 


Fmi  Fmj  6i  j 


( 1 1 1-27) 


Here  Fmi Fmj  are  components  of  the  right  Cauchy-Green  deformation 
tensor  [137].  Often,  the  only  deformation  of  interest  is  the  volume 
change.  In  that  case,  the  change  in  density  is  an  adequate  measure  of 


deformation.  Considering  the  deformation  of  the  solid  constituent 
s(U  in  a  binary  fluid-solid  mixture, 

p(1)  det  [Fij(1)]  =  P0(1)  ( 1 1 1-28) 

where  the  subscript  zero  denotes  the  initial  configuration.  For  small 
deformations 


detCF^'1’]  =  1  H-  ajjt1) 


( II 1-29) 


where  e..^  is  the  volumetric  strain.  The  change  in  density,  ther- 

J  J 

efore,  for  small  strain  is 


( 1 1 1-30 ) 


( o) 

For  the  fluid  constituent  sx  '  ,  the  deformation  is  completely  repres 
ented  by  its  density.  Following  (II 1-30) , 


(2) 


e  (2)  = 
jj  TT2 T 


n  <2>  p  (2) 
1-0  p0  , 
1  '  ,(2)  (2)*  -1 


p-  -  n'  -  p' 

Introducing  an  intrinsic  measure  of  deformation  [130] 


F..<2>*  .  [n<2>/„  (2>]1/3  F..<2> 
ij  0  J  ij 

Morland  [131]  also  proposed,  for  infinitesimal  strain  theory. 


*..<»*  -  eJ1*  -  [n<2>  -  „<*>]/[!  -  n  <2>] 
JJ  JJ  0  0 

p  (1)  -  n(l)  /  n  (1)  , 

-  ejj  -  0  '  po  -  1 

e,.^‘-!/3e  *"*  6, .  -  e,  -*1*  -  1/3  e  ( "  6,  . 
ij  mm  T  j  t  j  mm  tj 


( 1 1 1-31 ) 

( 1 1 1-32) 

( III -33 ) 
( II 1-34) 
(III -35 ) 
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and  [133] 


^(2)*  ^",2)  -  "o(2)^o<2) 


( 1 1 1-36) 


For  a  fluid-saturated  solid,  in  addition  to  e^^)  the 

porosities  n^\  n^)  have  been  used  as  state  variables.  An  intrinsic 


rate  of  strain  of  the  fluid 


n  (2)  n  (2) 

Dij<2)*  •  °u2)  * t  <-W  vr-  <-W  *u  i"1-37' 


was  introduced  by  Morland  [131]  assuming  the  deviation  parts 
(?)*  f2l 

of  D. .  '  ,  and  D. . '  '  to  be  the  same . 

1  J  J 

Another  measure  of  relative  deformation  sometimes  used  is,  e.g., 
Krause  [109]  the  change  in  fluid  content  of  a  fixed  volume  in  space. 
This  amounts  to 


n 


,<2>  •  <n<2> 


(in-38) 


Carroll  [34]  considered  a  fixed  volume  of  the  solid  in  the 
reference  configuration  following  Gibson's  [77]  approach.  The  total 
deformation  was  considered  to  be  the  sum  of  the  strains  of  the  solid 
particles  and  the  change  in  geometry  of  the  pores.  Writing  volume  of 
the  reference  set  of  particles  as 

V  =  yW  +  =  n^  V  +  n^  V  (II 1-39) 


the  deformation  (change  in  volume)  is  given. 


for  constituent  s 


by 


i.e.,  the  bulk  volume  strain  is 


AV  _  An^  t  AV^ 

~'r~7Ty  “7^ 

For  anisotropic  case,  Carroll  and  Katsube  [36]  wrote 


( 1 11-42) 


e,..-,'1),  <U%  n<2>e  <2>*.e  n<2>(e,  <2>*-e,  (111-43) 

as  the  strain  of  the  reference  material  volume  of  the  bulk  solid. 
Carroll,  Schatz  and  Yamada  [35]  assumed  the  solid  intrinsic  deform¬ 
ations  to  be  reversible  under  solid  pressure  and  the  pore  volume 
change  to  be  irreversible  under  effective  pressure.  This  would 
explain  cumulative  volume  change  under  cycling  of  load.  Kenyon  [106] 
used  the  specific  volume  of  the  solid  defined  as  the  reciprocal  of  the 
intrinsic  mass  density  as  a  state  variable  related  to  the  pore  fluid 
pressure  n  and  the  solid  deformation  gradient  F^1).  Assuming  spec¬ 
ific  volumes  to  be  independent  of  the  bulk  strain  deviation,  F—^1) 
would  be  replaced  by  the  volumetric  strain  e^1). 

Mindlin  [123],  considering  the  microdisplacement  of  the  solid 
particles  of  a  porous  solid,  introduced  relative  deformation  r .  and 
microdeformation  gradient  g^-  as  kinematic  variables  in  addition  to 
the  macroscopic  deformation  e^-.  Introducing  an  energy  function, 
quadratic  in  these  three  variables,  corresponding  mechanical  quant¬ 
ities  were  defined. 


3.3.  BALANCE  LAWS. 

3.3.1 .  Introducti on . 

Discussing  the  thermodynamics  of  d’Arcy's  law,  Mokadam  [126] 
wrote  the  equations  of  mass,  momentum  and  energy  balance  for  a  fluid 
flowing  through  a  porous  rigid  solid  medium.  Truesdell  [189]  post¬ 
ulated  equations  of  balance  of  mass,  linear  momentum,  moment  of  momen¬ 
tum  and  energy  such  that  the  form  of  the  equations  was  the  same  for 
each  constituent  and  for  the  mixture.  Apparently,  the  existence  of 
the  mixture  as  a  continuum  in  motion  was  implied  in  this  line  of 


thought.  To  accomplish  this  identity  of  form,  the  total  stress 
tensor,  the  total  heat  flux  vector,  and  the  specific  energy  supply  had 
to  be  specially  defined.  These  quantities  did  not  equal  the  sum  of 
the  correspondi ng  quantities  for  the  constituents.  The  energy  of  the 
mixture  was,  however,  equal  to  the  sum  of  energies  of  the 
constituents.  On  the  other  hand.  Green  and  Naghdi  [83,85]  postulated 
the  additive  property  of  stress,  heat  flux,  and  the  energy  supply  and 
derived  the  balance  laws  from  frame  invariance  of  a  rate  of  energy 
equality.  The  energy  density  of  the  mixture  was  seen  to  be  different 
from  the  sum  of  energy  densities  of  the  constituents.  This  was 
attributed  to  interaction  between  the  constituents.  Green  and  Naghdi 
[85]  compared  the  two  viewpoints.  In  the  following  section,  we 
summarize  both.  Recent  work  by  Gurtin  [90,91],  Oliver  [140,141], 
Sampaio  [154]  and  Williams  [197]  uses  a  different  approach  to  the 
theory  of  mixtures.  The  equations  of  balance  of  momentum  and  energy 
they  obtain  differ  from  those  of  Truesdell  [189]  and  Kelly  [104].  In 
addition  to  the  partial  stress  for  each  constituent,  they  obtained 
embedding  stresses  governed  by  additional  balance  of  force 
equations.  Beatty  [12]  defined  a  dynamical  process  as  one  in  which 
the  power  is  frame  independent  and  derived  the  balance  Taws  from  the 
frame  indifference  of  mechanical  power. 

3.3.2.  Truesdel 1 's  Theory. 

For  each  constituent  s(k)  of  the  mixture,  Truesdell  [189,191] 
postulated  the  balance  laws  in  the  point  form.  This  was,  using  the 
notation  of  [85], 


i .  Continuity  of  Mass. 


(Ill -44) 


f  k )  f  k  1 

where  pv  '  ,  cv  '  are,  respectively,  the  density  and  the  mass 
production  fraction  of  the  constituent  s^.  As 


■Tp(w  -r*  vkvk)  -^v0 


-  5  n ( k )  .  /.(k)  (k). 

-  at  p  +  {p  vi  '  ,i 


(III -44 )  can  also  be  written  in  the  form 


(k)  A  (k)  (k)  ,(k) 

15t—  P  +  P  vi,i  "  P  c 


( 1 1 1-45) 


Alternatively  [29],  (III -45 )  is  given  by 


1*-  (P(k)det  F  (k))  =  p  c(k)det  [  F1j(k)] 


( 1 11-46) 


For  no  chemical  reaction  i.e.,  c^k^  =  0,  this  gives  (II 1-28 ) . 


ii.  Balance  of  Linear  Momentum. 


/  (k)  (k) >  ,  (k)  (k)  (k),  „  (k)  m  (k)  (k).  (k)  .... 

(p  /viv  M+(p  'v.v  y)  j-tji  j  -P  ny  +P  bj  ( 1 1 1-47 ) 


Here  tj^k^,  p^b.  p  m^k^  are,  respectively,  components  of  the 
partial  Cauchy  stress  tensor,  the  body  force  vector  per  unit  mass  of 
the  constituent,  and  the  partial  momentum  supply  density  for  the 
constituent.  Using  (II 1-44) ,  (II 1-47 )  can  be  rewritten  as 


p(k)  f.(k)=  t^  j(k)  +  p  m^k)+  p^k)  b.(k^-  p  c(k)  v.(k)  ( 1 11-48) 


Bowen  [29]  proposed  to  replace  the  momentum  supply  term  m^k^  by  nyk' 
+  c^k)  Vi<k>  in  (III -47 )  to  get 


>(k)  fl(k)-tJi#/k>  .pm^)  +  p(k>  b^k) 


( 1 1 1-49) 


(k) 


Of  course,  p  m.v  '  in  (III -49)  is  different  from  the  quantity 
represented  by  the  same  symbol  in  (II 1-45)  and  (III -47 ) . 


'vyyyyy. 


iii.  Balance  of  Moment  of  Momentum. 


t.(k>  -  tJk>  -  p  M..(k> 
1J  V  TJ 


( 1 11-50) 


Here  p  M^'  '  are  components  of  the  skew-symmetric  tensor  describing 
the  partial  production  density  of  the  moment  of  momentum  of  the 
constituent  s^k^.  Bowen  [29]  derived  (III -50 )  from  the  axiom  of 
balance  of  momentum  applied  to  a  fixed  volume  in  space. 

iv.  Balance  of  Energy  Rates. 

P<k>  ^  p(k)r<k>=tf.(k,,(i>.)<k)  ♦  p4u<k>  (11,-51) 

where 

4u<k».e<k»-(uik)-^v.(k)Vl(k),c(k).v.(k)„i(k' 


kM..(k>v.  <k> 

7  ij  j,i 


( 1 1 1-52) 


Here  U^k^,  q^k^,  e^k)  are  respectively,  the  specific  internal 

energy,  components  of  the  partial  heat  flux  vector,  the  partial  energy 
supply  and  the  partial  energy  production  density  of  the  constituent 
s(k).  These  equations  represent  the  balance  of  elementary  volumes  in 
the  current  configuration. 

Bowen  [29]  wrote  the  energy  rate  balance  equation,  for  a  fixed 
volume  in  space  of  the  constituent  s^k\  in  the  form 

ft  /  p(k)  (U<kU  j,.<k>Vj<k>)  dV  + 

j  P<k>  (U<k>+  jv  <k>v  <k>)v.<k>n.dA 
A  ^  *  1  J  J 


-  /*  ***  . 
"  -  *  ^  O  '  ■ 


j  <t1j"1,.i(k,-qj('‘,>  "J  dA 


COm  (k) 


r  r  (k)  (k).  (k)  (k),  (kl  (k)  M  (k) 

+  J  [pv  'rv  '  +  'v-s  ;b.x  ;+  p  ex  +  p  v.x  ;m. 

\j  1  ' 


+  p  c 


<k>  (u<k>  ♦  >  «  <k>  »  <k>)3 


( 1 1 1-53) 


Application  of  the  divergence  theorem  and  substitution  of  equations  of 
balance  of  mass  and  linear  momentum  give  the  point  form 

p(k)  D^U  +  q  (k)  _  p(k)p(k)  =  t  (k)  v  (k)  +  p  e(k)  (in_54) 

Several  alternative  forms  of  the  above  equation  can  be  stated 
[29,189]. 

Truesdell  postulated  balance  equations  for  the  mixture  as  well. 
These  equations  can  be  derived  by  summing  over  the  equations  for  the 
constituents.  Truesdell  introduced  the  following  quantities  to  ensure 
that  the  resulting  equations  for  the  motion  of  the  mixture  have  the 
same  form  as  the  equations  of  motion  of  a  single  constituent.  We 
particularly  note  that  the  sum  of  the  specific  internal  energy  and  the 
kinetic  energy  per  unit  volume  of  the  mixture  is  assumed  to  be  equal 
to  the  sum  of  corresponding  energies  of  the  constituents  contained  in 
the  unit  volume.  The  quantity  u.^^  is  defined  by  (III-21). 


i .  Specific  internal  energy  of  the  mixture 


♦4ui<kVk') 


U  =  ~  l  PV  '( 
p  k  =  l 


0  (U  +  I  vi  vi)  =  i  v/V, 

k  _  1 


( 1 11-55] 


( 1 11-56) 


ii.  Total  stress  tensor 


t .  .  = 


(t  Jk>  -p<k>u  <k>  u  <k>) 


(111-57) 


.  V  /*  (k)  (k)  (k)  k). 

lij  ‘p  Vj  ‘  “ij  _p  V  VJ  > 


( 1 1 1-58) 


Iii.  Total  heat  flux  vector 


q,  ■  Ji(q,<k>+  p(k)  (U(k)*  \  uj(k)uj(k))  u!k)-  tj)(k)uj(k))  ( II 1-59) 

Or,  q,  +p  V,[U  -  t^Vj 

V  r  (kL  OO  ( k ) / ,,( k )  A  1  (k)  (kh  +  (k)  (k)-.  /TTT  cnN 

=  l  [q^  '+  pv  'vV  '(IT  '  +  j  vj  vj  )'  tji  vj  J  ( III -60 ) 

k  —  1 


and  iv.  Specific  energy  supply 


r=F  j1p(k>(p<k)tbi(k>Ul(k)) 

P  (r  +  b1v1)  -  J  p(k)(r(k)t  b,(k)  v(<k>) 
1  1  k=l  • 


( 1 11-61) 


( 1 1 1-62) 


Further,  it  was  assumed  that  the  mixture  is  thermodynamically 
isolated,  i.e.. 


I  c' 

k=l 


( 1 1 1-63) 


l 

k  =  l  ‘ 


( 1 11-64) 


1  M  -  0 
k  =  l  ,J 


( 1 1 1-65) 


i  e(k'-0 


( 1 1 1-66) 


Summing  the  equations  of  balance  of  the  constituents  over  all  of  them 
gives,  in  view  of  the  definitions  introduced  above  for  the  total  quan¬ 
tities: 


i.  Balance  of  Mass 


if  * 


■  "  pc<k>  -  0 
k=l 


ii.  Balance  of  Linear  Momentum 


at  (p  V1>  *  (p  Vj’.j 


■*Jid  *  jx  p(k>  bi'k) 


Bowen  [29]  introduced  a  body  force  for  the  mixture  through 
tion 


h  -  V  0<)  h  (k) 

P  bi  =  l  P  V  ' 

1  k=l  1 

Then  the  equation  of  balance  of  linear  momentum  is 


dt 


(p  V,)  ♦  (p  Vj)tj  =  tjfJ  *  p  b, 


An  equivalent  form,  obtained  using  (III -67 ) ,  is 

iii.  Balance  of  Angular  Momentum 


t.  .  =  t 
U  JT 


i.e.,  the  total  stress  tensor  is  symmetric.  As  v.-^)  v.^ 

J 

ric  and  so  is  v-jVj,  it  follows  from  (III -58)  that 


"it  <k>  =  y  t  <k> 

kil  ’J  kil  J' 


iv.  Energy  Balance 


DU  A  5  (k),  (k) ,  .  (k)  (k). 

0  i5t  +  qi,r  lij  vi,r  0  r  =  ki  p  (r  *  b<  V 


i  i 


( 1 1 1-67) 

( 1 11-68) 
the  defini  - 

( 1 11-69) 

( 1 11-70) 

( 1 1 1-71) 

( 1 1 1-72) 
is  symmet- 

( 1 1 1-73) 

( 1 1 1-74) 


Chao,  Sha  and  Soo  [39]  noted  that  this  approach  of  combining  the  bala¬ 
nce  equations  of  constituents  to  obtain  the  balance  equations  of  the 


mixture  can  lead  to  errors.  As  an  example,  they  cited  the  absence  of 
inertial  coupling  forces  in  the  momentum  equations  of  the 
constituents.  However,  their  own  work  has  been  criticized  by  several 
investigators.  The  review  in  reference  [14]  contains  a  discussion  of 
Chao  et  al.'s  theories. 

The  analysis  presented  above  gives  the  equations  of  balance  of 
mass,  linear  momentum,  and  energy  for  the  mixture  as  if  it  were  a 
single  continuum  moving  with  acceleration  derived  from  the  barycentric 
velocity.  We  note  that  the  analysis  above  is  founded  on  the  so-called 
fundamental  identity  [189]  involving  'material  derivatives  of  the  mean 
value'.  Whereas  these  can  be  accepted  as  hypothetical  entities  for 
simplification  of  analysis,  it  is  difficult  to  assign  a  physical  mean¬ 
ing  to  them.  This  is  because  the  mixture  consists  of  centers  of  mass 
and  is  defined  only  for  the  particular  instant  of  time.  It  cannot  be 
regarded  as  a  continuum,  consisting  of  a  set  of  non-penetrating  parti¬ 
cles,  in  motion. 

An  alternative  derivation  of  the  equations  of  mass,  momentum  and 
energy  balance  for  the  mixture  was  presented  by  Kelly  [104],  Truesdell 
[190],  Green  and  Naghdi  [87]  and  Bowen  [29].  They  wrote  the  equations 
of  mass,  momentum  and  energy  conservation  for  the  mixture  contained  in 
an  volume  V  bounded  by  an  arbitrary  fixed  surface  A.  The  conservation 
of  mass  was  expressed  by 

ft  /  l  P(k)dV  +  /  "  n  .  p(k)v  (k)dA  =  0  ( 1 1 1-75) 

oc  V  k-1  A  k=l  J  J 

This  equation  valid  for  any  V,  upon  use  of  the  divergence  theorem 
along  with  ( 1 1 1 -2 )  and  (III-19),  implies  (III -67 ) .  Thus  mass  conserv¬ 
ation  of  a  mixture  is  the  same  as  if  it  was  a  single  continuum  of 

density  p  moving  with  the  barycentric  velocity  v^ .  The  change  of 
linear  momentum  of  the  mixture  bounded  by  a  fixed  surface  A  is 

a  r  °  (k )  (k).„  (  r  (k)  (k)  (k)  .... 

tt  J  A  P  v.v  'dV  +  J  l  n  pv  'v.v  'v.v  'n  dA  ( 1 1 1-76) 

0  V  k=l  A  k=l  J  1  J  J 

Using  (III-2),  (III-21)  along  with  (III -67 ) ,  (III-76)  may  be  rewritten 


/  P]5tvidV+  /  nj  l  P  k  u,  (k)u.(k)dA 
/  vz  1  A  J  k=l  1  J 


( 1 1 1-77) 


and,  using  the  divergence  theorem, 


-vHp{jf.,  ♦  j  P(kVk)uj(k)),JOdv 


( 1 1 1-78) 


Equating  the  linear  change  of  momentum  to  the  total  force  exerted  on 
the  material  inside  the  region  V,  Green  and  Naghdi  [87]  obtained 


/'P5FV  (j  P(l<VkVk)>.j] 


-  /  l  P(kVk>dV  ♦  /  I  <k>dA 

V  k=l  1  A  k=l  J  J 


( 1 11-79) 


Again  using  the  divergence  theorem,  the  change  of  linear  momentum 
equals 


/  l  (P(k)b  (k)  ♦  t„  .(k))dV 
V  k=l  1 


( I 11-80) 


Hence  ( 1 1 1-79 )  can  be  rewritten  as 


I  C 


,  MI  P(k)u^k)u  (k))  ] 

1  k=l  1  J  ’J 


=  J  l  (P(k)bi(k)  +  t..  .(k))dV 
V  k=l  1 


(111-81) 


This  is  the  integral  form  of  the  equation  of  balance  of  momentum.  In 
point  form  the  balance  of  momentum  is,  [87], 


u  x  v  t  OO  (k)  (k) x 

0  ot  vf  ^  <P  V  ’"j  ),j 


l  p(k>b.<k> . "  t.. 

:=1  ’  k>l 


(III -82) 


Green  and  Naghdi  [87]  noted  that  with  the  definition  (II 1-57)  proposed 


by  Truesdell  ,  (III -82 )  is  the  same  as  (III -68 ) -  However,  if  we  assume 
that  the  total  stress  equals  the  sum  of  partial  stresses,  (III -82 ) 
reflects  the  fact  that  the  total  rate  of  increase  of  linear  momentum 
is  not  equal  to  the  barycentric  rate  of  increase  of  momentum  of  a 
continuum  of  density  p  moving  with  velocity  v^ .  Bowen  [29]  wrote  (IH- 
47)  in  the  equivalent  form  (III -49 )  and,  for  the  mixture,  simply 
regarded  (II 1-71 )  as  the  equation  of  linear  momentum  balance.  Similar 
formulation  was  used  for  the  other  balance  equations. 

3.3.3.  Green  and  Naghdi's  Theory 

The  theory  for  a  mixture  of  two  constituents  originally  presented 
by  Green  and  Naghdi  [83,85]  was  generalized  by  Mills  [122]  to  the  case 
of  multicomponent  mixtures.  Green  and  Naghdi  [83,85]  considered  the 
concepts  of  stress,  heat  flux,  and  energy  supply  assumed  to  be 
primitive  to  each  constitutent  to  be  primitive  to  the  mixture  as  a 
whole  as  well.  It  was  proposed  that  the  total  stress  and  the  total 
heat  flux  for  the  mixture  should  equal  the  sum  of  the  corresponding 
quantities  for  the  constituents.  Thus 


«,,  ■  I  t, 

( I I 1-83) 

1J  k=l  J 

«,  •  i  Ok) 

1  k=l  1 

( 1 1 1-84) 

V  00  (k), 
r  =  l  pv  'rv  7p 

( 1 1 1-85) 

k  =  l 

Green  and  Naghdi  [87]  stated  that  Truesdell’s  equations  were  correct 
but  they  had  difficulty  accepting  the  interpretations  associated  with 
some  of  the  quantities  which  occur  in  these  equations.  These 
interpretations  are  of  special  significance  if  boundary  conditions  on 
the  surface  of  the  mixture  involve  total  stress  and  total  heat  flux. 
Morland  [130]  supported  the  additivity  of  partial  stresses  on  the 
ground  that  tractions  are  additive  and  Cauchy’s  stress  principle 
should  hold  for  total  stress  and  total  traction  as  well  as  for  the 
constituents.  In  Green  and  Naghdi's  theory,  the  equations  of  mass  and 
momentum  balance  are  derived  from  the  material  frame  invariance  of  a 
rate  of  energy  equality.  For  a  volume  V  enclosed  by  a  fixed  surface 
A,  following  Truesdell's  contention  that  the  rate  of  energy  equality 


have  the  same  form  as  for  a  single  constituent.  Green  and  Naghdi  [87] 
wrote. 


It  y/  e<u  *  1  ,iv<>  dv  *  J  o"j  vj  <u  +  1  Vl>  dA 

■  J  I  p(k)(r(k)  +  b  <k)v <k>)  dV  +  /  (t^v,  -  q,)n,dA  ( 1 1 1 -86 ) 

Here  the  left  hand  side  represents  the  rate  of  change  in  energy  in 
volume  V  bounded  by  a  fixed  surface  A  plus  the  energy  flux  for  the 

mixture  across  the  boundary.  U  is  the  specific  internal  energy  per 
unit  mass  of  the  mixture  and  is  related  to  internal  energies  of  the 
constituents  by  Truesdell's  equation  (II 1-55 ) .  The  equality  may  be 
regarded  as  written  for  a  surface  moving  with  a  continuum  which  has  a 
velocity  field  equal  to  the  barycentric  velocity  v.j .  Then  the  left 
hand  side  is  equal  to  the  material  rate,  executed  on  the  mixture,  of 
the  sum  of  the  internal  energy  and  the  kinetic  energy  of  the 

mixture.  However,  this  line  of  thought  is  open  to  objection.  It 
would  assume  the  existence  of  mixture  particles  and  the  time  rates 

executed  on  them.  This  is,  in  general,  not  correct  as  the  barycentric 

velocity  is  not  particle  velocity  in  the  ordinary  sense.  Green  and 
Naghdi  [87]  accept  the  form  of  (1 11-86)  and  (II 1-74)  but  not  the 
interpretation  associated  with  some  of  the  quantities  occurring  in 
them. 

Green  amd  Naghdi  [83]  proposed  a  rate  of  energy  equality  in  the 
following  form: 
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(III -87 ) 


Here,  the  heat  fluxes  and  the  energy  supply  are  assumed  to  be 
additive.  U  is  the  internal  energy  per  unit  mass  of  the  mixture 
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.■fjv: 


.'aV-  . 


are 


A 


allowing  for  all  interactions  between  the  constituents.  t^k) 
components  of  tractions  on  the  surface  A.  Essentially,  (III -87 ) 
states 

Rate  of  change  of  energy  in  a  fixed  volume  V  +  energy  loss  due 
to  mass  flow  across  the  boundary  A  of  V  =  Rate  of  energy  supply 
to  the  volume  +  rate  of  work  supply  to  the  volume  +  rate  of  work 
supply  at  the  boundary  A  -  heat  loss  across  the  boundary  A. 

In  this  theory,  multipolar  stresses  and  externally  supplied  multipolar 
body  forces  were  excluded.  Green  and  Naghdi  [83]  also  made  no  attempt 
to  define  the  internal  energy  and  the  energy  supply  for  each 
constituent.  It  was  considered  unnecessary  for  a  complete  general 
theory.  Bowen  [29]  also  pointed  out  that  for  single-temperature 
mixtures,  explicit  use  of  an  energy  equality  for  the  constituents  in 
not  needed. 


In  a  later  paper,  [85],  the  role  of  interactions  between 
constituents  was  made  explicit  by  writing,  for  each  constituent,  the 
rate  of  energy  equality 


V 


If  /  P(kVk>  *  »,  <k>,  <kV  ♦  /  P(k)n,v  <k> 

Ot  y  C  l  l  ^  JJ 

(u^  ^v^y^dA 

=  /p(k>(r(k)  +  b <k>v <k>)dV+/  (t  <k>v  (k>-q(k>)dA 
V  A  1  1 

j  [(*  (k)v  (k)  +  x.  ,(k)  r.  ,(k)  )  +  4,(k)]  dV  -  J  q(k)  dA  (II 1-88) 
/  1  1  A 


where  a.^  ,  •  are  the  internal  force  and  coupling  acting  on  s^k) 

1  1 J  ( k  \  k  1 

due  to  interactions  and  4/  ,  q  represent  contributions  to  the 


balance  of  energy  arising  from  the  interactons  of  the  n  constituents 


( k  1  — ( k  1 

with  <j,v  '  being  a  volume  contribution  and  qv  '  a  surface  flux  of 


energy.  Also 


l  «f(k)-o 

k=l 


(II 1-89) 
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1  Vv  '  =  0 


q<k)=0 

k=l 


( II 1-90) 


( 1 1 1-91 ) 


Without  loss  of  generality,  a...  '  can  be  taken  to  be  [85]  anti- 

( k) 

symmetric  as  .  is  antisymmetric.  Bowen  [29]  and  Bedford  and 
Drumheller  [14]  used  essentially  the  same  form  in  writing  energy 
balance  equality  for  each  constituent.  The  rate  of  energy  equality 
for  the  mixture,  (III -86 ) ,  can  be  rewritten  using  the  following 
relationships : 

D(k)  (k)  (k)  (k)  a  (k)  A  (k)  (k)  A  (k)  (k) 
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=  pc'  'v.  ' v. '  '  +  2  pv  'v.v  'f.v  -(p  v. v  'v.'  'v.  ')  .  (III-93) 

11  ll  11J*J 

Assuming  the  quantities  p,  U,  smooth  in  the  space  and  time 
variables, 

£r  /  PU  dV  +  J  U  (  l  p(k)v.(k))  ndA 
V  A  k=l  J  J 


=  ~  /  pU  dV  +  /  PU  v.n.dA 
ot  V  A  J  J 
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%s  (ow+  u^tPUvj,j>  dv 

■  J  Ce  K  *  U  ♦  U  (p  Vj)tj]<.v 
and  using  (III -67 ) , 

=  /  (P  w  +  I  PUc(k))dV  ( 1 1 1-94) 

V  ut  k=l 

The  rate  of  energy  equality  upon  substitution  of  (III -84 ) ,  (III -85 ) , 
(III -93 ) ,  (II 1-94 ) ,  in  (II 1-87 )  and  noting  the  relationships  (III-2), 
(III-19),  and  ( 1 1 1-44)  is 
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P^-Pr  ♦  I  p^v.^f^-b.^ndV 
^“1 

+  J  l  pc(k)(U  +  |  vi(k,vi(k))  dV 
V  k=l  ^  1  1 

=  /  i  (t^V10 "  ^  dA 

A  k=l  1  1 


( III -95 ) 


Invariance  of  the  energy  equality  under  superposed  uniform 
translational  velocities  yields,  for  arbitrary  V, 

/  l  Pc(k)  dV  =  0  ( 1 1 1-96) 

V  k=l 

i.e.,  the  mass  elements  of  the  mixture  are  conserved.  By  (III -92 )  and 
( III -20 ) , 

w  *  pvi,i  '  0  <"1-97) 
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J  l  [  +  /  l  t/k)dA=0  ( 1 1 1-98) 

V  k=l  11  1  A  k=l  1 

Applying  (III -91 )  to  an  arbitrary  tetrahedron,  defining  partial 
tractions  on  surfaces  perpendicular  to  xk  axis  as  partial 
stresses  t^-^),  a  generalization  of  Cauchy's  stress  principle  is 
realized  i .e.. 


i  Vk)  -  i 

k»l  k=l  J 


( 1 1 1-99) 


Substituting  (II 1-99 )  in  (III -98 )  and  using  the  divergence  theorem,  in 
point  form 

l  [p(k)(b-(k)  -  f.(k))  -  pc<kVk>  +  t..  .(k)]  =  0  (III-100) 

k=l  1  1 

(II I -100 )  is  the  same  as  (II 1-71 ) •  the  equation  of  balance  of  linear 
momentum  for  the  mixture  postulated  by  Truesdell.  Recalling  (1 11-93) 
to  replace  f^*),  (III-100)  may  be  written  in  the  alternative  form 
[83] 

k=1  ’  ’  (III-101) 

This  is  the  equation  of  linear  momentum  balance  for  an  arbitrary  fixed 
surface.  Substituting  (III -97 ) ,  (II 1-99 ) ,  and  (III-101)  in  (III -87 ) , 
Green  and  Naghdi  [83]  obtained,  for  a  binary  mixture, 

/  [Pr  -p™  ♦  |  (P(1>(bi(1>-fi(1))-p^)(bi(2)-fi(2))>  (v/D-v/*)) 


♦  |  tij(vi(1)+  vi(2))>j]dV  +  |  /  (ti(1)-  ti(2))(v.(1)-  Vi(2))dA 


-  J  q  dA  =  0 


( 1 1 1-102) 


Following  Mills  [122],  define 
p.<k>  =  t!k)  -t. ,<k>n. 


Then  (III -99 )  implies 

l  Pt(k)  -0 
k=l 


( 1 1 1-103) 


( 1 1 1-104) 


Let  hj  be  the  heat  flux  across  Xj  planes  at  a  point  x^ .  Applying 
(II 1-95)  to  a  tetrahedron  bounded  by  the  coordinate  planes  at  x^  and 


by  a  plane  with  unit  normal  n  j , 

l  Pi(k)v.(k)  +  (h  n  -  q)  =  0 
k=l  1  J  J 


From  ( 1 1 1 - 104 ) 

Pi(n)=  -TPi(k) 

1  k=l  1 


Hence,  from  (II 1-105 ) 


OO  „  OO 


<  •  Yj  *  J,  7  vi 


.  7  P.<k>,.(k>  *  P. 

k=l  1  1  11 


V  P(<kvk)  -»t,n)  T  P,(k) 

k=l  1  1  1  k=l  1 


V  n  00/  OO  n, 

=  I  Piv  '(vO  '  -  v.v  ') 

k=l  1  1  1 

For  a  binary  mixture,  p.^  =  -  p..^.  Hence 

q  "  hjnj  =  "  vi^2^) 


-I  (P,'11  -P,(2>)  (.,»»  -v,«2>, 


( 1 1 1-105) 


( II 1-106) 


( 1 11-107) 


( II 1-108) 


( 1 1 1-109) 


=  t  .  (k)+  p(k)(b  (k)-f  (k))  -1  c(k)(v.(k)-v,(n))  ( II 1-116) 

J  *  »  1  C.  11 

For  a  binary  mixture,  (II 1-114)  is  identical  to  the  energy  equality 
proposed  by  Green  and  Steel  [84]. 

In  setting  up  the  rate  of  energy  equality,  it  was  stated  that  the 
internal  energy  per  unit  mass  of  the  mixture  may  not  be  equal  to  the 
sum  of  internal  energies  of  the  constituents.  We  note  that  Shi, 
Rajagopal ,  and  Wineman  [177]  assumed  the  internal  energies,  the 
entropies  and  the  Helmholtz  free  energies  of  the  constituents  to  add 
up  to  the  corresponding  quantities  for  the  mixture.  Green  and  Naghdi 
[87]  established  a  relationship  between  these  quantities.  Following 
[85]  we  write  the  rate  of  energy  equality  in  terms  of  the  energies  of 
the  constituents,  ignoring  external  body  couples,  as  follows: 

a  J  l  (P(k)u(k)+  i  P(k)vi(k)v1(l<))dv  ■ 
at  V  k»i  e  11 

♦  /  l  n.(p(k\,  .<k>u<k>*  lp<k>v<k>v  <k>v /k>)dA 

A  k-i  J  1  1  1  J 

=  J  l  p(k)(r(k)+  b.(k)v.(k))dV+  j  l  (t.(k)v.(k)-q(k))dA  ( 1 1 1-117 ) 

\/L_1  '  1  HI.  _  1  *  * 


V  k=l 


A  k=l 


Let  pi)  be  the  sum  of  the  energies  of  the  constituents,  i.e., 

_„*  .  ?  Jk)„(k) 


oU  •  i 
k»l 

Also  define,  [87] 

K  ■  1  <P<  , 

k=l  1 


( 1 1 1-118) 


( 1 1 1 -119) 


J  l  p(k>l/k>  d V  +  j  l  p(k)v.(k)l/k  V  dA 

O'-  w  i.  h  i .  _i  1  1 


V  k  =  l 


A  k  =  l 


V  k=l 


"  (PU*)  dV  +  /  l  p(k)(Vi+  Ui(k))  U(k)ni 


A  k=l 


=  /  [ft  (PU*)  +  I  (p(k)U(k)Vi  +  p(k)ui(k)U(k))>i]  dV 


=  VJ  [ft  (PU  )  +  vi  (PU  ),i  +  PU  vi,i  +  dV 


VJ  C^(pU*)  ♦  PU*  vn  t  +K]dV 


=  [p  W  +  K  +  u  (pvi,i  +  fit  )]  dV 


"  J  Cp?F  +  k+  u*  &  +  (pvi},i]  dV 


and  in  view  of  ( III -67 ) 


«  /  Cp  fif~  +  K]  dV 


( 1 1 1-120) 


If  we  define  the  total  energy,  U,  of  the  mixture  so  that 


DU  DU  .  „ 
p  fit  ‘  p  Dt  K 


^  [(  l  p(k)U(k))/p]  +  l  (p(k)u|k)  U(k)) 
ut  k=l  k=l  1 


( 1 1 1-121 ) 


the  rate  of  energy  equality  (II 1-117) »  using  (III -93 ) ,  becomes 


/  Cp  w  -  1  pv 
v  k*i 


(r<k>.  b.<kV<k>>  *  "  pc<k>(U  av'k>v  <k>) 

"1  I,  _  1  Cl1 


♦  "  p<k>»  <k>f  <k>]  dV  =  /  I  (t  <k>v  'k>-  q<k>)  dA 

i.  ■*  •  1  A  1,  _  1  *  ’ 


A  k=l 


( 1 11-122) 


(Ill -122 )  is  identical  to  (III -95 ) .  Equation  ( II 1-121)  defines  the 
relationship  between  the  specific  energies  of  the  constituents  and 
that  of  the  mixture  for  the  two  rate  of  energy  equalities  (II 1-95)  and 
(III-117)  to  be  equivalent. 

3.4.  OTHER  WORK 

3.4.1 .  Introduction 

Several  investigators  have  used  the  basic  concepts  introduced  by 
Truesdell  and  by  Green  and  Naghdi.  These  efforts  aim  at  simplifying 
the  description  of  motion  for  certain  special  cases.  For  instance,  in 
granular  porous  media,  the  total  deformation  can  be  viewed  as  made  up 
of  two  parts;  one  related  to  deformation  of  the  solid  particles  and 
the  other  to  their  rearrangement  i  .e.  changes  in  pore  geometry.  The 
volume  fraction  theories,  for  compressible  materials,  introduce  the 
volume  fractions  as  additional  variables.  Attempts  have  been  also 
made  to  separate  the  role  of  molecular,  kinetic  and  potential  energies 
in  setting  up  the  balance  equations  and  in  entropy  production 
inequalities.  In  this  section  we  outline  the  essential  features  of 
some  of  these  efforts. 


3.4.2.  Mass  Continuity  Equation  Using  Relative  Velocity 


( 1 1 1-125) 


^k>  ♦  (n<k>v.<k>)  .  -  0 
J  >  J 


For  a  binary  mixture,  (III- 123)  and  (III- 125)  give 


[v.fH  +„<2>(v.<2>  -  ,  =  0 

J  J  J  jJ 


( 1 1 1-126) 


Using  (II 1-23) ,  ( II 1-126)  can  be  rewritten  as  [106,109] 


('j'1’ +  wj>,j  ■ 0 


( 1 1 1-127 ) 


Hsieh  and  Yew  [95]  added  the  equations  (III -67 )  for  the 
constituents  of  a  binary  mixture  and  for  no  chemical  reation  obtained 


P  +  (P  vi(1))>i  +  (p(2)wi)J=  0 


( II 1-128) 


as  the  mass  continuity  equation  in  terms  of  relative  velocity.  This 
is  similar  to  Bowen  [29],  Assuming  small  deformations,  they  also 
wrote  the  continuity  equation  for  the  fluid  volume  contained  in  a 
fixed  volume  in  space  as 


ft  (P(U-P0(1)>  *  Po  ejjUt  l>0<2)  ft  <“i(2)-ui(1)>Ti-  0  C"-29) 


where  u^k^  are  components  of  the  displacement  vector  for  the 

'  /  i  \ 


(k) 

const. tuent  s'  '  .  Integrating  over  time,  for  compressible  fluid  and 


uniform 


(2) 


p  =  pn  (1  -  e. .(1))  +  Pp(2)  I  w,  ,  dt 


Qi  wi,i 


( 1 1 1-130) 


J0  -jj 

3.4.3.  Mass  Balance  in  Terms  of  Porosity  (Volume  Fraction  Theories) 


Fukuo  [68]  used  the  equations  of  mass  balance  to  set  up  equations 
in  terms  of  volume  fractions  of  the  constituents.  The  mass  balance 
equation,  (III -45 ) ,  assuming  no  chemical  interaction,  becomes 


*(k)  A  ,  (k)  (k)»  n 

Pv  1  +  (pv  ;  0^  =  0 


( 1 11-131) 
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it  .*. i  A'.,*  -...V  J 


ft 

m 

$\ 

f 


w 


Using  intrinsic  densities 

•(k)  p(k).  +  „(k)  -(k)*  t  („(k)  p(M*  W)  o  ,m.132) 

•  9  • 

•(!)*  (D* 

If  the  solid  skeleton  in  incompressible,  pv  '  =0,  and  p  .  '  =  0. 

9  « 

Hence,  (III -132)  gives 


;<»  ,  (n^>  ,  <»)  .  =  0 

‘  9  » 

If  the  fluid  is  incompressible, 

*  (n<2>  ,  <2>)  .  -  0 

*  9  » 


( III -133 ) 


( 1 1 1-134) 


Fukuo  [68]  used  Gibson's  [77]  approach  of  referring  to  the  fixed  set 
of  particles  in  the  reference  configuration  to  get 


n0<»  ^  p<‘>*  detCF,,] 


f 1 1 1-135) 


Thus,  for  small  deformation,  using  (III -29 ) , 

n'11  p<»*(l  ♦.„»>)  -n0<D  p0»>* 
Ill* 

If  pv  '  is  constant 

n(1,(l,ejJ(1)).»0(1) 


( 1 11-136) 


( 1 1 1-137 ) 


Then  (II 1-134)  can  be  written  as 

(1.  „(!))•-  n(D-  n0^) 


If  e  .  .(1)  «  1, 

JJ 

n„<»  .  <n<2> 


Srl  1  _,„(2).  (2), 


( 1 1 1-138) 


( 1 1 1-139) 


Li 

■sis. 

■s 

This  equation  combined  with  d'Arcy's  law  may  be  regarded  as  a 
generalization  of  Gibson's  equatio'  (11-42).  Hsieh  and  Yew  [95] 
considered  a  porous  solid  saturated  with  an  incompressible  fluid  and 
undergoing  small  deformations.  For  this  condition,  considering  a  unit 
volume  in  the  undeformed  state,  they  showed  that  the  fluid  content 
change  is 


„  .  .  („<2>.  „  (2),  .  „<2)e  ID 


( 1 1 1-140) 


The  time  derivative  of  (1 11-140),  noting  that  for  continued  saturation 
the  rate  of  change  in  fluid  content  must  equal  the  divergence  of 
relative  velocity,  is 


Hence 


-  If  cn[2)(i  +  ■  («(),, 

ft  <"(2>  •  "0<2)>  *  "o<2)  ft  ejj(11  *  "i,r  0 


( 1 1 1-141 ) 


( 1 11-142) 


This  gives  a  relationship  between  rate  or  porosity  change,  the  rate  of 
volumetric  strain  and  the  relative  velocity  vector.  These  quantities, 
in  this  theory  for  incompressible  fluids  and  no  thermal  or  chemcial 
effects  cannot,  therefore,  be  treated  as  independent  variables.  For 
compressible  fluids,  Morland  [130]  porposed  constitutive  equations  for 
porosity.  Hsieh  and  Yew  [95]  refer  to  (II 1-142)  as  the  compatibility 
condition  of  volume  change  of  a  fluid-saturated  porous  medium. 

3.4.4.  Energy  Balance  in  Terms  of  Porosity  (Volume  Fraction  Theory). 


Goodman  and  Cowin  [78]  postulated  the  equation  of  energy  balance 
for  a  porous  material  with  porosity  n  as 


( 1 1 1-143) 


Here  K,  s j ,  1  are,  respectively,  the  equilibrated  inertia,  components 


“Avs'W 


of  the  equilibrated  stress  vector,  and  the  external  equilibrated  body 
force.  This  equation  admits  an  additional  degree  of  freedom,  viz., 
the  volume  fraction.  A  kinetic  energy  term  is  associated  with  the 
rate  of  change  of  n.  Similarly,  rate  of  work  terms  are  associated 
with  the  rate  of  change  of  n  over  volumes  and  surfaces  using 
generalized  forces  Sj  and  1.  Invariance  of  this  equality,  as  in  Green 
and  Naghdi's  theory,  leads  to  the  equations  of  linear  momentum  and 
angular  momentum  balance.  The  equation  of  linear  momentum  balance,  as 
might  have  been  directly  written,  is 


" »  <V  bi>  *  ‘jfj 


( 1 11-144) 


Goodman  and  Cowin  [78]  also  postulated  an  equation  of  motion,  called 
the  equation  of  balance  of  equilibrated  forces,  for  the  variable  n,  as 


1st  v/  "  p*  u<|r+  1  4  ■  A/  sj"jdA 


( 1 11-145) 


Here  g  is  the  intrinsic  equilibrated  body  force.  The  local  form  of 
( 1 11-145)  is 


n  p*  K  n  =  s.  .  +  n  p*  (1  +  g) 

J  >  J 


( 1 1 1-146) 


Substituting  (III - 144 )  and  (III -146 )  in  (II 1-143 ) ,  the  local  form  of 
energy  equality  is 

"  It  U  '  V(1,J)  +  3  H  -  hj,j  *  "  <In-147> 

3.4.5.  Alternative  Form  of  Linear  Momentum  Balance  Equation. 


Hsieh  and  Yew  [95]  wrote  the  momentum  balance,  for  no  body 
forces,  using  (II 1-76)  and  ( 1 1 1 -79)  as 


VJ  k 


I  p<kVk>dv.j  I  p<kVkVkV -/  i  *ij’  ",dA 

=  1  A  K-l  J  J  A  k*l  J  (,„J-148) 


Using  the  divergence  theorem,  the  point  form  of  the  equation  is 


3_ 

at 


( 1 11-149) 


I 

I 

I 

i 

i 

i 

I 

R 


k=l 


(k) 

p  v. 


(k) 


n 

♦  I 

k=l 


/  (k) 

(p  v. 


(k)„  (X) 


)  ,  =  l  t 


(k) 


k=l 


Assuming  additivity  of  stresses  and  introducing  relative  velocity,  v^ 
v^(2)  _  for  a  binary  mixture,  they  obtained 


p  Dt  vi 


( 1 )  .  ( 2 )  r  a_  +  (2)  +vv  l-t 

p  Lat  i  j  i,j  j  i,jJ  ji,j 


( 1 1 1-150) 


Williams  [198]  summarized  the  work  by  Gurtin  and  De  La  Penha 
[90],  Oliver  [140],  and  Sampaio  [154].  Defining  the  force  exerted  by 
a  volume  C  of  constituent  s^J)  on  volume  D  of  constituent  s^1)  as 
F1-j(D,C)  for  all  D,  C  contained  in  the  body  B  if  i  is  different  from  j 
and  for  all  disjoint  D,C  contained  in  B  if  i=j,  if  prefix  a  denotes 
the  boundary  of  any  volume,  continuity  of  stress  implies 


F..(D,C)  =  /  T..n  dA  for  disjoint  D,C  ( 1 11-151) 

11  aD  n  ac  11 


F.  . (D,C)  +  F. . (D,C)  =  J  T..n  dA  for  disjoint  D,C  ( 1 11-152) 

1J  J1  3D  n  ac 

and 

=  /  -T.  .n  dA  +  /  (p..  +  p..)  dV  for  D  contained  in  C  (II 1-153 ) 

aD  n  ac  1J  d  1J  J1 


and  i  different  from  j. 


r 

E; 

i. 

£ 


Also, 


F  (D)  =  /  t.  dA  +  /  b.  dV 

1  aDraB  D 


( 1 11-154) 


The  equations  of  force  balance  were,  for  the  constituent  s^1')  denoted 
by  subscript  i, 


oiv-  Tti  *  X  pij  *  bi  = 


=  0 


oiv-  Tij  -  <pij  *  pji>  ■  0 


( 1 1 1-155) 
( 1 1 1-156 ) 


1 
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9. 


I 

I 


and 


+  I 
j*l 


( 1 1 1-157) 


=  "  (Sij 


+  Sji> 


( 1 1 1-158) 


where  s--  is  a  surface  force  on  aB.  exerted  by  aB.  e.g.,  a  capillary 

•  J  1  J 

force.  The  quantity  b.  includes  any  inertia  effects.  For  the  global 
body,  define 


T  =  II  T.  ( 1 11-159) 

i  <j  J 

Then 

Div.  T  +  b  =  0  in  B  ( 1 1 1-160) 

and 

T.n  =  t  on  aB  ( 1 1 1-161 ) 

where  b  is  the  sum  of  all  b^  and  t  is,  likewise,  the  sum  of  all  ti 
i.e.  the  body  forces  and  the  tractions  are  additive.  If  partial 
stresses  be  defined  by 


T'  ■  Tn  *  j,  Tn 


( 1 11-162) 


with  \ .  =  1  -  \ ,  then  T  is  the  sum  of  all  T^  and  the  interaction 

J  *0 

force 


pi '  ,L  ?  <pfj '  pji  5 


( 1 11-163) 


The  equations  of  balance  of  the  constituent  s^  are 


Div.  T.j  +  Pi  +  bi  =  0  ( 1 1 1-164) 

Div.  Ti j  -  (Pij  +  Pji)  =0  ( 1 1 1-165) 

and  sum  of  the  forces  p,;  vanishes.  The  definition  of  the  partial 
stress  is  thus  variable.  The  traditional  approach  ignores  ( 1 1 1 - 1 64 ) 
and  constitutive  equations  are  written  for  T.,  and  p^ .  Gurtin  [90]  and 
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Sampaio  [155]  postulated  these  equations  as  constitutive  identities. 
If  FiN(D,C)  is  defined  as  the  total  force  function. 


s 

B 

I 


I 


P 


FiN(D*C)  =  X.  Fij(D*C)  +  Fii(°  °  C’  Cn°) 

J 


( 1 11-166) 


where  D  o  C  is  the  interior  of  the  complement  of  D  in  C.  Additivity 
of  F^n(D,D),  for  a  binary  mixture  leads  to 


F1J(D,C)  =  -  Fij(C,0)  ( 1 11-167) 

which  seems  anti-intuitive  [91,141],  Similarly,  extension  of  the 
traditional  theory  for  single  materials  to  mixtures  by  simply 
replacing  the  forces  by  'total  force'  is  inadequate  to  express  balance 
of  forces  for  other  than  pure  constituents.  In  the  theory  of  balance 
of  forces  proposed  for  mixtures,  continuity  would  apply  to  all  force 
components  F^j(D,C)  and  the  balance  of  forces  would  require 

Fij(D’C)  +  Fji(C’D)  =  0  ( II 1-168) 

to  insure  balance  of  forces  on  any  part  of  B. 


3.4.6.  Other  Forms  of  the  Energy  Balance  Equation. 

Bowen  [29]  postulated  the  point  form  of  the  rate  of  energy 
equality  as 


P 

I 
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I 


ft  cu 
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Vi 


]  = 


{tji 


-  v.i 
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l 

k=l 


(k) 
pv  ' v 


(k)b.(k) 


( 1 1 1-169) 


For  the  case  of  no  relative  motion,  v^V=  for  an  the  constituents 
s^V  Hence 


n 

l  P 
k-1 


(k)b.(k) 


n 

l 

k=l 


(k). 

Pv  bj 


(k)  . 


=  P 


bivi 


( 1 1 1-170) 


In  the  case  when  bV)  is  the  same  for  all  the  constituents  sV,  we 
have  a  similar  results,  i.e., 
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,  / 


V  (k)  (k).  (k)  .  r  (k)  (k)  . 

I  P  v,v  'b.v  =  b-  I  pv  V.v  =  p  b-v. 

k=l  11  1  k=l  1  11 


( 1 1 1-171 ) 


Bowen  pointed  out  the  effect  of  certain  approximations.  If  all 
second  order  terms  in  diffusion  velocities  are  neglected,  (II 1-58 ) , 
(II 1-60 ) ,  and  (II 1-62 )  reduce,  respectively,  to  (II 1-83 ) ,  (III -84 ) , 
and  (II 1-85) .  If  we  also  neglect  products  of  the  diffusion  velocities 
with  their  time  and  space  derivatives,  the  energy  equality  (II 1-169 ) 
becomes 

of'  Vj.i  -"1.1  +pr  ♦  X  P(kVkVk)  U'I-172) 

We  note  that  the  interpretation  of  quantities  appearing  in  (II 1-172 ) 
and  (II 1-169 )  is  quite  different.  Green  and  Naghdi  [87]  wrote  the 
energy  balance  equation  in  the  same  form  as  (II 1-172) .  However,  their 
formulation  was  based  on  different  definitions  for  the  quantities 
associated  with  the  mixture  in  terms  of  those  associated  with  the 
constituents  rather  than  on  any  approximation. 

Using  an  interaction  force  vector.  Green  and  Naghdi  [88]  wrote 
the  energy  balance  equation  as 


DU 

p  Ut  '  p  r 


-"i.i  jt  <"i(kVk)*  tji(k)'i,j(k) 


x  1  00  (k)  (k). 

+  2  pcv  O. v  /uiv  ') 


where 


*  -  l  (p<k»Uj<k'u(k),  . 

k  =  l 


P  u-  l  P(k)u(k) 


t..  .(k)  +  p^)(b.(k)  -  f.f 
JT  .J  M  V  1 


( 1 1 1-173) 


( 1 11-174) 


( 1 1 1-175) 


-  pc'  ' u 


(M„  (k) 


( 1 1 1-176) 


In  an  earlier  paper,  [83],  Green  and  Naghdi  wrote  for  a  bii 


mixture. 


P«if-P  tj,(k\,j(k)^^(1)-i(2))  »,(1)-',<2)) 

1  ( 1 1 1-177 ) 


The  quantity  0..^-  0.^  represents  the  interaction  force  between  the 


two  constituents.  Using  (III-lll),  and  (III-100),  for  a  binary 
mixture 


>  <l>  ♦  e<2>  -  f  ec<k>v<k)  -» 

1  1  k=l  1 


( 1 1 1-178) 


Hence,  the  diffusive  force  is  given  by 


D  =  j  (©^  -  0^2^) 


=  1  (9d)  +9  (1)  -  l  Pc{k)Vi(k)) 
c  1  k=l  1 


=  ei(1)  l 

1  c  k=l  1 


( 1 1 1-179) 


and  again  using  (III-110), 


0,  ->  l  pc(k>vf(k> 

1  1  c  k=l 


'i 


( 1 1 1-180) 


Morland  [131]  stated  the  energy  balance  equation  for  the  mixture  as 


i(k) 


[p(k)  t5t~  u(k)  -  tij(k)vi,j(k)  +  qi,i(k)]  =  0  nn-i8i) 


This  is  the  same  as  ( 1 1 1-177 )  for  no  energy  supply  and  no  interaction 
force. 


Mokadam  [126]  considered  the  total  internal  energy  for  the 
mixture  to  consist  of  three  components, 


E  =  U  +  T  +  L 


( 1 11-182) 


where  U,  T,  L,  are  respectively,  the  molecular,  the  kinetic,  and  the 
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potential  energies  per  unit  mass  of  the  mixture.  Identifying  the 
diffusive  force  as  the  body  force  causing  mass  flow,  Mokadam 
postulated  equations  of  balance  of  energy  in  the  form 


(p  Evj  ♦  qj  -  j  ♦  DjV. 


( 1 11-183) 


where  Dj  are  components  of  the  diffusive  force.  Also 
-  (pLVjlj  -  Pfj*j 


( 1 11-184) 


For  no  chemical  reaction,  (III - 184 )  and  mass  conservation  imply 


OL  , 

Dt  '  fiv7 


( 1 11-185) 


Noting  that,  for  no  chemical  reaction,  mass  continuity  (III -45 ) 
implies 


p  m: =  itr  •  {p  Tvi  >,i 


and  also. 


Dt  =  (pfi  +  Di)vi 


( 1 11-186) 


( II 1-187) 


Mokadam  wrote 


«r=  -  (c  t#i  *  "i>.i 


( 1 1 1-188) 


3.5  COMMENTS 

Various  approaches  to  description  of  motion  of  the  constituents 
and  the  mixture  as  well  as  to  formulation  of  the  equations  of  balance 
of  mass,  linear  momentum,  angular  momentum,  and  rate  of  energy  have 

been  discussed.  We  note  that  in  many  theories  of  mixtures, 

deformation  is  referred  to  an  initial  configuration  for  each 

constituent  and  motion  to  the  place  coordinates.  Also  the  equations 
of  balance  are  written  for  a  fixed  volume  in  space.  This  appears  to 
be  the  correct  approach  for  mixtures  of  fluids  but  may  not  be 


convenient  for  multiphase  mixtures.  The  notion  of  motion  of  a  mixture 
as  a  whole  is  often  introduced.  Indeed,  Truesdell  would  require  the 
form  of  the  equations  of  balance  for  the  mixture  to  be  the  same  as  for 
a  single  material.  We  note  that  for  the  case  of  no  relative  motion 
between  the  constituents,  the  mixture  will  have  motion  and  deformation 
as  a  material  body.  For  that  case  the  development  of  equations  of 
motion  for  the  mixture  is  meaningful.  However,  if  relative  motion  is 
present,  the  mixture  does  not  satisfy  the  axiom  of  continuity  and  its 
corollary,  the  principle  of  impenetrability.  Hence,  the  'third 
metaphysical  principle'  stated  by  Truesdell  [162]  appears  to  be 
irrelevant.  We  note  too  that  the  sum  of  internal  energies  of  the 
constituents  does  not  equal  the  internal  energy  of  the  mixture,  (c.f. 
Biot's  theory  where  the  kinetic  energy  of  the  saturated  soil  is 
expected  to  be  a  quadratic  in  the  velocity  of  the  fluid  and  the  solid 
rather  than  the  sum  of  kinetic  energies  of  the  constituents).  The 
balance  equations  due  to  Truesdell  and  to  Green  and  Naghdi  have 

similar  form  and  are  essentially  equivalent  but  the  quantities 
appearing  in  the  two  sets  have  different  interpretations  based  upon 
the  relationships  postulated  between  the  quantities  associated  with 

the  constituents  and  with  the  mixture.  The  equations  of  energy 
balance  contain  scalar  products  of  corresponding  quantities.  This 
indicates  the  quantities  for  which  constitutive  relationships  would  be 
required.  These  will  be  discussed  in  the  next  section.  Some 

investigators,  considering  the  special  problem  of  fow  through 
deformable  porous  solids,  have  attempted  to  write  the  balance 
equations  in  terms  of  relative  motion  and  porosity,  which  is 
essentially  a  measure  of  relative  deformation.  It  would  appear  that  a 
theory  based  upon  balance  equations  written  for  a  reference  set  of 

particles  of  the  porous  solid  would  be  the  most  appropriate  for  this 
case.  For  the  one-dimensional  case,  Gibson  developed  such  a  theory 
for  the  quasi -static  problem.  This  theory  needs  to  be  developed 
further  to  include  inertia  effects  and  a  generalization  to  three- 
dimensional  problems.  This  would  require  a  convected  coordinate 
description  for  the  stresses  and  deformations  in  the  porous  solid 
along  with  a  referential  description  for  the  motion  of  the  fluid 
relative  to  the  solid. 


SECTION  IV 

CONSTITUTIVE  RELATIONSHIPS 


4.1.  INTRODUCTION. 

In  order  to  set  up  constitutive  equations,  it  is  necessary  first 
to  define  the  mechanical  quantities  for  which  such  relationships  are 
desired  and  to  identify  the  kinematic  or  state  variables  on  which 
these  quantities  might  depend.  For  a  mixture  of  several  constituents, 
the  stresses  in  the  constituents  are  obviously  the  primary  mechanical 
variables.  Truesdell  [189],  following  Maxwell  [118],  recognized 
diffusive  resistance  as  another  mechanical  variable,  reflecting  the 
interaction  between  the  constituents.  Green  and  Naghdi  [83,85] 
pointed  out  the  necessity  of  postulating  constitutive  relationships 
for  this  quantity.  The  rate  of  energy  equality,  e.g.,  (Ill -114) , 
contains  scalar  products  of  'corresponding'  quantities.  This 
indicates  the  need  for  postulating  constitutive  relations  for 
t[ij]^’  and  pi  ^  *  Again,  if  the  equations  of  mass, 
linear  momentum,  angular  momentum,  and  energy  balance  for  consti¬ 
tuents  be  regarded  as  equations  forp^,  v^^,  M-jj^),  and  T^) 
respectively,  constitutive  equations  are  required  for  the  other  quan¬ 
tities  appearing  in  those  equations,  viz.,  c^,  t^j^, 

pJM,  ij(M,  q-j^)  and  also  for  the  partial  entrophy  These 

constitutive  equations  are  subject  to  the  balance  equations  for  the 
mixture  and  to  an  appropriate  entropy  production  inequality.  There 
has  been  some  difficulty  in  defining  components  of  the  partial  stress 
tensors.  Williams  [198]  claimed  that  the  usual  interpretation  of  the 
partial  stress  tensors  is  misleading  and  stated:  "While  it  is 

indisputable  that  the  partial  stress  tensor  represents  a  force  exerted 
per  unit  area,  it  is  not  always  so  clear  exactly  upon  what  the  force 
is  exerted  and  by  what".  Truesdell  [189,190]  required  the  equations 
of  motion  of  the  mixture  to  have  the  same  form  as  those  for  the  motion 
of  a  single  constituent.  In  this  theory,  the  sum  of  partial  stresses 
was  not  equal  to  the  total  stress  in  the  mixture  except  in  the 
approximation  where  second  order  quantities  are  neglected.  Green  and 
Naghdi  [83,85,87],  on  the  other  hand,  started  with  the  assumption  that 
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the  sum  of  partial  stresses  must  equal  the  total  stress.  For  fluid- 
saturated  solids,  the  isotropic  fluid  stress  is  generally  considered 
to  be  the  stress  variable  in  addition  to  the  stresses  in  the  solid. 
Biot  [15]  regarded  the  total  stress  and  the  fluid  stress  as  the 
mechanical  variables.  The  definition  of  pore-fluid  pressure  used  by 
various  investigators  differs  considerably.  Traditionally,  for  a 
water-saturated  soil,  the  pressure  recorded  by  piezometers  inserted 
into  the  water-filled  pores  has  been  assumed  to  be  the  fluid  pressure 
acting  over  100  percent  of  the  area  of  internal  surfaces.  Biot  [21] 
pointed  out  that  the  generalized  forces  defined  by  divergence  of  the 
stresses  are  correctly  defined  by  the  virtual  work  of  microscopic 
stresses  per  unit  value  of  the  displacements  of  the  constituents  and 
not  as  the  average  of  the  microscopic  stresses.  However,  according  to 
Biot,  for  practical  purposes,  either  definition  would  be  acceptable. 
This  is  in  line  with  Mokadam  [125]  who  following  Guggenheim  [89], 
regarded  the  fluid  pressure  to  be  the  thermodynamic  property  such  that 


it  dV  =  dW 


(IV-1) 


where  dV  is  the  differential  change  in  the  intrinsic  fluid  volume  and 
dW  the  reversible  work  of  the  fluid  phase. 

Terzaghi  [187,188],  as  well  as  Green  [82-88]  assumed  the  partial 
stresses  to  act  over  the  entire  area  of  any  surface  element.  Further, 
assuming  partial  stresses  to  be  additive,  the  total  stress 


t..  -tJ1*  ♦t..<2> 

ij  ij  ij 

(21 

For  isotropic  fluid  stress,  t^'  =  %  6^. 


( I V— 2 ) 


(IV-3) 


Biot  regarded  t^^^  as  the  bulk  stress  acting  over  the  entire 
area  of  internal  surfaces.  In  his  earlier  work  [15],  there  was  no 
reference  to  the  area  over  which  the  fluid  pressure  acts.  Later  [16], 
Biot  assumed  the  fluid  pressure  to  act  only  over  the  pore  area.  This 
corresponds  to  the  notion  of  the  fluid  pressure  being  an  intrinsic 


quantity.  Thus 


+  n 


(2) 


**6ij 


( IV— 4 ) 


Garg  [70],  Morland  [130-132],  Pecker  [145]  and  Carroll  [34],  among 

i  k 

others,  introduced  the  notion  of  intrinsic  stresses  t.  .  '  for  each 

^  J 

constituent  leading  to 


t.  .^k)  =  n^  t.  .(k)*  ( IV- 5 ) 

1  j  1  J 

In  reference  [131],  the  coefficient  n^k)  being  area  fraction  could  be 
different  from  the  volume  fraction  and  was  denoted  by  m^k^. 

Tsien  [193]  divided  the  total  stress  into  stress  deviation  and  a 
hydrostatic  component  viz.  t^/3.  The  hydrostatic  stress  was  expected 
to  be  distributed  over  the  solid  and  the  fluid  in  proportion  to  their 
volume  fractions  n^k^,  and  the  solid  was  expected  to  take  the  entire 
stress  deviation.  Carroll  [34]  and  Morland  [130-133],  among  others, 
also  assumed  that  the  shear  traction  is  carried  by  the  bulk  solid  and 
is  related  to  changes  in  pore  geometry.  The  hydrostatic  stress 
components  in  the  solid  and  the  fluid  were  expected  to  cause  intrinsic 
volume  changes  in  the  constituents. 

In  the  theory  of  consolidation  originally  developed  by  Biot  [15] 
assuming  water  to  be  incompressible,  he  regarded  the  change  in  water 
content  of  the  soil  expressed  as  difference  from  the  initial  state  to 
be  the  kinematic  variable  for  setting  up  constitutive  equations. 
Later  [16],  while  extending  the  theory  to  anisotropic  solids  and 
compressible  fluids,  the  volumetric  strain  of  the  fluid  was  introduced 
as  a  kinematic  variable  in  addition  to  the  strains  in  the  solid.  In 
another  work  [24],  the  inflow  into  an  elementary  volume  was  regarded 
as  a  vector  and  its  gradient  was  the  quantity  conjugate  to  the  fluid 
pressure  in  making  a  contribution  to  the  energy  of  the  system.  This 
would  correspond  to  the  velocity  gradient-fluid  stress  pair  used  by 
Adkins  [2]  and  Green  [81].  Biot  [24]  explained  that  the  several 
measures  used  by  him  from  time  to  time  were  essentially  the  same 
quantity.  In  theories  of  mixtures,  densities,  deformat i ons , 


deformation  rates,  velocities,  and  vorticities  of  the  constituents  are 
taken  to  be  the  kinematic  variables.  Morland  [130-132],  Garg  and  Nur 
[70],  and  several  investigators  following  them  introduced  intrinsic 
densities  of  the  materials  as  well  the  porosity  (relative  volume 
fraction  of  the  pores)  as  variables.  Bowen  [32]  regarded  volume 
fractions  as  internal  state  variables  obeying  rate  type  constitutive 
relations. 

It  has  been  difficult  to  define  the  dependence  of  components  of 
the  partial  stress  tensors  upon  the  kinematic  variables  and  densities 
of  the  constituents  of  a  mixture.  In  his  earlier  work,  Adkins  [1,2] 
assumed  that  the  stress  in  each  constituent  depended  upon  the  density 
and  the  kinematic  quantities  associated  with  that  constituent  only. 
Green  and  Adkins  [81]  admitted  interdependence  of  partial  stresses  of 
each  constituent  upon  the  kinematics  of  all.  This  was  in  line  with 
the  principle  of  equipresence,  stated  by  Truesdell  and  Toupin  [189], 
according  to  which  the  mechanical  variables  for  any  constituent  are 
assumed  to  be  influenced  by  the  kinematics  of  all  constituents  unless 
precluded  by  considerations  of  thermodynamics  or  deliberate  assumption 
to  obtain  simplicity.  However,  to  make  the  independent  variables 
distinct  for  various  constituents.  Green  and  Adkins  [81]  stipulated 

that  the  partial  stresses  for  any  constituent  would  depend  upon  the 

densities,  velocities  and  antisymmetric  deformation  of  all 
constituents  but  only  on  the  symmetric  part  of  the  deformation 
gradient  of  the  constituent  itself.  Invariance  of  stress  undor 
superposed  rigid  body  motions  and  under  superposed  uniform  rigid  body 
angular  velocities  of  the  mixture  as  a  whole  showed  that,  for  a  binary 
mixture,  the  velocities  and  the  rotation  tensors  must  occur  as 

difference  terms  of  the  type  v.  =  v.  ^  -  v-^  and 

f  1 1  (2)  11  1 

A..  =  r .  .  '  -  Tj  .  in  the  set  of  independent  variables.  Bowen 

'  J  *  J  J 

[28]  restated  True$dell's  point  of  view  and  dropped  partial  stress 
tensors  from  the  discussion.  In  fact,  following  classical 
thermodynamics,  he  Jid  not  introduce  any  'partial'  quantities.  It  was 
shown  that  in  the  case  of  dilute  mixtures  of  ideal  substances  (dilute 
implying  preponderance  of  a  single  constituent  over  all  the  others), 

the  partial  stresses  can  be  defined  such  that  they  depend  upon  the 
density  and  kinematics  of  individual  constituents  and  their  sum  is  -he 


total  stress  tensor.  This  corresponds  to  Adkins'  earlier  view  [1,2] 
and  as,  in  general,  mixtures  may  not  be  dilute,  the  principle  of 
equipresence  would  appear  to  be  the  most  logical  basis  for  the 

development  of  constituive  theories.  Mori  and  [130]  postulated  that 
the  relation  of  'effective'  stress  to  'effective*  deformation  for  each 
constituent  should  be  its  own  constitutive  law  as  a  single 

continuum.  The  interaction  of  the  partial  constitutive  laws  would  be 
included  in  the  scaling  functions  relating  the  'effective'  and  'bulk' 
quantities.  These  scaling  functions  would  require  a  constitutive 

description  in  terms  of  macroscopic  variables.  For  a  binary  mixture 
of  a  liquid  and  solid,  porosity  was  taken  to  be  the  scaling 
function.  Mokadam  [125]  assumed  the  stress  tensor  to  be  a  linear 

function  of  a  velocity  gradients.  Truesdell  [189]  as  well  as  Biot 
[21]  considered  the  balance  of  momentum  equation  to  be  a  generalized 
form  of  d'Arcy's  rule. 

Thermodynamics  of  d'Arcy's  law,  for  fluid  flow  in  a  rigid  porous 
medium,  was  discussed  by  Mokadam  [125-127],  He  pointed  out  that 
d'Arcy's  law  is  valid  only  for  isothermal  flow  in  which  the  inertial 
and  viscous  effects  are  negligible.  Also  that  the  rate  of  entropy 
production  must  be  nonnegative  separately  for  terms  involving 
quantities  of  different  tensorial  ranks.  Thus,  regarding  stress  as  a 
flux  it  would  need  a  constitutive  equation  in  terms  of  velocity 
gradients  and  the  phenomenological  equations  for  heat,  mass  and 
molecular  flow  will  be  separately  formulated  as  in  ( 1 1 -5 ) .  Truesdell 
and  Noll  [192]  noted  that  a  'proper'  generalization  of  the  Clausius- 
Duhem  inequality  would  be  basic  to  development  of  constituive 
equations  for  diffusion.  Green  and  Naghdi  [83,85]  postulated  an 
entropy  production  inequality  for  the  entire  continuum  to  define 
restrictions  on  the  character  of  various  constitutive  relationships. 
This  was  applied  to  the  case  of  mixture  of  Newtonian  fluids  [83]. 
Crochet  and  Naghdi  [< 5]  applied  these  considerations  to  the  flow  of  a 
fluid  through  an  elastic  solid.  Atkin  [7]  explicitly  stated  the  form 
of  these  constitutive  assumptions  for  fluid  flow  through  an  elastic 
solid.  He  also  presented  an  alternative  method  of  deriving  the 
linearized  theory  of  elastic  solid-viscous  fluid  mixtures  and  the 
thermodynamic  restrictions  imposed  on  this  theory  by  the  entropy 
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production  inequality.  In  later  work.  Green  and  Naghdi  [88]  based  the 
thermodynamic  restrictions  on  the  behavior  of  each  constitutent  on  the 
requirement  that  suitable  combinations  of  the  equations  for  individual 
constituents  should  yield  a  single  entropy  production  inequality  for 
the  mixture  as  a  whole.  Bowen  [28]  noted  that  these  formulations  lead 
to  the  result  that,  in  equilibrium,  the  partial  free  energy  density  of 
a  given  constituent  Is  independent  of  the  deformations  of  the  other 
constituents.  Also  that  such  independence  fails  to  be  confirmed  by 
experiments  on  fluid  mixtures.  Muller  [136]  showed  that  if  gradients 
of  densities  of  the  constituents  were  Included  among  the  constitutive 
variables,  this  difficulty  would  not  arise.  Some  investigators  have 
proposed  use  of  an  entropy  production  inequality  for  each 
constituent.  Bowen  [28]  considers  this  to  be  too  special. 

In  this  section  we  describe  various  phenomenological  approaches 
to  development  of  constitutive  relationships  as  well  as  those  based  on 
theories  of  mixtures. 

4.2.  DIFFUSIVE  RESISTANCE. 


Diffusive  force  was  identified  as  the  interaction  force  by 
Truesdell  [189,190],  Green  [83-88]  and  Crochet  [45]  refer  to  it  as 
'diffusive  resistance'.  Biot  [22]  described  it  as  'disequilibrium 
force'.  The  words  'friction'  and  'drag'  have  also  been  used.  Maxwell 
[118]  defined  it  as  a  pair  of  equal  and  opposite  forces  acting  on  the 
two  constituents  in  a  binary  mixture.  The  force  was  proportional  to 
the  densities  of  the  constituents  and  to  the  relative  velocity.  For 
non-chemical ly  reacting  conti nua,  in  the  absence  of  inertia  effects, 
the  equilibrium  equations  for  the  fluid-saturated  solid  are: 


[tlj(1)  .t,  ‘*>J  ♦Cp(1V1)  ,p(V2)]  =  °  ( IV-6) 


♦ 

where  b. ^ 
constituent 

is  i 

s(k) 

(1) 

C'ij.l111  *  P(1)b,(1)]  -  -  *P(2V2)]  ( IV-7) 
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Each  side  of  the  equality  (IV-7)  representing  interaction  between  the 
constituents  is  set  equal  to  the  diffusive  resistance  Dj.  A  set  of 
single-constituent  stresses  in  equilibrium  can  be  added  to  the 
stresses  on  either  side  without  affecting  the  definition  of  D  j .  For 
instance,  let  ,  k=l,2  such  that  0  •  Then, 

evidently,  addition  of  t..^  to  t..^  in  (IV-7)  will  not  affect 

I  J  '  J 

the  value  of  the  quantities  on  either  side.  However,  superposition  of 

a  total  stress  state  t\ . ^  +  T.  satisfying  equilibrium  will 

■  J  *  J 

influence  Dj.  For  hydrostatic  fluid  stress. 


where 


v-<",*.>(2V2,>-W2) 

(TJ1(1)  +TJ1(2)).j  - 0 


( IV-8) 


( IV-9) 


If  inertia  effects  are  included,  for  a  binary  mixture,  (II 1-116 )  gives 
p1^k^  the  quantity  conjugate  to  the  relative  velocity  in  the  energy 
equality  (II 1-114)  as 

-P,  (IV-10) 

This  definition  is  somewhat  more  general  than  the  one  used  by  Green 
and  Naghdi  [85,86]  where  the  term  involving  c^  was  ignored.  In 
[87],  however,  following  Mills  [122],  a  general  form  of  (IV-10)  is 
stated  for  a  mixture  of  n-1  fluids  and  a  solid. 

Truesdell  [190]  examined  mathematical  theories  of  diffusion  and 
proposed  a  mechanical  theory.  Four  existing  theories  viz. 
kinematical,  hydrodynamical ,  kinetic,  and  thermodynamic  were 
recognized.  The  kinematical  theory  leads  to  Fick's  law.  Based  on 
analogy  to  the  flow  of  heat,  Fick  [64]  proposed  that  in  a  binary 
mixture  at  uniform  total  density 


p<k>  u.<k>  = 
M  l 


( I V— 11) 


where  is  the  density  of  constituent  s^k\  u.^  its  diffusive 

velocity  and  D  is  a  phenomenological  coefficient.  Onsager  [143] 
proposed  a  generalization  of  (IV-11)  taking  p^  u^k^  as  a  linear 


combination  of  the  density  gradients  of  all  constituents.  However,  as 
pointed  out  by  Truesdell,  the  purely  kinematical  presumptions  leading 
to  Pick's  law  are  unsupported  by  any  principle  or  method  of  mechanics 
or  physics  and  cannot  be  expected  to  explain  or  predict  motions  except 
in  specially  simple  circumstances.  Fick’s  law  is  applicable  for  no 
mass  exchange  between  constituents,  constant  pressure  and  total 
density  and  for  no  mean  motion.  Bowen  [29]  pointed  out  that  (IV-11) 
only  represents  a  special  case  of  a  more  general  mixture  theory.  He 
also  showed  that  (IV-11)  and  the  equation  of  mass  continuity,  (III-37) 
lead  to  Pick's  second  law  which  involves  rate  of  change  of  constituent 
density.  Stefan  [182]  generalized  Maxwell's  equation  and  proposed 

pCOf^VKy10-  \((k)-j  A.(<k>  P<kyj>  (u1(k,-u1(J,)  (iv-12) 

j  ~1 

f  kl 

for  a  mixture  of  n  constituents.  Here  uv  ’  is  the  hydrostatic 
pressure  (equal  to  partial  pressure  in  condition  of  equilibrium). 

for  given  i,  are  elements  of  a  matrix  of  phenomenological 
coefficients  assumed  to  be  independent  of  the  densities  of  the 
constituents.  for  k  not  equal  to  j  was  intended  though  not 

so  stated.  The  concept  of  equal  and  opposite  action  implied  A^2)  * 
A£^^.  For  mixtures  with  more  than  two  constituents,  Aji^  = 
was  assumed  without  comment.  This  is  due  to  the  assumption  of 

pairwise  interaction  of  the  constituents.  The  argument  of  equal  and 
opposite  forces  for  a  binary  mixture  does  not  extend  to  three  or  more 
components.  When  more  than  two  constituents  are  present,  the 
interaction  of  any  pair  may  be  affected  by  the  presence  of  the 
others.  Truesdell  [189]  cites  results,  for  the  Chapman-Enskog  process 
in  Maxwell’s  kinetic  theory  obtained  for  multicomponent  mixtures  by 

Hi  rschfelder ,  Curtiss  and  Bird  [94],  Truesdell  showed  that  the 

kinetic  theory  yields  a  quasi-static  approximation  to  the  Maxwell  - 
Stefan  theory,  with  the  addition  of  a  term  representing  thermal 
diffusion.  Also  that  the  coefficients  Aj-j^  depend,  in  general,  upon 
densities  and  molecular  constants  of  all  constituents  and  the  matrix 
Aj^),  for  given  i,  is  not  necessarily  symmetric.  According  to  the 
kinetic  theory,  diffusion  is  a  sum  of  binary  phenomena  in  first 

approximation  but  not  in  a  more  accurate  treatment.  Truesdell  [190] 


noted  that  in  the  kinetic  theory,  accelerations  are  not  considered  and 
hence  it  is  to  be  regarded  as  an  approximation.  The  Maxwel 1 -Stefan 
theory  is  more  general  in  that  accelerations  are  included.  However, 
Stefan's  relations  regarding  coefficients  Aj^*4)  arise  in  first 
approximation  of  the  full  results  at  the  first  stage  of  the  Chapman- 
Enskog  process.  Each  of  the  theories  has  aspects  of  generality 
wanting  in  the  other.  Eckert-Meixner  thermodynamic  theory  implies, 
neglecting  thermal  diffusion. 


pOOu.OO  ."l  L..<k>F. 

’  j-1  J 


( IV-13) 


where  u^k)  is  defined  by  (III-21)  and  Fj  are  'forces'  defined 
differently  by  different  authors.  An  additional  condition,  sometimes 
imposed  is 


l  F, 


00  . 


( IV-14) 


However,  if  no  terms  are  neglected,  there  is  no  reason  to  believe  (IV- 
14)  will  hold.  However,  from  ( 1 1 1 -19)  and  ( 1 11-21) 


Hence 


P<kVk).o 

j-1  1 


( I V- 1 5 ) 


l  l  V  Fj 

j=l  k*l  J  J 


F,  =  0 


( IV-16) 


The  coefficients  L j  0 k ^  are  not  uniquely  defined  but  must  satisfy 


1  Lji 
J-1  J 


^  =  L.  =  constant 


(IV-17) 


Further,  the  Onsager  relations  are  often  assumed  i.e.  LjOO  is 
assumed  to  be  symmetric  in  k,  j.  In  fact,  to  satisfy  the  non-negative 
production  of  entropy  all  that  is  needed  is  that  the  matrix  LjO*0  for¬ 
given  i,  be  non-negative  definite. 

Truesdell  [190]  proposed  a  more  general  mechanical  theory  of 
diffusion.  According  to  Truesdell  [190],  "diffusion,  being  a  change 
of  motion,  arises  from  forces;  the  motions  produced  by  these  forces 
must  conform  to  the  principle  of  linear  momentum  applied  to  each 


constituent  and  to  the  whole  mixture".  Thus,  for  constituent 
the  supply  of  momentum  is  defined  by,  (1 11-48), 


(k)  (k),-  (k)  .  (kk  .  (k) 

pm,'  •  p'  '(f,'  b, '  ')  -  tJ,  j' 


pc 


(k) 


( I V - 18 ) 


The  sum  of  supplies  of  partial  momenta  to  the  constituents  must 
vanish.  The  simplest  constitutive  equation  for  supply  of  momentum 
would  be  a  linear  dependence  of  the  partial  momentum  upon  the  quantity 
‘corresponding1  to  it  in  the  energy  balance  equation  (II 1-51 ) .  viz., 
the  relative  velocity  of  the  constituent  with  respect  to  the  others. 
For  an  isotropic  medium. 


p  m. 


(k)  . 


n 

■  I 

j=l 


*ji 


(k)  (k) 

v  j 


-  v 


(j). 


( I V- 19 ) 


The  scalar  coefficients  are  defined  uniquely  for  k  not  equal  to 

j  and  must  vanish.  For  each  constituent  as  perfect  fluid,  (IV- 

18)  and  (IV-19)  lead  to  Stefan-Maxwell  equations.  The  determination 
of  constitutive  equations  for  materials  has  nothing  to  do  with 
diffusion  in  linear  theories.  The  restriction  that  the  sum  of  partial 
moment  supplies  must  vanish  places  a  restriction  upon  the  coefficients 
Truesdell  showed  that  this  restriction  leads  to  the  necessary 
and  sufficient  condition 


"  (l-J10  -  Lki(j))  =  0  (IV-20) 
j=l 

For  a  binary  mixture,  this  implies  Considering  the 
rate  of  work  of  diffusive  drags,  it  is  necessary  that  diffusion  be  a 
dissipative  process.  Inversion  of  relations  (IV-19)  gives 


(k) 

p  v . 


(k)  . 


n 

•  l 

j  =  l 


ji 


(m 


mi(j)) 


( I V  -  2 1 ) 


where 


l  (Gji 

j=l  J 


(k) 


Si<j)>  ■ 


( I V  -  22 ) 


For  a  mixture  of  perfect  fluids,  t.^  =  Substituting  in  (IV-21) 


and  eliminating  using  the  equation  of  conservation  of  mass  of 

the  constituent  (II 1-45 ) ,  Truesdell  obtained  an  equation  for 

rate  of  change  of  density  of  a  constituent  in  terms  of  gradients  of 


pressures, 


Theories  of  diffusion  involve  variables  other  than 


pressures.  Assuming  the  existence  of  a  thermal  equation  of  state  for 
each  constituent,  Truesdell  showed  that  writing  gradients  of  pressure 
in  terms  of  gradients  of  density  and  temperature,  neglecting 
accelerations,  for  a  binary  mixture,  the  theory  reduces  to  Fick's  law 
under  isothermal  conditions  when  body  forces  are  absent  and  when  the 
total  density  is  uniform.  However,  the  general  form  of  the  theory 
holds  whether  or  not  there  exist  thermal  equations  of  state.  Truesdell 
showed  that  the  kinetic  theory,  the  hydrodynamical  theory  and  the 
thermodynamical  theory  were  all  specializations  of  the  general  theory 
of  mechanical  diffusion  proposed. 

Moxadam  [125]  proposed  constitutive  equations  for  heat  and  fluid 
flux.  These  essentially  amount  to  regarding  the  diffusive  force  to 
have  the  form 


„  a  b  t 

p<  “•7Vl'7  T-’ 


(IV-23) 


Crochet  [45]  showed  (we  give  a  summary  of  the  arguments  later  in 
this  chapter)  that  under  isothermal  conditions  and  in  the  absence  of 

chemical  reactions,  the  constitutive  relations  will  involve  F.-  •  ( ^ )  t 

(k)  J 

d.  .v  v.  and  A.  . .  For  the  linear  theory  of  irrotational  relative 

*  J  *  *  J 

motion,  in  the  absence  of  chemical  reaction  and  inertia  effects,  and 
non-Newtonian  behavior  in  a  binary  mixture  of  a  fluid  and  a  solid  this 
immediately  leads  to  an  expression  of  the  type 


p.  =  u  .  +  p(2)b  (2)  =  _  c  y 


(IV-24) 


The  inverse  form  of  this  equation  is  the  well-known  d'Arcy  flow  rule 
[146,  160].  Biot  [15,  18-25]  and  earlier  investigators  (e.g., 

Terzaghi)  [187]  had  used  this  rule  as  the  starting  point  for  their 
theories.  Biot  [16]  also  assumed  the  existence  of  a  dissipation 
function  so  that  the  permeability  tensor  was  symmetric.  Schiffman 
[166]  extended  d'Arcy's  law  to  non-isothermal  consolidation.  Further 


generalizations  based  on  Onsager's  principle  have  been  proposed  [e.g. 
25,  142,  156,  157,  159,  186].  Generalization  to  nonlinear  cases  has 
been  proposed  [128]  assuming  to  be  functions  of  porosity 

Ghaboussi  and  Wilson  [74]  pointed  out  that  Biot's  [21,  22]  formulation 
for  momentum  balance  may  be  regarded  as  a  generalization  of  (IV-24)  to 
include  inertia  effects  in  the  body  force  term  giving,  for  no  chemical 
reaction  and  isotropic  fluid  stress,  a  specialization  of  (IV-10) 


(IV-25) 


Regarding  (II 1-116 )  as  the  definition  of  p^ ,  a  general  relationship 
based  on  the  correspondence  of  p^k^  and  the  relative  velocity 
v.U)  -  Vi<2)  in  ( 1 1 1-114)  is 

Pi  =  +  P(2)(b.(2)  -  f.(2))  -{p  c(2)Vi  =  -  C^.v.  ( IV-26) 

Morland  [131]  stated  the  total  volume  flux  across  a  fixed  surface 

as 

?  (k)  (k) 

I  nv  'v.v  'n  - 

k-1  1  1 

For  a  binary  mixture  of  fluid  and  solid  this  reduces  to 


( IV-27) 


Thus,  the  volume  flux  of  the  fluid  relative  to  the  solid  surface  is 


This  was  used  as  the  justification  for  writing  d'Arcy's  law  in  terms 
of  relative  flux.  We  note,  however,  that  (IV-27)  does  not  follow  from 
mass  continuity  except  for  incompressible  constitutents.  [cf.  (HI- 
125)]. 


szm 


4.3.  STRESS-STRAIN  RELATIONS. 


4.3.1.  Effective  Stress. 


Terzaghi  [187]  introduced  the  concept  of  effective  stress.  This 
was  defined  to  be  the  stress  component  causing  deformations  of  the 
soil.  For  hydrostatic  pore-water  pressure  acting  on  incompressible 
soil  grains,  the  entire  deformation  of  the  soil  was  assumed  to  be  due 
to  changes  in  the  pore  volume  and  pore  geometry.  For  this  case, 
Terzaghi  [187]  termed  t-jj^  in  (IV-2)  as  the  effective  stress  related 
to  deformation  of  the  solid.  It  was  found  [e.g.  139]  that  the  fluid 
pressure  did  in  fact  influence  the  effective  stress-strain 
relationship  when  the  solid  grains  had  compressibility  comparable  to 
that  of  the  matrix  as  a  whole  and  the  fluid  was  not  incompressible.  To 
allow  for  this  the  effective  stress  related  to  deformation  was  defined 
as 


( IV- 28 ) 


+  (1  -  C)  ix6.  j 


(IV-29) 


where  c  =  1  implies  Terzaghi 's  definition  and  c  =  0  would  correspond 
to  total  stress  being  regarded  as  effective.  Suklje  [183]  discussed 
selection  of  appropriate  values  of  c.  Schiffman  [165]  expected  c  to  be 
between  n  anc|  ^  This  was  based  on  the  assumption  that  the  pore 
fluid  pressure  may  not  act  over  the  entire  area  of  surface  elements 
but  only  over  a  part.  This  fractional  area  of  action  of  the  fluid 
pressure  would  be  bounded  below  by  the  porosity  and  above  by  1.  It  is 
relevant  to  note  here  that  several  investigators  confuse  the  effective 
stress  with  the  partial  stress.  This  is  partly  due  to  the  dual 
definition  given  by  Terzghi.  As  used  in  this  review,  the  effective 
stress  is  the  stress  component  causing  deformations  of  the  solid  and 
is  thus  defined  completely  by  these  deformations.  The  partial  stress 
in  the  solid  cou^d  conceivably  be  related  to  quantities  other  than  the 
deformations  of  the  solid.  Considering  the  solid  grains  to  be 
compressible,  Nur  [139]  derived  the  equation 


-81- 


( I V- 30 ) 


V 


C  =  1  -  K/K. 


where  K,  K$  are,  respectively,  reciprocals  of  the  'bulk'  and  the 
intrinsic  of  the  solid.  For  incompressible  grains  K$  -*■<=  and  for 
highly  deformable  pore  space  K+0.  Hence,  for  either  of  these  cases 
K/K$  =  0,  and  Terzaghi's  definition  is  recovered.  Suklje  [183] 
suggested,  without  proof 


c  =  1  -  n 


(1)  K 


( I V- 3 1 ) 


Schiffman  [165]  gave  a  more  general  form  for  (IV-31)  allowing  fluid 
pressure  to  be  a  second  rank  tensor  and  c  a  fourth  rank  tensor,  i.e., 


t. .  =  t  . 
ij  iJ 


Aklijt 


i j  kl 


( IV-32) 


where  the  superscript  2  refers  to  the  fluid.  Carroll  [34]  carried  out 
a  similar  development.  These  approaches  were  based  on  superposition 
of  effects  of  the  hydrostatic  stress  and  the  stress  deviation. 
Carroll  [34]  determined,  for  the  linear  case,  the  relation 


tij  =  lij  '  EklijCklmn  11  6n 


( IV-33) 


where  E^-j  ^ j  are  components  of  the  elasticity  tensor  for  the  dry  solid 
material  and  C^  j  ^  ^  ^ s  ^  those  of  the  intrinsic  compliance  under 
hydrostatic  stress.  For  clays,  Mitchell  [124]  would  include 
attractive  and  repulsive  forces  within  the  material  in  defining  the 
effective  stress.  Biot  [22]  defined  effective  stress  as 


t, .  =  tid  -  a«*  6,j 


( I V - 34 ) 


and  assumed  t^1)  to  be  the  quantity  related  to  solid  deformation. 
Morland  [130]  criticized  Schiffman's  work  stating:  'Schiffman 

neglects  the  scaling  of  the  partial  solid  stress  to  define  an 
effective  stress  in  the  matrix.  As  a  consequence,  in  his  linear  soil- 
water  model  with  an  interaction  coefficient  given  by  the  water 
porosity  (n^),  for  example),  the  elastic  soil  and  water  responses  are 
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independent,  in  contrast  to  the  linearized  elastic-water  mixture 
equations  derived  later  from  the  present  model  and  to  the  Biot  [1956] 
equations'.  This  criticism  is  based  essentially  on  different  notions 
of  'effective'  stress.  Morland  defined  'effective'  stress  by 

the  equation 


<k>  *„<k> 


( IV-35) 


and  proposed  constitutive  equations  for  the  scaling  factors  m^k^. 
This  is  different  from  the  definition  in  ( IV-28)  above.  Later, 
Morland  [131]  used  the  term  'intrinsic'  to  replace  'effective'  to 
avoid  confusion. 

Garg  [70],  following  Haimson  [92],  proposed  a  dual  definition  for 
effective  stress.  For  strength  of  rock,  he  would  set  c=l  in  (IV-28) 
but  for  constitutive  relations  another  value  of  c  would  be  used. 

Carroll,  Schatz  and  Yamada  [35]  introduced  intrinsic  solid  stress 
on  the  solid  particles  and  an  effective  stress  influencing  deformation 
of  the  pore  space.  Using  the  equt.ion 


JJ 


■n^t 


OV 
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n<2>t 
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JJ 


( IV- 36 ) 


the  solid  stress 


^jj 


n<2>  n)/^ 


( I V- 37 ) 


Thus,  in  general,  t^j(^  would  be  a  function  of  the  volume  change  of 
the  solid  and  t..  =  t,,  -  n  would  be  a  function  of  porosity.  Kenyon 

J  J  J 

[105,106]  also  considered  the  effect  of  grain  and  fluid  compressibil¬ 
ities  and  introduced  material  parameters  to  characterize  this  depend¬ 
ence.  Contact  stress  in  the  solid  ana  tne  duik  stress  independent  of 
K$  were  used. 

For  large  deformation,  as  incremental  form  of  the  stress  tensor 

was  introduced  by  Biot  [24],  Carter  [37]  and  Prevost  [147]  used  the 

Jaumann  stress  rate  Tis  to  ensure  frame  indifference.  Its  relation  to 

J  (11 

the  Cauchy  stress  t..  is 


-83- 


( IV- 38 ) 


T  (1)  _  t  (1)  t  (Dr  (1)  t  (Dp  (1) 
*  '  tik  xkj  "  jk  1  ki 


ij 

(1) 


1J 


where  '  are  components  of  the  antisymmetric  rate  of  rotation 
tensor  for  the  solid  in  Eulerian  description  and  the  superposed 
denotes  the  material  time  deri"*tive.  The  effective  stress  rate 
was 

\j  *  Ti.j  -  *  61j  <IV-39> 


4.3.2.  Pore  Pressure  Parameters 

Skempton  [179,180]  proposed  a  simplification  of  the  theory  of 
consolidation.  He  suggested  that  the  fluid  pressure  be  expressed  in 
terms  of  the  total  stress  components  through  pore  pressure  para¬ 
meters.  A  general  statement  of  this  type  would  be 


K  1  x  k  i  1/2  K  ,  1/3 
n  -  bjjj  +  b2J2  +  b3J3 


( IV-40) 


where  b^  are  material  constants  and  are  the  invariants  of  t^,  the 
total  stress  tensor.  In  order  that  the  formulation  reflect  the 
difference  in  fluid  pressures  observed  during  loading  and  unloading 
cycles,  and  their  steady  increment  with  plastic  straining  and  at 
constant  strain  during  failure,  Lo  [113]  proposed  use  of  strain  invar¬ 
iants  1^  of  e.jj(l)  instead  of  .  The  coefficients  b..  have  been  found 
[100]  to  depend  upon  stress.  However,  both  these  approaches  are 
similar  in  character  and  as  stated  by  Skempton  [179],  are  semi- 
empirical.  Nur  [139]  argued  that  the  strains  caused  by  the  hydro¬ 
static  pore-water  pressure  and  the  total  stress  cancel  each  other 
i  .e. , 


Mij  11  +  Ck1  i  jtkl=  0 


( i v-4i : 


Consequently, 


-1, 


n  ::  “(Mij)  Cklijtkl 


( IV-42 ) 


For  isotropic  behavior,  this  would  be  a  special  case  of  (IV-40). 


4.3.3. 


Stress-Strain  Relations. 


4. 3. 3.1.  Phenomenological  and  micromechanical  approaches. 

Tsien  [193]  proposed  a  linear  elastic  isotropic  relation  for 
tjj^  in  terms  of  using  Terzaghi's  definition  i.e.  c=l  in  (IV- 
28).  Biot  [15]  assumed  a  quadratic  energy  function  in  0,  the  change 
in  water  content  per  unit  volume  of  the  solid,  and  e-jj^  leading,  for 
isotropic  linear  elastic  soil  and  incompressible  fluid,  to 


n  =  M  ekk  +  Ne 


tij  =  2  ^  6ij  +  X  ekk  6ij  +  M6&ij 


( I V -43 ) 


( I V - 44 ) 


In  later  work,  [21],  the  total  stress  in  (IV-44)  was  replaced  by  the 
solid  stress  t..^.  In  extension  to  anisotropic  elastic  [16]  solids 

1  J 

and  compressible  fluid  the  relationship  was  stated  as 


t  (1)=  C  e  (1)  +  M  e  (2) 
S'j  Lk  1  i j  ekl  +  Mij  ekk 


/  =  M..  ei.(1)  +  M  ekk' 


( IV-45) 


( I V -46 ) 


Similar  construction  was  used  for  viscoelastic  soils  [16].  In  [24] 
the  quantity  ekk^  was  again  replaced  by  9  as  used  in  ( IV-43)  and 
(IV-44).  The  same  concept  was  extended  to  the  case  of  finite  elastic 
deformation  [25].  Biot  [22]  defined  effective  stress  through  ( I V - 3 1 ) 
but  the  quantity  related  to  solid  deformations  was  taken  to  be  t^j^^ 
defined  by  ( I V- 3 ) .  This  leads  to 


TjU)  *  +  (l-‘>  '*  6*j  ”  E*,ij  eki(U 


( I V  -  4  7 ) 


Considering  u  to  be  given  by  (IV-46) 


Tj  ■  EXlij  eK,(l)  -  («-»  C«pq  'pq,’)  -  «ij  ( IV-48) 


,V»»P.WiTi 


-  •  ’*  •  V r  »»  '•*-*- 


A0-A174  895 
UNCLASSIFIED 


MECHANICAL  BEHAVIOR  OF  SATURATED  SOILS  -  A  REVIEM(U)  2/2 

OHIO  STATE  UN IV  RESEARCH  FOUNDATION  COLUMBUS 
R  S  SANOHU  JUN  85  0SURF-715927-85-1  AFOSR-TR-8C-2095 
AFOSR-83-0055  F/G  8/13  NL 


2) 

For  incompressible  fluid  and  incompressible  solid  matrix,  e,.v  '  =  0 

J 


a  equals  one.  Hence  t^1)  and  t..  coincide* 

*  J  1  J 


Lubinski  [114]  assumed  that  the  total  strain  of  a  bulk  porous 
solid  can  be  expressed  as  a  summation  of  the  strains  due  to  pore-water 
pressure  and  strains  due  to  stress  acting  on  the  solid  skeleton.  He 
proposed  relations  of  the  type 


t  (D  .  r  .  (1)  +  /n(D 
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where  y,  M,  N,  are  constants  to  be  determined  experimentally.  Krause 
[109]  added  terms  to  the  right  side  of  (IV-43)  to  reflect  linear 
dependence  of  the  fluid  pressure  on  the  deformation  rate  d.^2)  of  the 
fluid.  This  assumes  a  viscous  component  for  fluid  flow.  Adkins,  in 
his  earlier  theory  [2],  assumed  that  the  stress  in  each  constituent 
depended  only  on  the  density  and  the  kinematic  quantities  associated 
with  only  that  constituent.  Nur  [139]  assumed  the  effective  stress, 
given  by  (IV-28)  along  with  (IV-30),  to  be  related  to  e- j .  This 
admitted  a  certain  dependence  of  t-jj^  upon  the  fluid  pressure. 
Expl icitly , 


Hence 
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( IV-52) 


These  approaches  were  based  on  the  superposition  of  effects  of  the 
hydrostatic  stress  and  the  shear  stress.  Carroll  [34]  determined, 
that  for  the  linear  case,  (IV-52)  would  be 
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where  are  components  of  the  elasticity  tensor  for  the  dry  solid 
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material  and  C. '  those  of  the  intrinsic  compliance  under 
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hydrostatic  stress.  For  intrinsic  hydrostatic  fluid  stress  %  acting 
on  the  solid  particles,  the  strain  in  the  solid  phase  is 


material  and 
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( IV-54) 


The  strain  in  the  pore  space  is  due  to  the  bulk  solid  behavior,  i.e., 


eij(f)  =  Ckl i j  (tkl  "  11  6kl^ 


( IV-55) 


Therefore,  the  total  strain,  adding  (IV-54)  and  (IV-55),  is 
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Hence  (IV-53).  This  formulation  was  specialized  to  allow  for  the 
presence  of  internal  symmetries.  For  isotropy,  the  formulation 
reduces  to  Nur  and  Byerlee's  [IV-29].  Schiffman  [165]  proposed  a 
generalization  of  (IV-52)  in  the  form 


tij^1^  =  Ekl i j ei j ^ 1  ^  ’  ^6ik  6jl"  \lij^  Hl^  (IV-58) 


The  quantity  Aj jki  was  termed  the  soil  water  interaction  tensor.  Garg 
[70]  obtained  a  relationship  between  the  intrinsic  and  the  bulk 
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Dehavior  of  rocks  under  hydrostatic  stress.  Interpretation  of 
coefficients  appearing  in  his  equations  were  attempted  by  Biot  [20], 
Additions  to  this  discussion  have  appeared  from  time  to  time  [e.g. 
105-106,139,150]. 

In  extending  the  theory  to  the  nonlinear  case,  Westmann  [196] 
assumed  t^-^J  to  be  a  function  of  the  Cauchy  deformation  tensor  for 
the  solid  and  the  rate  of  deformation  (Eulerian  description)  of  the 
fluid.  The  fluid  stress  was  expected  to  consist  of  a  hydrostatic 
component  and  another  component  depending  upon  the  same  quantities  as 
t.jj(*).  It  was  noted  that  in  this  formulation  it  would  be  difficult 
to  design  experiments  to  evaluate  the  parameters.  A  simplification 
proposed  assumed  the  fluid  pressure  to  be  hydrostatic  and  related  to 
the  velocity  field  through  Darcy's  law.  This  is  similar  to  Sandhu's 
[156-158]  argument  that  the  constitutive  equation  for  diffusive 
resistance  is  a  sufficient  relationship  between  the  fluid  partial 
stress  and  kinematics  of  the  mixture.  Westmann  [196]  wrote  relative 
velocity  as  a  function  of  n,  D.,  and  the  Cauchy  deformation  tensor 

J 

for  the  solid.  This  would  reflect,  among  other  factors,  the 
dependence  of  permeability  on  the  porosity  of  the  solid. 

Bedford  and  Drumheller  [14]  have  pointed  out  that  fluid-saturated 
porous  media  fall  in  the  class  of  immiscible  mixtures.  The 

constituents  of  such  mixtures  remain  physically  separate  on  a  scale 
which  is  large  in  comparison  with  molecular  dimensions.  This 
immiscibility  has  two  important  consequences.  Because  of  physical 
separation,  in  some  local  sense,  each  constituent  will  obey  the 
constitutive  relations  for  that  constituent  alone.  Also,  the 

constituents  intrinsically  have  microstructure  defined  by  the 

interfaces  which  separate  the  constiuents.  To  set  up  macroscopic 
constitutive  relations,  one  approach  would  be  to  postulate  these 
relations  directly  as  described  above.  Another  alternative  would  be 
to  relate  macroscopic  behavior  to  intrinsic  properties  of  the 
constituents.  The  simplest  theories  Involve  the  volume  fraction  of 
the  constituents  in  addition  to  the  usual  variables.  Morland  [130] 

pointed  out  the  meaning  of  deformation  and  stress  associated  with  the 
continuum  model  of  each  constituent.  In  particular,  the  partial 
density  variation  is  not  the  density  variation  of  the  constituent 


since  the  mixture  postulate  eliminates  reference  to  the  actual  volume 
occupied  by  each  constituent  in  an  immiscible  mixture.  Morland  [130] 
and  Garg  [70]  assumed  thr  'effective*  (intrinsic)  stress  in  each 
component  to  be  a  function  of  the  deformation  of  that  constituent  only 
and  having  the  same  form  as  for  a  single  material.  Thus 


f(e  (k)*) 
mm 


( IV— 59 ) 


The  bulk  stresses  and  deformations  were  expected  to  be  related  to  the 
corresponding  intrinsic  quantities  through  scaling  functions.  Thus, 
the  bulk  stress  in  the  constituent  s^k^ 
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For  linear  isotropic  elastic  rock 
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where  KQk denoted  the  bulk  modulus  of  dry  rock.  The  bulk 
deformation  gradient  was  related  to  the  effective  deformation  gradient 

Fi j*  as 

F.j.*  =  [n(1)/n0{1)]1/3Fij(1)  (IV-62) 

The  relation  between  and  t^*)  was  expected  to  be  the  same 

function  f  as  in  (IV-59).  For  an  ideal  fluid,  Morland  [130]  used 

( IV -63 ) 
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where  p  is  the  fluid  pressure,  to  write  partial  stress 

‘tj'2’*  •"(2>cSy(det(Fpq(2)))]  p6ij 
0 

For  the  solid  constituent,  he  assumed  the  isotropic  relation 


(IV-64) 
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where 
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and  a.j  are  functions  of  the  invariants  of  .  The  corresponding 
partial  stress  tensor  was 

t  ( l )_  ( ^ ) rb  s  +  b  r  n(  i3  b  b  r  i3  p  ( i ) r  ( i ) -i 
*ij  n  [b0  6ij+  bl[  -TTT]  Bij  b2[  7TTT-1  Bik  Bjk  ] 

n0  n0  ( IV-67 ) 

Further,  it  was  assumed  that  the  volume  fraction  n^  is  a  dependent 
variable  for  which  a  constitutive  equation  is  required  to  be 
postulated.  For  fluid  stress,  assuming  ideal  fluid,  Morland  [130] 
proposed 
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( IV-68) 


For  a  linearly  viscous  compressible  fluid  with  zero  bulk  viscosity  and 
a  temperature  dependent  shear  viscosity  n(T),  Morland  [131]  wrote 

'ij'2’*  =  '  "*  6ij  +  2  “(2)  <Dfj(2)  '  T  vm,m(2)^f j >  <IV-69> 

Morland  [131]  assumed  the  stress  deviation  in  the  solid  to  be 
dependent  only  upon  its  strain  deviation.  Thus,  for  linearly  elastic 
isotropic  rock, 

(1)*-V-m*/3!’2>‘n(1)  (eij(1)-  6u  eJ1)/3>  (IV-70) 

The  stress-strain  relations  for  the  'bulk'  stresses  were  written  as 
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where  Kv‘;  is  the  compressibility  of  the  constituent  s'* 1 ' .  The 
isotropic  pressures  in  the  fluid  and  the  solid  were  assumed  to  depend 
upon  the  volumetric  strain  of  both  the  constituents,  i.e. 


t . .  =  a  e  ^  +  b  e  ' 
jj  mm  mm 
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(IV-73) 


( IV-74) 


However,  unlike  Biot,  the  existence  of  an  energy  function  was  not 
postulated  so  that  the  constants  b  and  c  in  (IV-73)  and  (IV-74)  do  not 
have  to  be  equal.  Morland  expressed  the  coefficients  in  terms  of 
compressibilities  of  the  fluid  and  the  solid,  and  the  bulk  shear 
modulus  of  the  mixture. 

Carroll  and  Katsube  [36]  postulated  the  following  relations  for  a 
fluid-saturated  solid: 

i.  Relation  between  stress  in  the  mixture  and  stresses  in  the 
constituents : 


1  t  .1  „(l)t  (1)  t  n(2)  * 
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ii.  Solid  stress-strain  law: 


(IV-75) 
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where  the  symbol  A  indicates  change  in  the  quantity  following  it, 


iii.  Effective  stress-strain  law; 
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Combining  these  relationships  they  obtained  the  bulk  relations  for  the 


mixture  as 
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Here  K'  '  is  the  bulk  modulus  of  the  fluid. 

4. 3. 3. 2.  Thermodynamic  considerations. 

Adkins  [3]  and  Green  [81]  admitted  interdependence  of  stress  of 
each  constituent  upon  the  kinematics  of  all.  This  was  in  line  with 
the  principle  of  equipresence  stated  by  Truesdell  [189].  In 
application  to  elastic  materials,  the  existence  of  an  energy  function 
for  the  mixture  was  assumed  by  Biot  [15,18-20,23-25],  This  has  been 
consistently  followed  by  numerous  investigators  [e.g.  7,8,32,45,51,81- 
88,177],  Sandhu  [158,  159]  pointed  out  that  as  the  mixture  could  not 
be  regarded  as  a  continuum  in  motion,  it  was  inappropriate  to  assume 
energy  functions  for  it  in  the  form  that  has  been  popular. 

Morland  [130-132]  did  not  assume  the  existence  of  an  energy 
function  for  the  mixture  but  still  admitted  interdependence.  This 
gave  relations  of  the  type  postulated  by  Biot  (e.g.,  equations  (IV-45) 
and  (IV-46)  with  i  in  the  second  equation  replaced  by  say  H.  .  not 
necessarily  equal  to  in  the  first  equation.  This  implies  a 

nonsymmetric  constitutive  relationship  of  the  type  proposed  earlier  by 
Schiffman  [165].  For  c  =  n^  the  constitutive  equations  for  t.^1) 

and  u  become  uncoupled.  Sandhu  [156-159],  Westmann  [196]  and  Morland 
[130-132]  have  followed  Adkins*  [2]  original  idea  that  the  stresses  in 
each  constituent  depend  upon  the  kinematics  of  only  that 
constituent.  However,  Morland  [130-132]  used  this  for  the  intrinsic 
rather  than  the  bulk  stresses.  This  brings  back  some  dependence  of 
upon  the  fluid  pressure  because  the  porosity  was  postulated  to 
be  a  linear  function  of  the  partial  stresses. 


i.  Entropy  Production  Inequality. 


Mokadam  [126,127]  wrote  the  entropy  production  rate,  per  unit 
volume,  of  a  multicomponent  fluid  as 
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( IV-82) 


Here  m^k^  is  the  density  of  diffusive  flow  of  the  constituent  s^k^ 
and  is  given  by 


(k)  (k)  (k) 
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( IV-83) 


Fi  is  the  total  interaction  force  between  the  porous  solid  and  the 

(kl 

multicomponent  fluid,  and  '  is  the  chemical  potential  of  the 
constituent  s^k^.  The  multicomponent  fluid  and  the  porous  solid  were 
assumed  to  be  at  the  same  temperature.  Irreversible  changes  occur  in 
the  system  due  to  viscous  flow,  heat  conduction,  diffusion  and 
movement  of  center  of  mass  of  the  liquid  phase  in  the  presence  of 
forces  of  friction.  The  first  term  on  the  right  side  of  the  equality 
(IV-82)  is  the  scalar  work  of  terms  of  second  order  while  the  last 
three  terms  are  scalar  works  of  tensors  of  the  first  order.  Hence, 
there  is  no  relation  between  the  first  and  the  other  three  terms. 
This  argument  of  separability,  of  the  first  and  the  last  three  terms 
on  the  right  side  of  (IV-82),  was  the  basis  for  Mokadam's  generalized 
d'Arcy  flow  law  ( 1 1 -5 ) .  Considerations  of  mass  conservation  and  non¬ 
negative  entropy  production  rates  led  to  certain  restrictions  on  the 
coefficients  appearing  in  ( 1 1 -5 ) .  (IV-82)  was  derived  from  Gibb's 
equation.  Truesdell  [190]  objects  to  the  use  of  this  equation  which 
essentially  is  the  condition  that  the  thermodynamic  state  be 
homogeneous.  It  is  questionable  to  apply  it  to  any  problem  involving 
departures  from  equilibrium. 

Green  and  Naghdi  [83,85]  postulated  an  entropy  production 
inequality  for  the  entire  continuum  to  define  restrictions  on  the 
character  of  various  constitutive  relationships.  Crochet  and  Naghdi 
[45]  applied  these  considerations  to  the  flow  of  a  fluid  through  an 


elastic  solid.  The  starting  point  is  the  requirement  that  the 
Clausius-Duhem  inequality  be  satisfied,  i.e.,  [83] 


/  P  (Sr  -  t)  dV  +  /  ydA  >  0 
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( IV-84) 


Assuming  q  =  hj  n  • ,  (II 1-105 ) ,  and  t^  =  tj^nj,  (II 1-99 )  •  the  above 
inequality  upon  use  of  the  divergence  theorem,  gives 


J  Cp  (T5F  '  T 1  +  dV  >  0 


( IV-85) 


Here  S  is  the  entropy  per  unit  mass  and  T  (>0)  is  the  temperature  of 
the  mixture.  Seeking  constitutive  equations  for  internal  energy, 
entropy,  heat  flux  vector,  partial  stresses  and  diffusive  resistance 
vector,  the  point  form  of  (IV-85)  was  written  as,  [84], 
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Substituting  for  r  from  the  energy  equality  (II 1-114) ,  an  alternative 
form  of  the  entropy  production  inequality  is, 
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Shi  [177],  considering  the  binary  mixture  of  a  liquid  and  a  solid. 


added  the  term 


(p(k)u.(k)s(k): 


( I V - 88 ) 


to  the  left  hand  side  of  (IV-86).  Here  Uj(k)  is  the  velocity  of  the 
constituent  s^k^  relative  to  the  barycentric  velocity  v^  of  the 
mixture.  This  gave  the  local  form 
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Dunwoody  and  Muller  [63]  used  the  inequality 
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Here  H^k)  is  a  constitutive  variable.  For  and  T^k^  =  T  for 

all  k,  the  local  form  of  (IV-90)  is  (IV-89).  Truesdell  [191] 
proposed,  for  the  constituent  s^k^,  the  inequality 
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Bowen  [29]  showed  that  (IV-92)  and  (IV-90)  are  the  same  if 
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Green  and  Naghdi  [85]  wrote  the  entropy  flux  balance  equation  for  a 
fixed  volume  bounded  by  surface  A  as 

|p  /  P(k)S(k)dV  -  /  [p[k)  +  pc(k)s(k)]  dV 
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where  S^k)  is  the  entropy  of  the  constituent  s^k^  per  unit  mass'  The 


entropy  production  inequality  was  assumed  to  be 
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where  4>x  ,  zv  '  represent  volume  and  surface  interaction  terms.  The 

local  form  of  the  inequality,  using  results  of  [84],  was  shown  to  be 
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where 


q(k )  +  z(k)  =  nj  (n^k^  +  g^k^ ) 


( IV-97 ) 


( I V -98 ) 


Further,  combining  the  inequalities  for  the  two  constituents  in  a 
binary  mixture,  (IV-91)  was  established.  The  entropy  of  the  mixture 
was  related  to  that  of  the  constituents  by  the  relation 

pT§=  \  [p(k)T(k)  ♦  pc(k)T,k)S(k)  -  5(k)]  CIV-99) 

k-1 

where  is  related  to  </k^,  g-^,  h-^,  hi  and  p^.  In  [86], 

J  J  J  * 

the  Helmholtz  free  energy  was  introduced  as 


A  =  U  -  TS 

Then 

A  =  U  -  TS  -  TS  (IV-100) 

(IV- 100 )  was  substituted  in  the  energy  equality  to  reduce  the  entropy 
inequality  (IV-87),  for  a  binary  mixture,  to  the  form 


-  P  (A  +  ST)-pi(vi(1)-vi{2))+^  tji(k)vijj(k)-  |  h.T'.  >  0  (IV-101) 

i.  Constitutive  Equations. 

Seeking  constitutive  equations  for  internal  energy,  entropy,  heat 
flux  vector,  partial  stresses  and  diffusive  resistance  vectors  p^, 
Q.  and  the  quantity  q  -  hini  in  the  case  of  mixture  of  two  Newtonian 
fluids.  Green  and  Naghdi  [83]  assumed  these  to  depend  upon  the 
densities  of  the  constituents,  the  velocities,  the  gradients  of 
velocities,  and  the  temperature.  For  the  heat  flux  vector,  the 
temperature  gradient  replaced  the  deformation  gradients.  For  the  case 
of  fluid  flow  through  an  elastic  solid,  Crochet  and  Naghdi  [45] 
started  by  admitting  fairly  general  constitutive  assumptions  in  line 
with  Truesdell's  [189]  principle  of  equipresence,  and  then  determined 
the  restrictions  placed  upon  these  general  constitutive  relations  by 
thermodynamic  considerations .  This  approach  is  similar  to  that  used 
by  Noll  [137]  and  Coleman  and  Mizel  [44].  The  following  functional 
dependence  was  assumed. 


U  =  U(T’  Fmn(k)  ’  dmn(k)  ’O 
mn  mn  m 

( IV- 102 ) 

S  =  S(T,  F  (k),  d  (k',  v  ) 
v  mn  mn  nr 

( IV- 103 ) 

tr  =  tr  -/k)(T  F  (k^,  d  (k),  v  ) 

(ij)  (ij)  ,  mn  mn  m' 

( IV - 104 ) 

tr.  =  tr.  . -i ( T,  F  (k),  d  (k),  v  ) 

[ij]  [ij]  mn  mn  m; 

( I V  - 1 0  5 ) 

p.  =  p.  (T,  F  (k) ,  d  (k) ,  v  ) 
i  '  mn  mn  m 

( IV- 106 ) 

q.  =  q.  (T,  F  (k) ,  d  (k) ,  v  ,  T  ) 
m  mn  mn  m  ,m 

( IV- 107 ) 

Here  U,  S,  p^ ,  q^  are,  respectively,  the  internal  energy,  the  entropy, 
components  of  the  diffusive  resistance  vector  and  the  heat  flux 
vector.  This  interaction  force  is  equal  and  opposite  for  the  two 
constituents  in  view  of  the  symmetry  of  the  total  stress  tensor  for 
the  mixtur^  in  the  absence  of  body  couples.  T  is  the  temperature 
f^eld,  Fmp  and  dmn(k)  are  C0mp0nents  of  the  deformation  gradient 


and  the  deformation  rate  defined  in  section  III,  and  are  components 
of  relative  velocity  given  by 


V.  -  V.(2)  -  v  (!) 
l  i  i 


( I V - 108 ) 


Assuming  U  to  be  smoothly  invertible  in  T,  the  variable  T  in  ( IV- 102 ) 
to  (IV- 107)  can  be  replaced  by  U.  Crochet  and  Naghdi  [45]  showed  that 
the  restriction  of  nonnegative  entropy  production  require  that  S  and  U 
be  independent  of  the  deformation  rates,  relative  velocity  and  the 
temperature  gradient.  Thus 


S  =  S(U,  Fi(j(k)) 


( I V - 109 ) 


and  as  U  =  U(T,Fj  _-(k)) 

T  =  T(U,Fij^k^)  (IV-110) 

Defining  t( , j j (0) (k 5 ,  t[ij](0)(k) ,  pi(0)  as  values  of  t(ij)(k>, 

t[ij]^^’  and  pi  ^or  the  e(luOibrium  state  of  no  relative  motion  and 

vanishing  rates  of  deformation,  i.e.,  d-jj^  =  °*  =  0.  and  assuming 

fkl  m 

pj  continuous  in  v^  at  v^  =  0  and  t  ^  ^ v  ',  t^jj  ,  p^  all 
continuous  in  d-.^at  d..^  =  0,  Crochet  and  Naghdi  [38]  also 

*  J  *  J 

established  the  following  relationships: 


t  00  _  1  F  (k)/  &A 

t(ij)(0)  “2phim 

jm 


(IV-111) 


t[ij](0) 


1  r  (1)/  aA  3  A  \ 

■  *  P  ™  ‘rrm  -  ^-TTT> 

jm  mj 


F.  <Z>  (-5* 
jm 


( I V- 1 1 2 ) 


D  _  (D  3A  (2)  .(2)  3A  (1)  ,  . 

Pun\  -  P  - 7TV  Ln  1  -  p  - nr  1  (  IV-113  ) 
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where  the  Helmholtz  free  energy  A  is  a  function  only  of  the 
temperature  T  and  the  deformation  gradients  F^1^  because  U  and  S  are 
dependent  only  upon  T  and  Equations  (IV-110)  thru  (IV-112) 

are  quite  general  and  can  be  specialized  to  binary  mixtures  of  various 
classes  of  materials. 

4. 3. 3. 3.  Fluid  flow  through  an  elastic  solid 

Noting  that  the  rate  of  deformation  d^j^  is  uniquely  related  to 

F^(k)  and  the  spatial  velocity  gradient  j ^ k ^ ,  and  Ln- j ^ k ^  can  in 
turn  be  expressed  in  terms  of  D-.(k)  and  r.  . the  variables  d^^) 
can  be  replaced  by  and  Superposed  rigid  body  angular 

velocity  Q  would  change  r^. ^  to  +  Q  leaving  the  other  vari¬ 

ables  unchanged.  Thus,  for  material  objectivity  to  be  satisfied,  the 
vorticity  tensors  must  occur  as  the  difference  ^.=1^.^  -r..^. 

For  the  fluid,  the  deformation  can  be  represented  by  its 
density.  For  the  solid,  we  can  replace  F^1)  by  the  strain  measure, 
(II 1-27 ) ,  as  the  independent  variable  in  (IV-102)  through  (IV- 

107).  In  the  light  of  the  foregoing,  and  in  view  of  the  thermodynami c 
restrictions  on  the  energy  and  the  entropy  functions.  Crochet  and 
Naghdi  [45]  reduced  the  general  relationships  admitted  by  (IV-102)  to 
( IV— 107 )  to  the  following  form  for  application  to  fluid  flow  through  a 
solid. 


A  =  A(T,p(2),  e..) 

( IV- 1 14 ) 

u  j) 

-  f  (|<)  (J  _(2)  g  n  (1)  p  (2)  \ 

’  u’  p  ’  emn ’  mn  ’  mn  •  Amn’V 

( I V- 1 1 5 ) 

pi 

=  Pi  (T,  p^2),  emn»Dmn^^  ,Dmn^2^  ,Amn’Vm) 

( I V- 1 1 6 ) 

t[i  j] 

=  tr.  (T,  p^  ,  e  ,D  W  ,D  (2) ,  A  , v  ) 

L  t  J  J  mn  mn  mn  mn  m ' 

( IV-117) 

<*1 

=  a  (T  d^2^  e  D  ^  D  ^  A  v  T  ) 

°*i  '  *  p  *  mn’  mn  ,umn  ’  mn,vm’  ,m 

( I V- 1 18 ) 

For  the  equilibrium  state,  i.e., 


vanishing  Di j(k),  A^,  wi  ,Tj  , 


keeping  in  view  the  restriction  that  for  p'  0  the  constitutive 
relation  must  reduce  to  that  for  P-jj^  and  for  p^=  0,  P-jj^  must 
reduce  to  the  form  for  the  solid  only.  Crochet  and  Naghdi  [45]  wrote 


+  ( 1 ) _  1  / _a  c  \  /8A  ^  5A  ^ 

'(«)(*>  * p  '  ^ 


t  (2), 


(2)  &A  . 

PP 


-  A 1)  aA  „  (2)  12)  5A 


aejm  jm,i 


t[ij](0)=  0 


(IV-119) 


( I V- 1 20 ) 


( IV- 121 ) 


( IV- 122 ) 


To  allow  for  static  equal  and  opposite  stresses  in  the  consti¬ 
tuents  of  a  binary  mixture,  it  is  only  necessary,  [86],  to  add  a 

functional  F(p^,  T)(v.^-  v.  *  to  the  Helmholtz  free  energy. 

fki  1 

This  will  add  F(pv  ;  ,T)  to  t .  ...v  '  and  subtract  the  same  amount 

m\  •  J  / 

;  .  The  diffusive  resistance  p^  is  changed  by  Fi#  This 

additional  term  in  the  Helmholtz  free  energy  is  necessary  only  if 

equations  for  partial  stresses  and  diffusive  resistance  are 

required.  Kenyon  [105]  Introduced,  for  a  binary  mixture  of  a  liquid 

and  a  solid. 


K.  =  l  Ch  <k>  -  TX  <k>  -p'k>A<k>„  <">] 
J  k=l  1  1 


( IV- 123 ) 


where  h^k),  are,  respectively,  the  components  of  the  heat  flux 

and  the  entropy  flux  and  A^k)  is  the  Helmholtz  free  energy  associated 
with  the  constituent  Following  Bowen  [30],  Kenyon  assumed 

h,(k)=T  so  that 


K.  >  f  p'k>A<k>u <k> 
J  1--1 


( IV- 124 ) 


The  additional  term  in  the  expression  for  stress  was  in  ti{(^  and 

m  5wi  J 

an  equal  and  opposite  term  in  t^fc>.  It  is  easy  to  show  that 


PU^  =  P  (v.^2)  -  vi)  =  p^  wi 


(IV-125) 


IWWSVOW 


PH,'4'  -  p  (Vi'1'  -  Vi)  -p'1'  w. 


and  pK,  =p(‘>  p<2>  wi(a(2>  -  *(»)) 

.  Further,  setting 


Here  w^  =  v. 


l  r  r  1 

(2)  .  M  (1) 

2 


(IV-126) 
( IV- 127 ) 


'  Vi 


A-  |  p<k>A(k) 
k=l 

Kenyon  [105]  wrote 

.  (2)  .  (2)  r  (2)  8A(2)  .  (1)  SA(1),  , 

*U(0)  '-p  [(J 

Here  is  the  ‘bulk  area-averaged1  fluid  stress, 

intrinsic  hydrostatic  fluid  stress  n  by  the  equation 

n(2)  *  (2)  f  (2)  dA<2)  ‘  (1)  BA*1* 

"  s  Cp  ^  P  ^T1 

Kenyon  [105]  wrote  the  chemical  potential  as 


( IV- 128) 


( IV- 129) 
Defining  the 

( IV— 1 30 ) 


( I V- 1 31 ) 


and  the  balance  equation  for  the  fluid  as 


p(\i-p(2>  b,<«.o 


The  stress  in  the  solid 


(1)  (1)  * 
t..v  '  =  contact  stress  +  nv  '  n 
U 


( IV-132) 


(2)  aA(2) 

p 

*  -U)  bA11’ 

p 

* 

It 

dp(2)  j 

5Fik 

jk 

P  BF.. 

ik 

jk 

a  j  J 

(2)  SA(2) 

r 

.  (l)  aA(1) 

c 

„d) 

★ 

Fjk 

p 

Fjk 

-  n 

*  6-j  - 

J  6ij 


where 
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1  *  i  j 
( IV-133) 

( I V- 1 34 ) 


-V., 


are  functions  of  the  intrinsic 


For  negligible  adsorption,  A^),  p^* 

★ 

fluid  pressure  n  only.  Hence 


t. . 
U 


(1)  * 


AD  aA 
p  aTT 


(1) 


ik 


,(1) 


ij 


* 

Tt  K  6 


lj 


( I V- 1 35 ) 


For  incompressible  mixture  <=0  and  the  above  equation  simplifies  to 


(1) 


ij 


=  P(1) 


a  A' 


(i) 


bf. 


ik 


jk 


-  n 


(1) 


ij 


( IV-136) 


4. 3. 3. 4.  Linear  theory 

In  the  development  so  far,  apart  from  the  functional  dependence  of 
A,  t^  t[ij](k^,  and  upon  thermal  and  kinematic 

variables,  there  is  no  restriction  as  to  the  form  of  the  relationship 
admissible.  Thus,  the  relationship  may  be  linear  or  nonlinear  and  the 
constituents  may  be  rate-type  materials  or  materials  with  memory. 
Smoothness  and  invertibi  1  ity  of  U  and  S  have  been  assumed  as  also  the 
admissibility  of  e^j  as  a  measure  of  deformation.  Saturated  state, 
i.e.,  all  points  of  the  system  being  occupied  by  both  the  constituents 
is  another  assumption  made  throughout  the  theory  of  mixtures.  The 
relationships  for  the  linear  theory  can  be  directly  obtained  as  a 
specialization. 

For  infinitesimal  deformation,  the  strain  measure  for  the  solid 

constituent  is,  (II 1-26 ) ,  e- •  =  uM  For  the  relationships  for 

m  rn  J  ,J 

t(ij)  ’  t[ij]L  »  and  pi  be  '’1’near'  the  Helmholtz  free  energy  has  to 
be  quadratic  in  T  and  the  measures  of  strain  in  the  constituents. 
Thus  for  an  initially  isotropic  solid,  in  the  absence  of  body  couples 
(i.e.,  symmetric  total  stress),  the  'equilibrium'  values  of  the 
mechanical  variables  are  given  by  the  following  equations: 


‘(UHO)'1’  ■  X1  *  ^1  eij  *  *  >(Z)  6ij  -  «  T  6ij  (IV'137: 


t, .  =  ( -K  p(2)  +  Y  e  -  p  T)  6-  . 

(ij)(0)  y  H  '  m  r  >  -ij 


Pi  =  0 


( IV- 138) 
( I V- 1 39 ) 
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0 


( IV- 140 ) 


[ij](0)  =  u 

where  X^,  n^t  y,  and  K  are  constants  depending  upon  initial  conditions 
and  temperature.  Symmetry  of  the  constitutive  relations  for  the 
stresses  follows  from  the  existence  of  the  Helmholtz  free  energy 
function . 

It  may  be  noted  here  that  the  equations  (IV- 137 )  and  (IV- 138 )  are 

essentially  the  same  as  those  postulated  by  Biot  [15]  using  fluid 

strain  as  the  additional  kinematic  parameter.  The  conceptual 

difficulties  associated  with  Biot's  theories  are  now  obvious.  Biot 

regarded  the  irreversible  flow-deformation  process  as  a  series  of 

equilibrium  states,  independent  of  all  rate  parameters.  This  would 

conflict  with  the  constitutive  equation  of  flow  (d'Arcy's  law)  which 

relates  the  partial  stress  in  the  fluid  to  the  relative  motion  rate  of 

the  solid  and  the  fluid.  In  the  theory  proposed  by  Crochet  and  Naghdi 

[45],  the  mechanical  component  may  depend  upon  all  kinematic 

variables.  For  the  linear  theory  this  dependence  becomes  additive. 

The  total  values  of  these  quantities  would  consist  of  an  equilibrium 

and  a  non-equilibrium  component.  For  the  'excess'  stresses  and  other 

quantities  in  addition  to  the  'equilibrium'  values,  we  would  expect 

(k) 

dependence  only  upon  D.  ,  A..,  and  v-,  which  vanish  at 

*  J  •  J  J 

equilibrium.  Thus,  for  isotropic  materials, 

t, .  .  J^=t,.  .un^^+  X-D  ^6.  .  +  2^-0.  MK  X.D  ^6.  -+2|i.D. 

(ij)  ( i J ) ( 0 )  3  mm  ij  ^3  t j  4  mm  ij  ^4  ij 

( I V- 14 1 ) 

tM -/2)=tHiWn/2)+  M  ^6- •  +  2pKD.  X9D  ^2)S. -+2u9D 


5  mm  uij  ^5  i j 


2  mm  i  j  ^2  i  j 


Pi  =  C3  wi  +  c4  eijk  Ajk 


'[ij]  =  C1  Aij  +  C2  eijkWk 


q.  =  -KT  .  +  c,  v.  +  cc  e.  ..  A., 
l  ,i  5  l  6  ijk  jk 
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( I V - 1 42 ) 


( IV-143) 


( I V - 144 ) 


( IV- 145 ) 


1 


H.J'* 

■«>l 


vis 

•A 

I 


■m 


f 

I 


where  e^jk  is  the  permutation  symbol.  For  initially  isotropic 
materials,  c2,  c4,  Cg  vanish.  To  satisfy  the  Clausius-Duhem 

inequality  for  arbitrary  velocities  and  velocity  gradients,  the 
following  must  hold  [45]: 


Cj,  Cy  K,  must  be  non-negative 

\2  +  |  P2)  >  ° 


(X3  +3  n3)  >  0 


(p4  +  P5)  <  4  P2  ^3 

c5  <  4(kTq)c3 


( IV- 146 ) 
( IV-147 ) 
( I V- 148 ) 
( IV-149) 
(IV-150) 


[X4+  ^  +  ^  (^4  +  <  4  (x2  +  ?  ^  +  "3  ^  ( I V-151 ) 

Combining  the  equation  for  resistance  to  diffusion,  (IV-143),  with  the 


definition  of  diffusive  resistance,  ( 1 11-116) ,  and 

ignoring  chemical 

m 

action,  we  obtain 

c  w  .  P<2)  (f  <2>  -  b  (2))  -  t  (2) 

C3  wi  p  Di  ;  ji,j 

( IV- 152 ) 

xmtm 

Neglecting  inertia  effects  and  assuming  v^^  « 

(2) 

v .  ,  the  above 

■J  7 

km  '  • 

equation  reduces  to 

c3  »,(2)-(2)  N,2,-jl,'2) 

( IV- 153 ) 

Wy 

or, 

I 

',,2)-<c3r‘  (p(2)  S,2> 

( IV- 154 ) 

which  is  d’Arcy's  equation  of  fluid  flow.  In  ( I V- 1 52 ) ,  ignoring  body 
forces  b.j(2)  an(j  writing  f^)  as  u.^  , 


This  equation,  for  isotropic  fluid  pressure,  has  the  same  forms  as 
Garg's  [71]  statement  of  Biot's  equation  of  momentum  balance,  allowing 
for  interaction. 

Green  and  Naghdi  [87]  extended  Crochet  and  Naghdi's  [45]  work  to 
anisotropic  solids  and  to  include  initial  stresses.  In  the  absence  of 
chemical  action  and  for  uniform  temperature,  the  Clausius-Duhem 
inequality  to  be  satisfied  is,  for  a  binary  mixture, 


H^Wk%jVVkVkVVkVkS(kS>>W,J 


>  o 


( IV-156) 


where  0  is  the  change  in  temperature. 

The  explicit  form  of  the  dependence  of  stress  in  the  solid  or  the 
fluid  upon  the  strains,  strain  rates  and  history  of  deformation  will 
not  be  discussed  here.  For  illustration,  the  materials  have  been 
assumed  to  be  linear  elastic.  However,  these  relations  may  be 
nonlinear  elastic,  viscous,  viscoelastic,  hypoelastic  or  elasto- 
viscoplastic.  A  large  volume  of  literature  on  the  subject  exists.  As 
examples  of  recent  work  on  a  variety  of  models  we  cite  Desai  [53-55], 
Prevost  [149],  Digby  [59],  Dragson  [60],  Christofferson  [40],  Coon 
[41],  and  Mullenger  [135].  In  developing  plasticity  theories  for 
saturated  porous  media,  Carroll  [35]  has  pointed  out  the  difficulty  in 
defining  plastic  strains.  Shi  [177]  derived  the  explicit  form  of 
stress-strain  relations  following  Crochet  [45]  and  Green  [87].  We 
note,  however,  that  Shi  assumed  the  Helmholtz  free  energy  for  the 
mixture  to  be  the  sum  of  Helmholtz  free  energies  of  the 
constituents.  This  would  imply  uncoupled  constitutive  equations  for 
the  'equilibrium'  stresses  as  contended  by  Sandhu  [156-158].  However, 
this  assumption  could  lead  to  other  difficulties  because  the 
derivation  of  balance  laws  in  Green  and  Naghdi's  theory  was  based  upon 
a  different  set  of  assumptions. 


4.4.  CONSTITUTIVE  EQUATIONS  FOR  POROSITY. 


Gibson  [77]  treated  porosity  as  a  function  of  effective  stress 


and  proposed  compliance  relationships  in  the  form 


n(2)  =  fltjj.Xj.t) 


( I V- 1 5  7 ) 


Walsh  [195]  regarded  the  pore  space  or  the  volume  fraction  of  the 
pores  to  be  a  function  of  the  solid  stress.  Thus,  as  implied  in  (HI- 
39), 


v<2>  .  „<2>v 

Consequently, 

,(2) 

ITT 


dV' 


=  V 


dn 


(2) 


dt . . 
JJ 


dtJJ 


nr 


+  n 


(2)  dV 


dt . . 
JJ 


ITT 


( IV- 1 58 ) 


( IV - 1 59 ) 


Hence 


dn 


(2) 


dt . . 
JJ 


nr 


I  dV 
V 


(2) 


n^  dV 


1  dV 


(2) 


V<2>  dV 


dt . . 
JJ 


m 


dt . . 
JJ 


nr 


vdtjj 


nr 


For 


V<2>  <  < 
dn<2> 


V2, 


dt . . 
JJ 


ITT 


1  dV 

V 
0 


(2) 


dtjj 


ITT 


( IV- 160 ) 


( IV- 161 ) 


Garg  ar d  Nur  [70]  pointed  out  that  Walsh's  analysis  was  acceptable  for 
very  dense  rock  i.e.,  n^'  nearly  equal  to  1.  For  the  general  case, 
they  introduced  bulk  and  intrinsic  solid  compressibilities  and  set  up 
more  general  expressions.  They  proposed  constitutive  equations  for 
porosity  based  on  the  existence  of  an  energy  function. 

Aifantis  [4]  assumed  effects  of  changes  in  fluid  pressure  and  the 
solid  stress  to  be  additive  and  proposed  a  compliance  relationship 


A  n(2)  =  a  A  w  +  b  a  tkk(1)  ( IV-162) 

Morland  [130]  introduced  scaling  functions  for  density  and 
stress.  Hence  constitutive  equations  would  be  needed  for  these. 
These  factors  could,  in  general,  depend  upon  the  deformation  history 
of  all  the  constituents  and  thereby  define  an  interdependence  between 
the  partial  constitutive  relations.  However,  a  simplifying  assumption 
was  made  that  the  area  factors  are  uniquely  related  to  volume  scaling 
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**<  „  ’  •  4 


ft 


factors  i.e.  volume  fractions.  Then  only  the  constitutive  equations 
for  volume  fractions  are  needed.  Further  assuming  that  the 
relationship  is  linear  and  the  same  for  all  constituents,  Morland  set 
the  area  fraction  equal  to  the  volume  fraction,  a  result  postulated 
directly  by  Biot  [21],  Morland  [130]  assumed  the  porosity  n^2)  to  be  a 
function  of  the  current  deformations  of  the  constituents.  Assuming  a 
linear  relationship  (i.e.  using  only  the  first  invariants  of  the 
infinitesimal  strain  tensor) 


n<2>  .  n<2>  |lu,  0)  ♦„«,  <2>) 

J  J  J  J 


( IV- 163 ) 


where  c,  d  are  of  order  unity  or  less.  Equivalently,  the  relationship 
could  be  stated  in  terms  of  partial  stresses  instead  of  strains  as 


n<2>  -n<2>  (H.t, .<»  ,b.) 

J  J 


( I V— 1 64 ) 


To  include  dilatancy,  the  relationship  was  generalized  further  to 

nfZ)  *  n0(2)  (1  +  a..  t^.^U  b  n)  (IV-165) 


Of  course,  temperature  dependence  of  porosity  can  be  included.  Also, 
to  allow  for  the  influence  of  shearing  stresses  in  the  solid  on 

porosity,  Morland  [131]  included  the  second  invariant  of  the  solid 
stress  in  the  equation. 

Noting  the  dependence  of  t-,  n  upon  e--^  Morland  [110], 

J  U  J  >J 

following  Garg  and  Nur  [70],  proposed  a  constitutive  equation  for  the 
ratio  of  n^1^  to  nQ^^  as  an  expansion  in  and  e-^2).  Garg  and 

Nur  [70]  assumed  the  existence  of  an  energy  function  in  terms  of 

porosity  and  the  volumetric  stress  defined  by  (II 1-31)  to  get 

symmetric  expressions  for  hydrostatic  stresses  in  the  fluid  and  the 

solid.  Morland  [133]  extended  this  concept  to  cyclic  loading  using  a 
hypoelastic  model  for  the  deviatoric  response.  The  change  in  porosity 
was  expected  to  have  the  form  [133] 

n(Z)  =  H(n(Z),  K)  =  H^n^)  (K)  >  0  ( IV-166) 


-107- 


.gru+AT A-  mJL ,  ZaLZuOk fawftfc mU  L wtll  ‘ -a  JlAa 


Morland  [133]  suggested  a  polynomial  representation  for  He  also 

described  experiments  to  determine  the  coefficients  in  ( I V- 164 )  and 
( IV-165) . 

4.5.  CONSTITUTIVE  EQUATIONS  FOR  SPECIFIC  VOLUME. 

4.5.1.  Kenyon's  Theory. 

Kenyon  [105]  considered  the  fluid  density  to  be  a  function  of  the 
intrinsic  fluid  pressure  and  the  bulk  volume  measure  J,  (III-10). 


,(2)  . 


=  1  -  n 


(1)  -W*»(2)  (p(2)*.  J) 

p 


( I V- 16  7 ) 


would  be  a  constitutive  equation  for  n^.  For  an  incompressible 
mixture  n^2)  "must  depend  upon  J  as  fluid  must  replace  solid  when  the 
latter  changes  its  macroscopic  volume.  In  addition,  at  constant  J, 
n^2)  will  generally  depend  on  the  amount  of  fluid  present  since  fluid 
can  be  used  to  compress  the  cells  or  granules  independent  of  any  bulk 
deformation  of  the  solid".  In  choosing  the  functional  form  for  (IV- 
167),  Kenyon  assumed  that  the  bulk  strain  deviation  will  not  change 
pore  volume  i  .e.  only  the  bulk  volume  change  of  the  solid  represented 
by  J  can  cause  microscopic  volume  changes. 

Kenyon  [106]  regarded  velocity  of  the  fluid  to  be  a  function  of 
the  fluid  pressure  and  the  velocity  of  the  solid  a  function  of  the 
fluid  pressure  and  the  solid  stress  i.e.. 


-V(<2>  (,*) 


( I V— 1 68 ) 


( IV-169] 


Then,  differentiation  leads  to  the  incremental  relationships 


;,<*>  /  » <2> 
^  3tc  1 


•a)  *  v  (i) +  ;,  o) 


( IV- 170) 


(IV-171) 


Si 


Kenyon  proposed  using  these  incremental  equations  to  identify  the 
functional  relationships  ( IV- 168)  and  ( I V- 169 )  by  experimental 
observations. 


4.6.  COMMENTS. 

For  fluid-saturated  soils,  indentif ication  of  quantities  for 
which  constitutive  equations  are  required  along  with  various 
approaches  to  development  of  constitutive  relationships  has  been 
reviewed.  In  addition  to  the  partial  stresses  in  solid  and  the  fluid, 
constitutive  relations  have  been  proposed  for  diffusive  resistaznce 
and  for  porosity.  Micromechanical  approaches  have  been  concerned  with 
relating  the  macroscopic  behaviour  of  a  saturated  porous  solid  to  the 
intrinsic  behaviour  of  the  solid  matrix  and  the  fluid.  Thermodynamic 
considerations  based  on  the  Clausis-Duhem  inequality  and  the  existence 
of  a  free-energy  function  provide  certain  guidelines.  However,  when 
writing  the  specific  form  of  the  Helmholtz  free  energy  function  for  a 
mixture,  it  is  important  to  keep  in  view  ( 1 1 1 -121 )  which  expresses  the 
relationships  between  the  internal  energies  of  the  constituents  and 
the  internal  energy  of  the  mixutre. 

It  is  clear  form  (II I -42 )  that  for  a  complete  description  of  the 
kinematics,  if  the  solid  grains  are  compressible,  it  is  necessary  to 
introduce  an  additional  quantity.  This  has  led  to  the  constitutive 
equations  for  porosity. 

The  difficulty  with  the  constitutive  modelling  of  saturated  soils 
lies  in  design  of  procedures  for  determination  of  the  actual  form  of 
the  constitutive  relationships.  The  mechanical  quantities  e.g. 
effective  and  partial  stresses,  diffusive  resistance  etc.  need  to  be 
carefully  defined.  Mori  and 's  work  appears  to  offer  the  most  promising 
approach  in  this  direction. 
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SECTION  V 
DISCUSSION 


Theories  of  mechanical  behaviour  of  saturated  soils,  including 
various  approaches  to  description  of  motion,  derivation  of  balance 
laws  and  development  of  constitutive  relationships  have  been 
reviewed.  These  include  early  work  on  soil  consolidation  and  wave 
propagation  in  saturated  soils  as  well  as  modern  approaches  based  on 
theories  of  mixture^. 

Early  work  on  quasi-static  behaviour  of  saturated  soils  was  based 
largely  upon  phenomenological  observations,  'reasonable'  extensions 
and  generalizations  and  certain  ad  hoc  assumptions.  More  recent  work 
has  examined  the  thermodynamic  basis  of  the  mechanics  of  mixtures. 

Linear  dependence  between  nominally  one-dimensional  waterflow  and 
the  fluid  pressure  gradient  observed  by  d’Arcy  was  generalized  to 
three-dimensional  flow,  to  thermoelastic  flow  and  to  multicomponent 
flow  in  presence  of  thermal  and  chemical  potential  gradients.  It  was 
shown  to  be  a  special  case  of  the  more  general  Onsager  relations 
between  fluxes  and  Prigogine  forces.  Similar  constitutive  formula¬ 
tions  were  proposed  considering  relative  motion  in  the  case  of  a 
moving/deforming  solid  matrix. 

Terzaghi  introduced  the  notion  of  effective  stress.  The  dual 
description  he  used  quite  appropriate  to  his  problem  has  been  the 
source  of  considerable  confusion.  He  defined  the  effective  stress  as 
the  stress  related  to  soil  deformation  and  also  as  the  partial  stress 
in  soil.  He  assumed  additivity  of  the  partial  stresses  and  considered 
the  fluid  pressure  to  be  acting  over  100  percent  of  the  area  of  the 
soil.  This  assumption,  combined  with  defining  the  pore  pressure 
measured  by  a  piezometer  as  the  fluid  pressure,  in  the  case  of  zero 
fluid  velocity  leads  to  the  soil  being  'buoyant'  in  water  with  a 
reduced  density  and  corresponding  effective  stress. 


Biot  [15]  while  developing  the  equations  of  three-dimensional 
consolidation,  assumed  the  irreversible  process  of  soil  consolidation 
to  be  a  sequence  of  reversible  processes  and  assumed  existence  of  an 


energy  function  in  terms  of  the  soil  strains  and  the  fluid  content  or 
the  change  in  fluid  density.  Existence  of  a  dissipation  function  was 
assumed  as  well  as  the  additivity  of  partial  stresses.  Gibson's 
theory  [77]  of  one-dimensional  consolidation  considered  a  reference 

volume  of  the  soil  and  fluid  flow  relative  to  it.  This  essentially 
amounts  to  use  of  convected  coordinates  for  the  solid  and  a 
'referential  description'  of  fluid  kinematics. 

For  wave  propagation  in  soils,  two  approaches  have  been  used. 
The  'Engineering  Approach'  introduced  by  Seed  was  initially  based  upon 
the  total  stress  analysis  using  methods  of  structural  dynamics 
assuming  the  saturated  soil  to  be  a  single  material  i.e.  relating  the 
total  stress  to  kinematics  of  the  soil  skeleton.  Later  development 
used  effective  stress  vs.  soil  strain.  Change  in  fluid  pressure  was 
assumed  to  be  related  to  volume  changes  of  the  soil.  Fluid  inertia 
was  ignored  and  fluid  velocity  relative  to  the  soil  was  part  of  the 

process  of  'dissipation'  of  'excess'  fluid  pressure  built  up  due  to 
soil  volume  changes.  The  procedure  did  not  consider  simultaneous 
coupled  motion  of  the  soil  and  the  pore  fluid.  The  theory  apparently 
was  successful  in  predicting  liquefaction  in  layered  soil  systems 
subjected  to  shear  excitations  applied  to  the  base  as  in  seismic 

event.  However,  the  approach  cannot  predict  post-liquefaction 
distribution  of  pore-water  pressures,  does  not  give  reliable 
acceleration  and  displacement  response  histories,  cannot  be  extended 
to  two-  and  three-dimensional  cases  and  its  use  requires  considerable 
experience  and  'judgment'  in  addition  to  extensive  laboratory 
testing.  For  these  reasons,  the  theory  is  not  likely  to  be  useful  for 
application  to  study  of  soil  systems  under  blast  loading. 

Biot's  formulation  of  the  coupled  problem  was  the  basis  of 

theoretical  solutions  by  8iot  [21,22],  Deresiewicz  [52],  Garg  [69]  and 

Chakraborty  [38]  among  others.  Ghaboussi  [73],  Prevost  [147]  and 

Sandhu  [157]  implemented  this  theory  in  finite  element  computer 

programs.  In  this  theory,  in  addition  to  the  ad  hoc  assumptions  on 
which  Biot  had  set  up  the  constitutive  relationships,  he  assumed  the 
kinetic  energy  to  be  a  quadratic  function  in  the  velocities  of  the 

fluid  and  the  solid  including  a  coupling  term  rather  than  simply  the 


sum  of  the  kinetic  energies  of  the  constituents.  A  dissipation 
function,  quadratic  in  relative  velocity  was  introduced.  Garg's 
formulation  of  the  dynamic  problem  is  identical  to  Biot's.  However, 
he  assumed  the  kinetic  energies  of  the  constituents  to  be  additive  and 
related  the  coefficients  appearing  in  the  equations  of  the  motion  to 
be  properties  of  the  constituents  and  the  volume  fractions.  While 
developing  finite  element  solution  procedures  to  Biot's  equations, 
Ghaboussi  [73]  introduced  an  arbitrary  Rayleigh's  type  viscous  damping 
term  to  reflect  the  intrinsic  damping  in  the  saturated  soil  in 
addition  to  the  dissipation  of  relative  motion  based  on  permeability 
and  porosity. 

Theories  of  multicomponent  mixtures  and  their  applications  to 
dynamics  of  saturated  soils  were  reviewed  with  a  view  to  evaluate  the 
'reasonableness'  of  the  various  assumptions  by  reference  to 
fundamental  principles  of  mechanics.  Several  approaches  have  been 
used  to  describe  the  motion  of  constituents  in  a  multicomponent 
mixture.  One  is  to  refer  the  motion  of  each  component  to  its 
reference  configuration.  Another  is  to  refer  the  motion  of  all 
constituents  to  the  reference  configuration  of  one  of  them.  Yet 
another  is  to  refer  all  motion  to  the  current  configuration  which  is 
same  to  all  constitutents.  Superposition  of  relative  diffusive  motion 
of  the  constituents  upon  the  mean  motion  of  the  mixture  as  a  whole  is 
also  used.  For  a  binary  mixture,  most  investigators,  following  Biot, 
describe  the  motion  of  the  solid  with  respect  of  its  reference 

configuration  but  the  motion  of  fluid  is  described  as  relative  to  the 
solid.  Another  procedure,  used  by  Gibson,  is  to  refer  to  a  material 
region  in  the  reference  configuration  consisting  of  a  certain  set  of 

particles  of  one  of  the  constituents  so  that  the  bounding  surface  of 
this  set  varies  with  time.  This  essentially  represents  the  use  of 
convected  coordinate  description  for  the  solid  in  a  binary  fluid-solid 
mixture.  In  most  theories  the  deformation  of  each  constituent  is 

referred  to  an  initial  reference  configuration  but  the  motion  is 
referred  to  the  place  coordinates  i.e.  the  current  configuration.  For 
granular  porous  media,  some  investigators  regard  the  total  deformation 
as  made  up  of  two  parts;  one  related  to  deformation  of  the  solid 
particles  and  the  other  to  their  rearrangement  i.e.  change  in  pore 
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geometry.  These  approaches  are  of  assistance  in  understanding  the 
relationships  between  the  ‘bulk’  and  'intrinsic'  behavior  of 
constituents.  Volume  fraction  theories  for  compressible  materials 
introduce  the  volume  fractions  as  additional  variables. 

The  notion  of  a  mixture  as  a  whole  is  often  introduced.  For  the 
case  of  no  relative  motion  between  the  constituents,  the  mixture  will 
have  motion  and  deformation  as  a  single  body.  However,  if  relative 
motion  is  present,  the  mixture  does  not  satisfy  the  axioms  of 

continuity  and  its  corollary,  the  principle  of  impenetrability. 
Therefore,  the  mixture  cannot  be  regarded  as  a  continuum  in  motion, 
the  particles  of  which  are  assumed  to  be  constructed  by  superposition 
of  constituent  particles,  except  in  the  case  of  no  relative  motion 
between  the  constituents . 

Atkin  [8]  pointed  out  that  the  mixture  having  a  mass  density 
calculated  as  the  weighted  average  of  the  constituents  cannot  be 

associated  with  a  material  in  the  physical  sense.  The  mixture  so 

defined  consists  only  of  centers  of  mass  and  not  a  set  of  particles. 
It  should  be  kept  in  mind  that  the  mixture  density  is  purely 

mathematical  entity  having  no  physical  interpretation  except  in  the 
special  case  of  no  relative  motion.  Many  investigators  introduce  a 
material  rate  for  the  mixture.  This  rate  again  has  no  physical 
meaning  being  executed  not  on  a  particle  but  a  center  of  mass  having 
only  instantaneous  existence. 

Truesdell  postulated  equations  of  balance  of  mass,  linear 
momentum,  moment  of  momentum  and  energy  for  each  constituent  and 
required  that  the  form  of  the  equations  be  the  same  for  each 
constituent  and  for  the  mixture.  Apparently,  the  existence  of  the 
mixture  as  a  continuum  in  motion  was  implied  in  this  line  of 
thought.  To  accomplish  this  identity  of  form,  the  total  stress 
tensor,  the  total  heat  flux  vector  and  the  specific  energy  supply  had 
to  be  specially  defined  and  did  not  equal  the  sum  of  the  corresponding 
quantities  for  the  constituents  as  shown  in  (III -57 ) ,  (  1  1  1-59)  and 
(II 1-61 ) .  The  energy  of  the  mixture,  however,  equalled  the  sum  of  the 
energies  of  the  constituents  (III -55 ) .  The  additivity  of  stresses, 
heat  fluxes  and  specific  energy  supply  could  be  restated  using  special 
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definitions  for  these  quantities  (III -58)  „  (III -60 ) ,  (II 1-62) .  For 
component  velocities  relative  to  the  mean  velocity  of  the  mixture 
approaching  zero,  the  stresses,  heat  fluxes  and  the  specific  energy 
supply  approach  additivity.  On  this  basis,  Bowen  [29]  regarded 

Green's  theory  to  be  a  special  case  of  Truesdell's.  This  is 
incorrect.  As  the  mixture  cannot  be  regarded  as  a  continuum  in 
motion,  Truesdell's  definitions  of  the  total  stress,  heat  flux  and 
specific  energy  supply  are  artificial  and  unnecessary. 

Green  [83]  postulated  the  additive  property  of  stress,  heat  flux 
and  the  energy  supply.  Mori  and  [130]  supported  the  additivity  of 

partial  stresses  on  the  ground  that  tractions  are  additive  and 

Cauchy's  stress  principle  should  hold  for  total  stress  and  total 

traction  as  well  as  for  the  constituents.  Williams  [197]  presented 
similar  arguments.  The  sum  of  internal  energy  densities  of  the 
constituents  did  not  equal  the  internal  energy  density  of  the  mixture 
in  this  theory  (III-121). 

Green  derived  the  balance  laws  form  the  invariance  of  a  rate  of 
energy  equality.  The  form  of  the  balance  equation  is  the  same  as  in 
Truesdell's  theory  but  the  terms  appearing  in  these  equations  have 
quite  different  meaning. 

Some  investigators,  considering  the  special  problem  of  the  flow 
through  deformable  porous  solids,  have  attempted  to  write  the  balance 
equations  in  terms  of  relative  motion  and  porosity,  which  is 
essentially  is  a  measure  of  relative  deformation.  It  would  appear 
that  a  theory  based  upon  balance  equations  written  for  a  reference  set 
of  particles  of  the  porous  soil  but  in  the  current  configuration  would 
be  the  most  appropriate  to  the  problem.  For  the  one-dimensional 
quasi-static  case  this  would  reduce  to  Gibson's  theory  [77]. 

Fukuo  [68]  used  the  equations  of  mass  balance  to  set  up  equations 
in  terms  of  volume  fractions  of  the  constituents  and  specialized  these 
for  incompressible  solid  skeleton  and  also  incompressible  fluid. 
This,  based  on  use  of  d'Arcy's  law  in  conjunction  with  reference  to  a 
fixed  set  of  particles  in  reference  configuration,  may  be  regarded  as 
a  generalization  of  Gibson's  equations.  Hsieh  and  Yew  [95]  set  up  a 
relation  between  rate  of  porosity  change,  rate  of  volumetric  strain 


and  the  relative  velocity  for  a  porous  solid  saturated  with  an 
incompressible  fluid  and  undergoing  small  deformations  in  absence  of 
thermal  and  chemical  effects.  Then  these  quantities,  in  the  case 
considered,  cannot  be  treated  as  independently  variable.  Goodman 
admitted  the  volume  fraction  as  an  additional  quantity  in  the  equation 
of  energy  balance  for  a  porous  material.  Rate  of  work  terms  were 
associated  with  the  rate  of  change  of  porosity  over  volumes  and 
surfaces  using  generalized  forces. 

In  setting  up  constitutive  relationships,  in  addition  to  the 
components  of  total  stress  on  saturated  soil,  the  existence  of  fluid 
pressure  and  the  Interaction  force  between  components  has  been  long 
recognized.  Skempton  [179]  eliminated  the  need  for  a  constitutive 
equation  for  the  fluid  pressure  by  assuming  it  to  be  a  function  of  the 
stress  components  or  the  invariants.  Later  an  explicit  dependence 
upon  soil  strains  was  empirically  introduced.  Biot  assumed  the 
existence  of  an  energy  function  in  terms  of  the  soil  strains  and  fluid 
content/change  in  density.  In  theories  of  mixtures,  the  quantities 
for  which  constitutive  relations  are  required  are  determined  on  the 
basis  of  number  of  variables  vs.  the  number  of  equations  of  kinematics 
and  balance  or  the  quantities  appearing  in  the  rate  of  energy  equality 
or  entropy  inequality  which  are  not  defined  by  kinematics.  The 
equation  of  mass  continuity  for  soils  with  compressible  grains  involve 
porosity.  This  has  been  the  motivation  for  setting  up  constitutive 
equations  for  porosity.  The  existence  of  Helmholtz  free  energy 
function  has  been  often  assumed.  However,  if  the  energies  of  the 
constituents  be  considered  additive,  there  would  be  no  coupling 
between  their  constitutive  relations.  This  is  contrary  to  Biot’s 
assumptions  and  is  in  line  with  Adkins'  earlier  view  [1,2]  and  with 
Morland  [130].  However,  if  Green's  theory  be  followed,  (II 1-121 ) 
would  indicate  the  type  of  coupling.  Crochet  and  Naghdi's  theory  [45] 
would  have  to  be  interpreted  in  that  light.  This  does  appear  to  lead 
to  Biot's  equations.  Mori  and 's  theory  [130]  admitting  coupling 
between  constitutive  equations  for  partial  stresses  in  the  solid  and 
the  fluid  without  assuming  the  existence  of  energy  function  for  the 
mixture  appears  to  be  worth  further  examination. 
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